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tomorrow 


for Progress 


Open ears lead to the open thinking that 
made possible today’s rocket powerplants. 
They are alert ears—ears acutely sensitive 
to the unsolved problems that must be 
conquered in designing powerplants for the 
vehicles of the future. 


RMI puts a premium on just such thinking. 
Its engineers and scientists form a talented, 
alert team, efficiently and effectively 
meeting the rocket power needs of tomorrow. 


Engineers, Scientists— Perhaps you, too, 
can work with America’s first rocket family. 
You'll find the problems challenging, the 
rewards great. 
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her Tough Job 


Above you see an extremely complex body section of a 
missile being finish machined on a horizontal boring mill. 
You are looking at another intricate example of Diversey 
skill and craftsmanship. A visit to our plants would show 
you many comparable jobs that would astonish you. 

At Diversey you have the LARGEST FACILITIES exclu- 
sively devoted to your missile metal machining problems. 
You work with fast, precise, and progressive technical 
people who know what works and what won't. Bring your 
tough jobs to Diversey. . 


HYDROSPINNING 
NOW AVAILABLE 


A new Hydrospinning Division has 
been formed at Diversey which uses 
the latest and largest equipment to 
produce intricate missile parts. 
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Scope of JET PROPULSION 


This Journal is a publication of the Ameri- 
can Rocket Society devoted to the advance- 
ment of the field of jet propulsion through 
the dissemination of original papers disclosing 
new knowledge or new developments. As 
used herein, the term ‘‘jet propulsion” 
embraces all engines that develop thrust by 
rearward discharge of a jet through a nozzle 
or duct; and thus it includes air-consuming 
engines "and underwater systems as well as 
rockets. JeT PROPULSION is open to con- 
tributions dealing not only with propulsion 
but with other aspects of jet-propelled flight, 
such as flight mechanics, guidance, telemeter- 
ing, and research instrumentation. In- 
creasing emphasis will be given to the scientific 
problems of extraterrestrial flight. 


Information for Authors 


Manuscripts must be as brief as the proper 
presentation of the ideas will allow. Ex- 
clusion of dispensable material and concise- 
ness of expression will influence the Editors’ 
acceptance of a manuscript. In terms of 
standard-size double-spaced typed pages, a 
typical maximum length is 22 pages of text 
(including equations), 1 page of references, 
1 page of abstract, and 12 illustrations. 
Fewer illustrations permit more text, and vice 
versa. Greater length will be acceptable 
only in exceptional cases. 

Short manuscripts, not more than one 
quarter of the maximum length stated for full 
articles, may qualify for publication as 
Technical Notes. They may be devoted 
either to new developments requiring prompt 
disclosure or to comments on previously pub- 
lished papers. Such manuscripts are usually 
published within two months of the date of 
receipt. 

Sponsored manuscripts are _ published 
occasionally as an ARS service to the indus- 
try. A manuscript that does not qualify for 
publication according to the above-stated 
requirements as to subject scope or length, 
but which nevertheless deserves widespread 
distribution among jet propulsion engineers, 
may be printed as an extra part of the Journal 
or as a special supplement, if the author or 
his sponsor will reimburse the Society for 
actual publication costs. Estimates are 
available on request. Acknowledgment of 
such financial sponsorship appears as 4 
footnote on the first page of the article. 
Publication is prompt since such papers are 
not in the ordinary backlog. 

Manuscripts must be double spaced on one 
side of paper only with wide margins to allow 
for instructions to printer. Include a 100 to 
200 word abstract. State the authors’ 
positions and affiliations in a footnote on the 
first page. Do not type equations; write 
them in ink. Identify unusual symbols or 
Greek letters for the printer. References are 
to be grouped at the end of the manuscript 
and are to be given as follows: for journal 
articles: authors first, then title, journal, 
volume, year, page numbers; for books: 
authors first, then title, publisher, city, edi- 
tion, and page or chapter numbers. Line 
drawings must be clear and sharp to make 
clear engravings. Use black ink on white 
paper or tracing cloth. Lettering should be 
large enough to be legible after reduction. 
Photographs should be glossy prints, not 
matte or semi-matte. Each illustration must 
have a legend; legends should be listed in 
order on a separate sheet. 

Manuscripts must be accomplished by 
written assurance as to security clearance in 
the event the subject matter lies in a classified 
area or if the paper originates under govern- 
ment sponsorship. Full responsibility rests 
with the author. 

Submit manuscripts in duplicate (original 
plus first carbon, with two sets of illustra- 
tions) to the Editor, Martin Summerfield. 
Professor of Aeronautical Engineering, Prince- 
ton University, Princeton, N. J. Preprints 
of papers presented at ARS national meetings 
are automatically considered for publication. 
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RECENT ADVANCES IN CERMETS’ 


BY ROBERT STEINITZ 


American Electro Metal Div. of Firth Sterling, Inc., Yonkers, N. Y. 


The author, born in 1909 in Germany, 
came to this country after completing uni- 
versity training at Goettingen where he 
obtained a doctorate in physics. Since 
coming to this country in 1936 he has been 
intimately associated with problems in 
solid state physics and metallurgy. As 
Director of Research of the American 
Electro Metal Division of Firth Sterling, 
Inc., Dr. Steinitz has contributed exten- 
sively to the powder metallurgy and high 
temperature material fields. Co-author of 
several books and contributor of about 20 
papers to the scientific literature, Dr. 
Steinitz now focuses his attention on 
cermets and problems associated with 
their practical use. 


Introduction 
Competition with Superalloys 


URING the last two years, the 

somewhat vague and wishful think- 
ing of the past about applications of 
cermets as high temperature materials 
has come to a more definite conclusion 
as to where these materials could and 
should be applied. The shortcomings of 
cermet materials—either based on mix- 
tures of metals with true oxide ceramics 
or on combinations of refractory hard 
metals with so-called binder materials— 
have, first of all, almost excluded appli- 
cations below 1700 to 1800 F. The 
development of special nickel base alloys 
(1)? has stretched the usable range to 
above 1700 F, with long-time strength 
properties close to those of the most 
common cermet material, titanium 
carbide with nickel. It does not seem 


1 Director of Research 
* Numbers in parentheses indicate Ref- 
erences at end of paper. 


worthwhile to struggle with the brittle- 
ness and fabrication difficulties of 
cermet materials, when it is possible to 
sacrifice only a small amount of strength 
and to use instead a material which at 
least can be cast, if not fabricated by 
rolling and other common forming 
procedures. However, the limit of the 
nickel base materials has been reached, 
as the melting point is being closely 
approached. The excellent high tem- 
perature properties of cermets (2) be- 
come more apparent when the tempera- 
ture considered for their application is 
so high that no materials based on iron, 
nickel, or cobalt alloys are usable any 
longer. The development of molyb- 
denum has to a limited extent influenced 
the applicability of cermet materials. 
While molybdenum has phenomenal 
strength at high temperature, its lack 
of oxidation resistance is a serious handi- 
cap. 

At present, the use of cermet mate- 
rials on the basis of carbides, borides, 
nitrides and silicides has, therefore, ex- 
perienced a serious re-evdluation. This 
is not due to any failure of cermet mate- 
rials or to any unfulfilled promise by 
them. Properties of cermets determined 
earlier are still correct and have even 
been improved lately, but the develop- 
ment of nickel base materials has forced 
the cermets into a new field, where they 
can perform, while all other materials 
fail. The use of such extremely high 
temperatures is not yet common prac- 
tice, as a large number of serious diffi- 
culties confront the designer. To be 
able to make certain highly stressed 
parts out of cermet materials for use at 
temperatures above 3000 F, is not suffi- 
cient for the designer who has to take 
the behavior of the complete assembly 
into account. 


Pure Hard Metals 
for Still Higher Temperatures 


The word “cermet” should be some- 
how re-defined in the light of the above. 
If titanium carbide with a nickel binder 
has so far been considered a typical 
representative of this kind of materials, 
it should be understood that the mate- 
rial would have only limited applica- 
bility in the field of very high tempera- 
tures because the melting point of the 
binder, or rather the eutectic tempera- 


ture between the metal and the hard 
phase which is often considerably lower, 
limits the applications. The useful tem- 
perature is even lower than this point. 
Titanium carbide-nickel cermets cannot 
be employed much above 2000 to 2200 
F. Boride base cermets with a chro- 
mium or chromium-molybdenum binder 
have a somewhat higher range of pos- 
sible application, but the gain in tem- 
perature range is not large. 

The range in which these materials 
presently are being tried lies consider- 
ably higher. In such cases, the amount 
of binder is decreased to a very small 
minimum, or the refractory hard metals 
are used by themselves without the 
addition of another metal. If a liquid 
phase during sintering is desired, this 
can in some cases be obtained by the 
eutectic between two or more hard 
metals, or by the addition of a lower 
melting hard metal to the main con- 
stituent. Actually these materials no 
longer deserve the name ‘“cermet’’— 
they contain neither a ceramic nor a 
pure metal—but common usage right 
now is to include these pure refractory 
metallic materials in this category. 

The use of temperatures from 2500 to 
5000 F requires the highest melting ma- 
terials known (3). In addition to the 
refractory metals, such as molybdenum, 
only two groups of material can be con- 
sidered at all. The first is that of the 
ceramic oxides, where melting points 
around the above mentioned tempera- 
tures are quite common. The other 
group is that of the refractory hard 
metals, as a number of carbides, borides 
and silicides melt within this range. 
Applications include rocket motors and 
combustion chambers in high speed 
planes. In addition, the problem of 
aerodynamic heating in supersonic flight 
has opened a new field for these mate- 
rials. At present, cooling of the exposed 
parts is considered the easier method. 
New production procedures had to be 
devised to make it possible to apply the 
high melting compounds to such parts. 
These manufacturing procedures or 
coatings will be discussed later. An- 
other field of application which has just 
started to open its doors lies in mate- 
rials for nuclear reactors (4). Some of 
the hard metal materials have desirable 
nuclear properties for certain applica- 


(Continued on page 68) 


* This paper is the first of a monthly series of invited articles that will describe new developments in specialized areas of jet and 


rocket engineering. 
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this base uses Mitchell cameras for recording test operations. 


Tracking at 40,000 feet, this film frame from a Mitchell 
camera operating at 96 frames per second provides an 
accurate record. 


High speed 
Mitchell 
Camera in 
operation on 
tracking 
telescope 
mount during 
test run at 
Salton Sea. 


At the Atomic Energy Commission‘s Salton Sea Test Base this special tracking mount uses a 35mm high speed Mitchell camera. Operated by Sandia Corporation, 


Salton Sea Test Base Uses Mitchell Cameras 
to Capture High Speed Action of Dummy Bombs 


High speed flight and laboratory tests, which hitherto have been difficult or 
impossible to view with the human eye, are today providing revealing informa- 
tion through high speed film recordings. 


Typical example of the widespread use of high speed cameras is the Salton 
Sea Test Base in Southern California, where drop testing of “dummy” bombs is 
a major activity. In testing carried on there, by Sandia Corporation for the Atomic 
Energy Commission, as many as 20 Mitchell high speed cameras may record dif- 
ferent angles in the flight of an experimental weapon “shape” from drop aircraft 
to impact area. 


Operating at 48 to 100 frames per second, the Mitchell cameras film accurate, 
steady images with maximum uniformity—even under difficult and complicated 
filming conditions. 


Mitchell cameras play a growing role in today’s research and development— 
just as 16mm and 35mm Mitchell cameras have become the leading professional 
motion picture equipment used by industry, television, and film studios through- 
out the world. Write on your letterhead for further information on the uses of 
Mitchell cameras in the field of military and industrial research. 


Photogram- 
metric mounts 
for 8 Mitchell 
cameras 
determine the 
position in 
space of a 
nuclear 
explosion. 
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666 WEST HARVARD STREET 
GLENDALE 4.CALIFORNIA 
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Inertial Guidance for Rocket-Propelled Missiles 


WILLIAM T. RUSSELL’? 
The Ramo-Wooldridge Corp., Los Angeles, Calif. 


The basic elementsof inertial guidance systems primarily 
intended for rocket-propelled missiles are explained. Pre- 
cision gyroscopes, accelerometers and stabilized platforms 
are described. Guidance schemes appropriate to ballistic 
missiles are discussed and a trajectory for the German V-2 
rocket is presented as an example. The problem of gravity 
computation is considered and the effect of instrument 
errors is explored. A missile steering computer and a 
rocket motor shut-off computer are given for the V-2 
rocket which is used as an example. Numerical values of 
coefficients are calculated for the sample trajectory. 


Introduction 


REVIOUS published papers on inertial guidance have 
considered the problem of navigating at constant altitude 
over the earth’s surface in aircraft, ships or submarines 
(1-5). Although one of the first applications of inertial 
guidance was in the German V-2 missile, the use of inertial 
systems for rocket-propelled missiles has apparently had 
little discussion. In such missiles, great altitudes may be 
reached at a rapid rate, making necessary guidance systems 
which differ considerably from those for constant altitude 
use. Since most rocket-propelled missiles are weapons, 
details of specific inertial systems or the performance of the 
latest instruments are for the most part classified security 
information. However, the basic ideas and principles of 
inertial guidance form a rather simple and obvious applica- 
tion of classical mechanics, and are no more classified than 
are the principles of rocket propulsion or of radar. In this 
paper, these basic ideas are presented. The operation of 
typical inertial components is explained, together with some 
discussion of systems problems. For illust-ative purposes, 
a simple guidance system is conceived for the V-2 rocket. 
Inertial navigation is a form of navigation akin to “dead 
reckoning.’’ The position, velocity, time and orientation 
of a missile are known at the start of a flight. Subsequent 
velocity and positién is then obtained by measurements made 
solely within the missile. This fact eliminates all dependence 
upon clear weather for star-sighting or upon any form of 
electromagnetic radiation which is subject to interference 
or possibly to line of sight limitations. On the other hand, 
inertial navigation is subject to the disadvantage of dead 
reckoning, that is, the longer the duration of the flight, the 
greater are the navigation uncertainties. 


Basic Elements 


The basic principle is extremely simple. The missile ac- 
ccleration relative to a known reference frame established by 
an initial orientation is measured inside the missile. This 
acceleration is then integrated to give velocity and position. 
A guidance or steering function may be added to the naviga- 
tion or position determination function by comparing present 
velocity and position with stored values defining the desired 


' Presented at the ARS Semi-Annual Meeting, San Francisco, 
Calif., June 10-13, 1957. 

* Manager, Inertial Guidance Dept., Space Technology Labora- 
tories. 

* Numbers in parentheses indicate References at end of paper. 
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flight. After suitable computation, steering and thrust 
signals can be sent to the missile control system. The basic 
conceptual elements of all inertial systems can be identified 
as an attitude reference, an accelerometer, a clock, a computer, 
and a memory. 

The attitude reference is a device which is oriented at the 
start of the flight and then maintains or remembers this 
orientation. In effect, it remembers “which way is up,” 
and which way is North. Philosophically, there are many 
ways that this can be done, but most practical schemes 
employ gyroscopes physically mounted on a stable platform. 
Most accelerometers are single-degree-of-freedom instruments, 
that is, they are only sensitive to one component of accelera- 
tion. For rocket-propelled missiles, accelerations may be 
large in all three space coordinates and, in general, three 
accelerometers, one for each coordinate axis, must be used. 
For submarine, ship or aircraft in essentially constant altitude 
flight, the accelerometers are usually oriented in the local 
horizontal plane and only two are required, the vertical 
component being monitored by an altimeter. The computer 
and clock are needed for the process of integration to obtain 
velocity and position, for computation of gravity, for cal- 
culation of the earth’s rotation and for computation of steering 
commands. The memory unit stores the vehicle’s prescribed 
path or destination. These five basic elements are con- 
ceptually distinct but each element is not usually a separate 
piece of hardware. For instance, accelerometers may per- 
form integration and contain the memory in addition to 
measuring acceleration. 

Associated with the inertial guidance system is equipment 
at the missile launch point which is used to orient the attitude 
reference, to establish initial conditions and to store informa- 
tion in the memory unit. Closely related to the inertial 
system, but by definition excluded from it conceptually, is 
the vehicle control system which physically steers and sta- 
bilizes the missile, and speeds it up or slows it down. 

In the traditional form of shipboard navigation, the attitude 
reference is given by the position of the sun, moon and stars. 
The accelerometer becomes a bubble level or the horizon. 
The clock is a chronometer. The computer is an almanac 
and a pencil and paper. The memory is written in ship’s 
orders or contained in the brain of the captain. The control 
system is the wheel and rudder and the throttles on the en- 


gines. 


Gyroscopes 


The heart of an inertial guidance system is the set of gyro- 
scopes and accelerometers. It is only because of recent major 
improvements in the precision of these instruments that in 
ertial guidance has become feasible. A leading group in this 
field is the Instrumentation Laboratory of the Massachusetts 
Institute of Technology, directed by Dr. C. 8. Draper. In 
particular, he developed the type of single-degree-of-freedom 
gyroscopes known as HIG gyros (Hermetic Integrating 
Gyros). Other types of gyros (6) can be used for inertial 
systems, but precision HIG types are available commercially 
and therefore will be described and used as examples here. 
A detailed description with many references is given in (7). 

Fig. 1 is a schematic drawing of the HIG type, single- 
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Fig. 1 HIG gyroscope 


degree-of-freedom gyroscope originally developed by the 
Instrumentation Laboratory, MIT. Various commercial 
versions differ in mechanical detail and in type of signal 
generator and torque generator. Basically, the gyro consists 
of a spinning wheel driven by an electric motor, mounted on 
preloaded ball bearings and contained in a _ hermetically 
sealed can or float with shaft extensions. The float is com- 
pletely submerged in a viscous fluid which has the same 
average density as the float and shaft. In other words, 
the float is in neutral buoyancy and no radial forces are car- 
ried by the jeweled pivots at either end of the shaft. This 
serves to reduce friction about the axis defined by the pivots 
to a negligible value and makes the instrument a rugged 
design. Coaxial with the shaft is a signal generator which 
gives a voltage proportional to the angular displacement of 
the float relative to the case, and a torque generator which 
can be used to apply torque to the float. Means may be 
provided for temperature control, for balancing the float 
and for conducting power to the spin motor. 


The basic principle of operation of the gyro can be ex- 
plained in terms of the three axes shown in Fig. 1. The 
spin reference axis lies along the angular momentum vector 
(spin axis) of the wheel when the signal generator output is 
zero. The output axis is normal to the spin reference axis, 
and is the axis about which the float is free to turn. The 
input axis is normal to the output axis and spin reference axis. 
The input quantity to the instrument is an angular motion 
of the case, relative to inertial space, about the input axis. 
The resulting output is a movement of the float relative to 
the case which results in a voltage from the signal generator. 
The operation is explained by the familiar physical fact that 
when a torque is applied to a spinning wheel so as to change 
the direction of its spin axis, the spin axis tends to align 
itself with the torque vector. Conversely, when the axis 
of a spinning wheel is forcibly precessed or rotated, the wheel 
through its bearings exerts a torque about an axis perpen- 
dicular to the axis of forced rotation. In the HIG gyro, 
movement of the case about the input axis causes a forced 
precession of the gyro wheel about this axis. The gyro 
wheel thus exerts a torque on the float about the output 
axis. Initially, this torque accelerates the float, but as the 
float gains angular velocity, the viscous shear torque reduces 
the acceleration to zero and the float reaches a steady angular 
velocity. The gyro torque is then balanced by the viscous 
shear torque and the output angular rate is proportional to 
the input angular rate. 


Bo, = Ho; 
where 


@, = output angular rate 
@®; = input angular rate 
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= fluid shear coefficient 
H = angular momentum 


and for small angles 
H H 


The fact that the output signal is proportional to the integral 
of the input angular rate is the source of the term “‘inte- 
grating gyro.” The gyro thus serves as an attitude reference. 
Once the wheel is spinning, and provided no signals are ap- 
plied to the torque generator, an angular movement about 
the input axis results in an output signal. Over its operating 
range, the instrument is linear, so that a reverse movement 
about the input axis returns the output signal to a null. 

Torques other than the gyroscopic element, viscous dray 
and float inertia torques can act about the output axis. They 
may arise from two sources, intentionally applied through the 
torque generator, and unintentionally applied by various 
disturbances. These torques result in output signals which 
are indistinguishable from those caused by input angular 
rates. They therefore act to change the reference orientation 
of the gyro at a rate proportional to the torque. If the 
torques are uncertainties caused by such things as float 
unbalance, signal generator reaction, fluid convection currents, 
etc., the resulting output signal looks like an input angular 
rate which is called the drift rate of the gyro. This drift 
rate is the basic performance figure of merit for inertial 
navigational use. The lower the drift rate, the better the 
attitude reference is maintained, and the more accurate the 
guidance system can be. 

HIG gyros act as precision angular motion sensors, rather 
than sources of torques to overcome friction or unbalances. 
Normally, they have an operating range of only a few degrees 
and to prevent various cross coupling errors, input angles 
should be kept small. Therefore, for inertial guidance use, 
gyros are usually mounted on a platform or base which is 
servo-driven to maintain the gyro output at or near a null. 
The platform thus remains fixed in orientation relative to 
inertial space, or is rotated at a rate determined by torque 
generator input. 


Accelerometers 


The accelerometer is probably the most critical instrument, 
in so far as accuracy is concerned, for inertial guidance sys- 
tems for rocket-propelled missiles. This is in contrast to 
the constant altitude case in which the gyros are usually the 
most critical item. For applications considered here, the 
accelerometers must work over a wide dynamic range with 
great precision and usually in the presence of severe environ- 
ments. Fig. 2 is a schematic of a class of accelerometers 
suitable for inertial systems (8). 

The instrument operates by applying force to a mass 
inside the instrument case so that the mass moves with the 
case. The mass is termed the seismic or proof mass and 


CASE PICKOFF 
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Fig.2 Basic accelerometer 
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except for instrument dynamics, its acceleration is equal to 
the acceleration of the case. The force which is needed to 
accelerate or move the seismic mass is measured in terms 
of some physical quantity which is proportional to the force. 
For instance, if the force generator is a simple spring, the 
force is measured by the deflection of the spring. A common 
foree generator is obtained by applying current to a coil 
in a magnetic field. As shown in the figure, the current is 
obtained by feedback from a pickoff between the seismic mass 
and the case. The basic equation for the instrument is 


= g 
M 
where 
ad» = acceleration of seismic mass along input axis relative 
to inertial space 
Ff = all nongravitational forces acting on seismic mass along 
input axis 
MM = mass of seismic mass 
¢ = component of gravity along input axis 


This equation clearly shows that an accelerometer does not 
measure acceleration itself, but rather the difference between 
acceleration and gravity. This difference which is equal to 
the quantity #/M has been called thrust acceleration a,. 
This nomenclature stems from the fact that F/M for the 
accelerometer is equal to the corresponding quantity for the 
complete missile, that is, to the sum of the nongravitational 
forces acting on the missile, divided by the instantaneous 
mass of the missile. Most of the nongravitational force 
acting on a rocket-propelled missile is thrust, hence, the term, 
thrust acceleration. To obtain acceleration itself, gravity, 
which depends upon position, must be calculated in the 
computer and added to thrust acceleration. It should be 
noted that when the acceleration itself is accurately known, 
as it is on the earth’s surface, a calibrated special acceler- 
ometer can be used to measure gravity, or conversely, gravity 
can be used to calibrate an accelerometer. It is sometimes 
convenient to take the viewpoint that an accelerometer 
measures acceleration relative to a freely falling coordinate 
frame, that is, to a frame accelerating at the value of g ex- 
isting at the instrument. It is important, of course, that all 
the forces acting on the mass along the input axis be measured, 
and that only these forces be measured. An example of the 
former type is friction between the seismic mass and the case. 
An example of the latter is change in geometry caused by 
distortion so that the instrument becomes sensitive to cross 
acceleration. The above discussion assumes that the ac- 
celerometer is not rotating. If rotation is present, various 
correction terms must be included. 

Fig. 3 is a schematic of a practical type of instrument based 
on the HIG gyro construction. The seismic mass exists in the 
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Fig.3 Constrained pendulum accelerometer 
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Fig.4 Pendulous gyro accelerometer (PGA) 


form of a pendulum and the force generator is the torque 
generator. The pickoff is the signal generator. Flotation 
virtually eliminates uncertainty friction torques at the pivots. 
Because of the pendulous mass, acceleration of the instrument 
along the input axis creates a torque about the float pivot axis. 
This torque causes rotation of the float and a consequent sig- 
nal generator voltage proportional to @,. This voltage is 
used to generate a current which is applied to the torque gen- 


-erator to give a torque which “constrains’’ the pendulum and 


keeps 9, small. The current, J, is thus proportional to ac- 
celeration along the input axis. The gain of the feedback sys- 
tem must be kept quite high so that deflection of the pendu- 
lum under high input acceleration is small. Otherwise, a 
“cross talk’’ torque is developed which is proportional to the 
product of the acceleration along the pendulous reference 
axis and the sine of the deflection angle 6,._ Instruments of this 
type, called force feedback pendulums or constrained pendu- 
lums, are available commercially. They are made with a 
wide variety of dynamic range and frequency response. 

Velocity and position are the quantities of interest in navi- 
gation, rather than acceleration. Increased accuracy and 
reliability may sometimes be obtained by performing integra- 
tion in the accelerometer. Basically, this is done by making 
the force acting on the seismic mass proportional to a rate of 
some kind. For instance, the current fed to the torque gen- 
erator in Fig. 3 might be applied in pulses of constant area 
but variable rate. Pulse rate is then proportional to accelera- 
tion and total number of pulses to velocity. Another useful 
device and one used by the Germans in the V-2 is the pen- 
dulous gyro accelerometer (PGA). This instrument is a 
pendulum in which the force generator is a gyroscopic ele- 
ment. Fig. 4 is a schematic of such an instrument based on 
the HIG gyro construction. The output of the signal genera- 
tor is fed to the servo motor which rotates the gyro case about 
the input axis at a rate such that the torque developed by the 
gyro element just equals the pendulous torque. Omitting dy- 
namics, the torque balance about the float axis is 


H6; = mla 


where 


ml = pendulosity 
a = acceleration along input axis 


hence 


ml ml 
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Fig.5 Stabilized platform 


The angular rate of the gyro case is thus proportional to ac- 
celeration and the total angle turned by the gyro case is pro- 
portional to velocity. As with the instrument of Fig. 3, the 
gain of the feedback loop must be kept high so that 6, is very 
small. The torque generator of the gyro can be used to apply 
additional torques to the gyro float which add to the pendu- 
lous torque. In this way, gravity may be added to the thrust 
acceleration to give true acceleration and velocity. 


System Mechanization 


The two basic instruments of inertial guidance systems have 
been briefly described; namely, gyros to maintain an atti- 
tude reference and accelerometers to measure thrust acceler- 
ation. These instruments must be related so that the accel- 
erometer orientation relative to the attitude reference is 
known. A practical method for doing this consists of mount- 
ing three gyros on a frame so that their input axes are orthog- 
onal (9). This frame, usually called a stable platform, is 
mounted in a gimbal structure attached to the vehicle. The 
gimbal structure serves as a ball and socket joint and allows 
the vehicle to rotate in any way about the platform. Output 
signals from the gyros are fed to torque motors on the gimbal 
axes which serve to overcome gimbal bearing friction and 
keep the gyro outputs at, or close to, a null value. In doing 
so, the platform upon which the gyros are mounted maintains 
its orientation relative to inertial space irrespective of vehicle 
attitude. Three single-degree-of-freedom accelerometers can 
be mounted on the platform with their input axes orthogonal 
but not necessarily parallel to the gyro axes. The accelerom- 
eters are thus nonrotating and their input axes define the 
coordinate system fixed in the attitude reference frame. 

For full angular freedom, four gimbals are required and slip 
rings on the gimbal axes are usually necessary to carry the 
various signal and power leads. For vehicles, such as the 

’-2 rocket, which have limited angular motions, three gimbals 
are sufficient and flexible cabling from the vehicle to the stable 
platform can be used. 

A schematic drawing of a gimbal system with gyros and 
accelerometers (sometimes called an inertial measurement 
unit—IMU) is shown in Fig. 5. Any small angular motion 
of the platform relative to inertial space is sensed by the set of 
gyros. The gyro output signals are used to generate currents 
which are applied to torque motors mounted at the gimbal 
bearings. These rotate the gimbals and platform so as to 
return the platform gyros to a null position. In the figure, 
the output of the z axis gyro is always applied to the pitch 
axis gimbal torque motor. For pitch angles other than that 
shown in the figure, the « and y axes are not parallel to the 
gimbal yaw and gimbal roll axes. Therefore, the x and y 
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gyro signals must be resolved through a synchro resolver on 
the gimbal pitch axis to give proper signals to the yaw and roll 
gimbal torque motors. There is room for much mechanical 
ingenuity in the design of gimbals and stable platforms to 
achieve minimum weight and size and maximum stiffness and 
environmental protection. 

Some means must be provided for orienting the platform 
at the start of a flight. One possible way is simply to set the 
three gimbal angles to the required value. This then orients 
the platform relative to the missile frame. Another more 
accurate way employs the navigation accelerometers or two 
separate pendulums to position the platform relative to the | 
launch point horizontal. Accurate azimuth orientation can 
be done by optically sighting on a mirror mounted on the 
platform from a theodolite external to the missile. 

Computers for inertial guidance systems can be either 
digital or analog. The former, of course, can be made as 
accurate as desired at the expense of added weight and com- 
plexity. The choice of type of computer and method of 
computer mechanization depends upon over-all system char- 
acteristics such as range, accuracy and rocket thrust char- 
acteristics; upon the computing task to be done (called the 
guidance scheme); and upon available computer hardware. 
Much can be done in making approximations and in combining 
computing functions to minimize weight and errors and to 
maximize reliability and flexibility. 


Guidance Schemes 


Rather than discuss the broad subject of guidance schemes 
in general terms, a specific example will be considered. Fig. 
6 shows a trajectory for the V-2 missile which is based upon 
data given in (10). 

Fig. 7 details the trajectory during the powered _por- 
tion of the flight. The missile, under standard conditions. 
burns for 70 see and achieves a velocity of 5860 fps. It then 
follows a ballistic or free-fall path until it impacts on the earth 
at about 230 miles range. For simplicity in making the cal- 
culations, it is assumed that the earth is flat and the gravity 
field constant. Vertical flight of the rocket lasts 8.5 sec 
at which tirne it is assumed the rocket instantaneously ac- 
quires an attitude of 10 deg from the vertical. From this 
time until motor shutoff, the missile follows a gravity turn, 
i.e., the missile is steered so that the thrust vector is paralle! 
to the velocity vector. In making the calculations the aero- 
dynamic forces were neglected; in any case, they are very 
small compared to thrust and gravity. 

There are many ways in which this missile could be guided 
to impact at its intended target. The simplest scheme con- 
sists of aiming the missile at the target and firing it as a gun 
is fired. Guidance is then done only while the missile is con- 
strained by the launcher. A scheme which is closer to the 
usual meaning of the word, guidance, consists of installing 
a programmer in the missile which generates a missile atti- 
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Fig.6 V-2 trajectory 
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Fig.7 V-2 powered flight trajectory 


tude angle as a function of time. This command signal is 
then compared with the actual missile attitude as measured 
by the pitch gimbal angle of an inertial platform such as that 
shown in Fig. 5. The difference between the command angle 
and the actual angle is then used as an input to the missile 
control system to steer the missile in pitch. Yaw and roll 
angles are similarly held to programmed values which are 
usually zero. A clock could be used to shut the rocket engine 
off at 70 see or, alternatively, enough fuel could be placed in 
the missile so that it would burn out at 70sec. This scheme is 
essentially an extension of the previous one in that the launcher 
or gun barrel is greatly lengthened. While simple, guid- 
ance probably would not be accurate because of uncertainties 
in missile parameters. For instance, if the rocket engine 
thrust was, say 10 per cent more than the preplanned or stand- 
ard value, the missile would accelerate more rapidly than 
standard and the velocity at 70 see would be much greater 
than standard. Furthermore, the missile would be closer to 
the target than standard at motor shutoff. These effects 
would cause an impact far beyond the target. Conversely, 
a rocket engine thrust less than standard would result in an 
impact short of the target. Other nonstandard flight pa- 
rameters, such as winds, would also cause target miss. 

The guidance scheme actually used by the Germans in the 
V-2 (11) had an important refinement over the simple scheme 
just described. single integrating accelerometer was 
mounted with its input axis along the missile longitudinal 
axis, and the motor was shut off when the accelerometer output 
(thrust velocity) equaled the preset value. For the example 
shown in Figs. 6 and 7, this would be 5860 fps plus the time 
integral of the gravity component along the trajectory. The 
\-2 guidance system essentially measured only one component 
of velocity. With such a scheme, it is possible to have con- 
siderable target miss because the other two components of 
velocity at motor shutoff may be nonstandard. A further 
improvement in guidance can be made by using a stable plat- 
form and three accelerometers similar to that shown in Fig. 
5. In this case, all three components of velocity are obtained. 

The guidance schemes just discussed are all of the type 
which might be labeled gun barrel guidance. They all guide 
the missile during powered flight and terminate rocket thrust 
so that the missile follows a ballistic path to the target. After 
motor shutoff, the only nongravitational forces acting on the 
missile are aerodynamic. These can be calculated prior to 
flight for standard conditions and their effects included in 
the instructions placed in the system memory. Some target 
miss will be caused by such uncertainties as winds acting on 
the missile after motor shutoff. The inertial system could be 
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used to measure and correct for these uncertainties, but they 
might be smaller than the inertial instrument errors. Of 
course, in this case, guidance after motor shutoff would be a 
mistake. If no atmosphere whatsoever existed, the missile 
would be in a free-fall state after motor shutoff and the 
accelerometers would have zero output. The airborne com- 
puter would calculate gravity, velocity and position, based 
only on the initial values existing at the time of motor shutoff. 

For a nonrotating, spherical earth, the missile trajectory 
lies in a plane. The actual earth does not have a perfectly 
spherical gravitational field but the variations are negligible 
for a missile such as our example and will not be considered 
further here. For inertial guidance applications, it is conven- 
ient to visualize trajectories in terms of inertial or nonrotat- 
ing coordinates. From this viewpoint, the target and launch 
point are moving eastward over the surface of the earth with 
a velocity 


V, = cos 


where R is the earth’s radius, w, is earth angular rate and \ 
is the latitude angle. The missile at launch then has an 
initial horizontal velocity due to the motion of the launch 
point. If the launch point and target both lie on the equator, 
the trajectory will lie in the equatorial plane. If the missile 
is fired over an earth’s pole, the target is moving normal to the 
trajectory plane and hence the missile must be ‘“aimed’’ at 
the location of the target at the predicted time of impact. 
Furthermore, the missile initial velocity is normal to the 
trajectory plane. Therefore, the missile must be aimed 
slightly to the West of the free-fall trajectory plane so that the 
velocity vector at motor shutoff will be in the trajectory plane. 
Intermediate latitudes and azimuths, of course, have both 
lateral and range effects from the earth’s rotation. 


Gravity 


If every missile flight were exactly standard or as precal- 
culated, there would be no reason to calculate gravity during 
the flight. Instead, the effect of gravity could be precal- 
culated and guidance could be done in terms of thrust acceler- 
ation alone. However, nonstandard trajectories will result 
in gravity accelerations different from standard. To 
account for this, either the actual value of gravity can be 
‘alculated in flight or some correction to the standard pre- 
planned value can be made. The manner in which gravity 
is handled depends upon the range, accuracy desired, trajec- 
tory used, and the method of mechanization. 

The components of gravity are a nonlinear function of 
position, and are somewhat difficult to caleulate exactly 
(12-14). Therefore, if possible, some approximation should 
be used. Of course, the simplest approximation considers 
gravity a constant in direction and magnitude. A step more 
sophisticated is to assume a linear variation of gravity with 
displacement from some reference point. Consider, for 
example, a nonrotating coordinate system with origin near 
the midpoint of the powered flight trajectory, with y axis 
vertical and with x axis in the trajectory plane. Then linear 
approximations to the gravity components are 


290 
(Ro + y)? 


where Ry is the distance from the earth’s center to the origin 
and go is the magnitude of gravity at the origin. 


Jo 

2% 
Ry 

These are the first order terms of a Taylor’s series representa- 

tion of the gravity components. Higher order terms could be 

included if necessary. 
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Fig.8 Gravity computer 


As stated earlier, the basic vector acceleration equation in 
coordinate form is 


Q, = Az + Gz a; = acceleration 
dy = Ay t+ Gy a4 = thrust acceleration 
a, = az + gz = gravity 


Block diagrams for the computation of acceleration assuming 
linear gravity terms are shown in Fig. 8. No computation for 
z or lateral direction need be made for a missile such as our 
example of Fig. 6 because no appreciable z position coordinate 
develops, i.e., g, can be considered a constant. 


Improved System 


The improved V-2 guidance system mentioned earlier will 
now be considered in more detail for the trajectory of Figs. 6 
and 7. For simplicity in explanation, it will be assumed that 
the earth is spherical and nonrotating. A complete inertial 
measurement unit as shown in Fig. 5 will be used. The plat- 
form and the missile are oriented at the launch point so that 
the platform z axis is normal to the trajectory plane, the y 
axis is vertical, and the x axis is horizontal and lies in the tra- 
jectory plane. The trajectory plane itself passes through 
the earth’s center, the launch point and the target. Except 
for gyro drift, the platform will maintain its initial orientation 
during flight. The missile is roll-stabilized so that the missile 
pitch axis is parallel to the platform z or pitch axis. By 
means of a programmer, the missile pitch attitude will be 
made a function of time. This is identical to the system actu- 
ally used in the V-2. Any movement of the missile out of 
the trajectory plane caused by winds or thrust misalignment 
will give a signal from the z or lateral accelerometer. This 
signal will be integrated and sent to the missile control sys- 
tem which will steer the missile so that the lateral velocity 
remains essentially zero. A block diagram showing the essen- 
tials of this simple steering scheme is given in Fig. 9. In 
practice, provision for adequate dynamic stability and response 
are necessary. 

The problem now remains as to how to obtain a motor shut- 
off signal so that the miss at the target in the plane of the 
trajectory, called the range miss, M,, will be zero. What is 
done is to precalculate the effect on target miss of small dis- 
placements in the missile burnout position and velocity from 
the standard values. The total miss is then the sum of the 
separate effects. In other words 


M, = C,Ax + C,Ay + C;A% + CjAy 
where the C’s are constants and the A quantities are the differ- 
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ence between the actual and standard values. The above 
equation is simply the linear approximation to a Taylor’s 
series expansion of range as a function of position and ve- 
locity about the standard motor shutoff point. 


oR oR 
R — R, = M, = — (¢ — 2.) + — (y — ys) + 
or oY 


dit 

where FR designates range, the s subscript designates stand- 
ard conditions and the derivatives are evaluated at the 
standard shutoff point. For the example of Figs. 6 and 7, 
the values of the coefficient are 


ft ft 
C,=1- Cz; = 300 — 
ft fps 
C = 29 ft 
fps 


This means, for example, that if the « component of velocity 
is 1 fps greater than standard, then the range will be 300 ft 
greater than standard. 

A block diagram of a computer which continuously cal- 
culates target miss is shown in Fig. 10. The rocket motor 
thrust is terminated when the value of miss, .W,, becomes 
zero. 

Note that gravity computers similar to Fig. 8 are included. 
The location of the target is defined by the initial values, 
Le, Ye; Zs, Ye and by the values of the C’s. These values are 
constants which are set into the computer prior to launch. 

The actual nonstandard missile will have motor shutoff at 
a point near the standard point. For instance, if thrust is 
low, the missile position and velocity will not be as high as 
standard at a given pitch angle (because pitch attitude angle 
depends upon time only). The missile trajectory will there- 
fore lie below the standard trajectory and motor shutoff will 
probably occur at a greater x value than standard. The 
actual trajectory and shutoff point can, of course, be calculated 
in the same manner as was done for the standard case. 


System Errors 


The magnitude of target miss due to various approxima- 
tions and instrument errors can be calculated quite easily, 
at least in principle. The two principal approximations in 
our example are found in linearizing the gravity computation, 
and in linearizing the caleulation of target miss. Errors due 
to these sources will be very small provided that the rocket 
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Fig. 10 Motor shutoff computer 


thrust and pitch program are reasonably close to standard 
values. In fact, it is quite possible that the gravity computa- 
tion itself could be eliminated. 

To consider the effects of accelerometer errors and gyro 
drifts, an important characteristic of the basic acceleration 
equation must first be explored. A coordinate system iden- 
tical to that used for Fig. 8 will be considered. The accelera- 
tion equations are then 


i= + 
where #, 7 are alternate symbols for a, and a, 
= dy + Oy 


Now a small change or perturbation in the value of # is equal 
to a perturbation or error in the value of thrust acceleration 
plus an error in the computed value of gravity. The gravity 
error is composed of any computer mechanization errors and 
error due to using an erroneous position in calculating gravity. 
In what follows, only the latter will be considered. We thus 
have 


Az = + Ag, 


Ay Ag, 

By taking a differential of the gravity expression, we obtain 
expressions for the dependence of gravity error on position 
error. 


Jo 
Ag, = — — 
Ry 

290 
Ag, = +—A 
Jy 


Substituting these expressions into the above leads to the 
differential equations 


Aé + = Aa 


290 
Aj - Ro Ay = Aay 
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The solutions for constant values of thrust acceleration errors 
(or perturbations) are 


A Da, 
cosh —1 
290/Ro Ro 


For positions near the surface of the earth, the numerical 
value of ~/ Ro/go is about 84 min. 

These linearized error equations illustrate a basic limitation 
of pure inertial guidance for long times of flight. That is, 
errors in the vertical direction increase exponentially with 
time. Errors in the horizontal direction, however, are oscil- 
latory. This is one reason that inertial systems for aircraft 
and ships employ altimetérs to measure altitude. For powered 
flight duration of a few minutes, such as that for our example, 
the equations can be closely approximated by the simple 
expressions. 


9 


Az 


Ay Ay = Adtyt 


which are those that would be obtained by ignoring the feed- 
back error from the gravity calculation. 

For our example, the x accelerometer, for instance, must 
have a dynamic range of about 0 to 6 g’s (see Fig. 7). A 
zero offset in this accelerometer of 0.006 g’s (0.1%) will 
produce a velocity error at motor shutoff of about Az = 
(0.006) (32.2) (70) = 13.5 fps. The miss coefficient for + 
was 300 ft target miss per fps variation in x velocity. Hence, 
the range miss at the target caused by this accelerometer 
error would be about 4000 ft. The accelerometer error would 
also cause position error of about 500 ft [Ar = Aa,,(t?/2)] 
so that the total miss would be about 4500 ft. If the z 
accelerometer had a similar zero offset error, the missile 
would have a lateral velocity error of 13.5 fps at motor 
shutoff. This would cause a lateral miss of about (13.5 fps) 
(300 sec) = 4050 ft. In a similar manner, the effects of 
accelerometer scale factor errors or of integrator errors can 
be evaluated. 

The effect of gyro drift about the pitch axis can be calcu- 
lated by expressing accelerometer error as a function of atti- 
tude reference error. However, gyro drift in the azimuth 
direction is probably more critical than in the pitch or roll 
direction. In the azimuth case for instance, a constant gyro 
drift rate of 0.3 deg/hr would result in about 1 milliradian plat- 
form azimuth misalignment at the end of 70 sec of flight. This 
would cause the missile velocity vector (actual trajectory 
plane) to point somewhat less than 1 milliradian away from 
the standard trajectory plane. As an upper bound, this drift 
could be considered as occurring at the start of flight, i.e., as an 
initial misalignment of the trajectory plane. It would then 
lead to a target miss of about (0.001 radian) (230 miles) 
(5280 ft/mile) = 1300 ft. 

It is apparent that to achieve good accuracy, precision com- 
ponents are needed. It is indicative of the state of the art that 
components with performance better than the hypothetical 
values assumed here are advertised commercially. Even 
making allowance for all the other errors in a guidance sys- 
tem, it would appear that these commercial components would 
give an accuracy better than was obtained by the German V-2 
missile of World War II. 
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A Deceleration Probe for Measuring Stagnation Pressure 
and Velocity of a Particle-Laden Gas Stream 
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The measurement of the stagnation pressure and, from 
this, the velocity of the gas phase, in a stream of gas laden 
with liquid droplets or solid particles is investigated 
theoretically and experimentally. The three requirements 
of a successful design are as follows: (i) To prevent flood- 
ing of the probe and pressure lines with particles, the rear 
of the probe must be vented by an orifice of appropriate 
size. Moreover, pressurized manometers with a special 
operating technique must be employed. (ii) To minimize 
the effects of aerodynamic interactions outside the probe, 
the latter must be made sufficiently small. (iii) To mini- 
mize the effects of aerodynamic interactions within the 
probe, the internal pressure tap must be situated close to 
the nose of the probe. The phenomena associated with 
the momentum interchanges between gas and particles 
are analyzed, and quantitative evaluation is made of the 
associated error recorded by a total pressure probe. For 
convenience in design the magnitude of the error is pre- 
sented graphically in terms of the probe dimensions 
and the properties of the gas and particle flows. The 
theoretical findings are checked by the testing of a series of 
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experimental probes of varying geometry in an air stream 
carrying water droplets. Satisfactory agreement exists 
between theory and experiment. It is shown that a prac- 
tical probe can be built to give a negligible error in the 
measurement of stagnation pressure in the range of vari- 
ables investigated. Further it is found possible, with 
simple modifications, to adapt the experimental probe con- 
veniently to the local measurement of other stream proper- 
ties, namely, the mass rate of flow per unit area of particles, 
the particle size, and the particle velocity. 


Nomenclature 

A = cross-sectional area 

Cp = drag coefficient of spherical particle, based on relative 
velocity 

d = diameter of droplet 

D = characteristic diameter of probe or tube 

e = local collection efficiency: ratio of actual capture rate of 
particles to capture rate if trajectories were straight 
lines 

i = vV-1 

K = ratio of average gas velocity inside probe to free stream 
velocity 

m = droplet mass 

M = Mach Number 

p = static pressure 

R = probe radius, D/2 

Rey = local Reynolds number of droplet in gas stream (Eq. 9b) 
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V = velocity 

Var = component of gas velocity in x-direction 

Vay = component of gas velocity in y-direction 

Vaco = gas velocity at upstream infinity, assumed in .c-direc- 
tion 

Ver = component of water droplet velocity in x-direction 

Vuy = component of water droplet velocity in y-direction 

View = water droplet velocity at upstream infinity 

w = mass rate of flow 

r,y = Cartesian coordinates 

r = distance measured downstream from probe inlet in 
axial direction 

y = distance measured normal to axis of symmetry 

mM = coefficient of dynamic viscosity 

p = mass density 


Subscripts 


a = refers to gas 

av = average value 

meas = measured quantity 

0 = true stagnation condition of gas phase 
stag = stagnation point 

stat = static 

w = refers to particle flow 

= component parallel to x-direction 

y = component parallel to y-direction 


= at upstream infinity 


Introduction 


N SEVERAL fields of technology one encounters the prob- 
lem of measuring the velocity and total pressure of an aero- 
sol. This aerosol may consist of liquid droplets in a gas as in 
atomization and combustion problems, water injection, flow 
of wet stream, etc.; or it may consist of solid particles, as in 
smokes, pneumatic conveyors, dust separators; or mixtures 
of solid and liquid particles as found in connection with wash- 
ing and scrubbing equipment. 

This paper considers the case when the size of the particle 
is small compared with that of the measuring probe. It en- 
deavors (1) to offer a solution to the practical difficulties asso- 
ciated with the physical presence of the particles in the stream 
such as probe plugging, and (2) to rationalize and evaluate, 
both theoretically and experimentally, the dynamic interac- 
tion taking place between the particles and the gaseous 
stream. The ultimate objective is to present a measuring 
probe of practical operation and capable of measuring the gas 
phase total pressure of the aerosol with a known degree of 
accuracy. 

The experimental work was conducted with an air stream 
laden with small (i.e. in the range of 5 to 50 microns diam) 
water droplets. Accordingly, the particular case of water 
droplets will hereinafter be discussed. The solutions found, 
however, are also valid for small solid particles in a gas stream. 


Preliminary Consideration of Problems of Probe 
Design 


The use of an impact tube gives rise to two basic difficulties. 
Firstly, the instrument and the pressure transfer passages 
may rapidly become fouled with the accumulation of cap- 
tured particles. Secondly, the measured reading will be in- 
fluenced by the momentum interchanges between gas and 
particles both upstream of and within the instrument. 


The Plugging Problem 


Flooding of the probe may be circumvented by venting the 
interior through a suitably-proportioned purging orifice at the 
downstream end. Fig. 1 is typical of the arrangement. For 
water air mass-flow ratios up to 0.25, a smoothly tapered exit 
passage of throat area between three and five per cent of the 
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Fig. 1 Schematic of experimental probes 


R OD 


Fig. 2 Gas streamlines and droplet trajectories near probe 
mouth 


inlet area was found to scavenge the captured water success- 
fully and yet introduce a stagnation pressure error of 0.25 
per cent, at most, of the free stream dynamic pressure. With 
the experimental arrangement shown, internal static pressures 
are recorded through wall taps. 

Fouling of the pressure transfer lines is avoided by using 
pressurized manometers. The manometer pressure is initially 
held above the probe pressure and equilibrium is achieved by 
flow from the manometer to the tube, thus circumventing the en- 
trance of water into the pressure transfer passages. 


Dynamic Interacting External to the Probe 


There exists in front of the probe a complex three-dimen- 
sional gas-flow field (Fig. 2). A droplet, by virtue of its com- 
paratively great inertia, tends to continue in a straight-line 
trajectory with unchanged speed. Upon crossing the gas 
streamlines, however, the droplet is subject to a drag force 
(owing to the vector difference in velocities), and, to a lesser 
degree, to a buoyancy force (owing to the pressure gradient). 
Thus there are imposed on the droplet both a retarding force 
and a radially outward force. The droplet not only decelerates, 
but also migrates away from the axis, with the resulting tra- 
jectory lying between a straight line and the gas streamline 
(Fig. 2). The significant effects of these disturbances in the 
droplet motion are two in number: 

First, because of the divergence of the trajectories, the 
capture or collection efficiency of the probe is less than one. 

Second, any momentum change of the droplets as the latter 
traverse the gas streamlines must cause the gas pressure in 
the probe to be different from what it would have been if the 
gas had been decelerated without aerodynamic forces having 
been exerted on it by the droplets. Generally, the water loses 
momentum ahead of the probe, and this causes a rise in gas 
pressure over and above that which would occur if the particles 
were absent. 

For fixed properties of the gas and particle flows, the dimen- 
sionless error (i.e., the overpressure ratio) depends only on the 
probe diameter D. Some insight into what is important in the 
probe design is gained by considering the two limiting cases of 
very small and very large probe diameters, respectively: 

(a) When the probe is very small. (D — 0), the extent of the 


25 


| 
Inlet 
Pressure 
its, Probe 
ns, 4 
and a 
| 
| 
| 
4 
| 
ri- 
re 
N |_| 


Probe Diameter, D 


Fig. 3 Illustrates schematically the effect of probe diameter on 
error in measured stagnation pressure 


probe disturbance field is correspondingly small. Therefore 
the droplets are exposed to aerodynamic forces for only a very 
short time; accordingly, they suffer only small changes in 
momentum; from momentum considerations, consequently, 
the gas pressure is then relatively unaffected by the droplets. 
Hence, as the probe becomes very small, the error in mouth 
pressure tends to vanish. 

(b) When the probe is very large (D ~ ©), the size of the 
disturbance field is correspondingly large. Because of the 
long time during which aerodynamic forces can act on the 
droplets, the latter tend to follow closely the motion of the 
gas, and the capture efficiency tends toward zero. In the 
limit the gas-droplet mixture may then be thought of as a 
quasi-homogeneous fluid having the average density p.,,[1 + 
(w./Wa)]. Bernoulli’s equation applied to the stagnation 
streamline of an incompressible flow then gives, for the case 
of very large probe diameter 


From the foregoing considerations it follows that the effect 
of probe size on probe error must be generally as portrayed in 
Fig. 3. Although one may conceive of interpreting the probe 
reading at either extreme of such a curve, in practice it seems 
easier to strive for a design in which the probe is so small that 
it reads the gas stagnation pressure directly, with negligible 
error. For design purposes, then the important question is, 
how smal] must the probe be (for specified free-stream prop- 
erties) in order that the error may indeed be negligible? 

Later it is shown that the overpressure curve of Fig. 3 de- 
pends not of course on D alone, but on two dimensionless 
groups (the Reynolds Number and Obedience Number) 
which involve as well the properties of the gas and particle 
flows. 


Dynamic Interactions Within the Probe 


After entering the probe, the droplets find themselves mov- 
ing through relatively stagnant air, and they decelerate 
rapidly. Momentum considerations necessitate a concurrent 
rise in pressure in the gas (Fig. 4), with maximum pressure at 
the point where the velocity of the droplets is reduced to that 
of the gas. For this reason the pressure tap used for measur- 
ing stagnation pressure must be located as near the inlet as 
feasible, provided that gas stagnation conditions prevail at 
that point. How near this must be is analyzed quantitatively 
later. 
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Fig. 4 Longitudinal pressure distribution inside probe (sche- 
matic) 


A Practical Design for a Stagnation-Pressure Probe 


Reviewing the requirements outlined previously, a success- 
ful design requires the following features: 

(a) The captured liquid should be disposed of by some such 
means as a vent. 

(b) The external dynamic interactions should be minimized 
by making the outside diameter small and by suitably design- 
ing the supporting shank and the shape of the inlet. 

(c) The internal dynamic interactions should be minimized 
by placing the internal pressure tap very close to the inlet. 

A practical design for meeting requirements (a) and (c) is 
illustrated in Fig. 5. The probe is self scavenging. Its inlet is 
smaller than the maximum diameter, thus tending to minimize 
external dynamic effects. The ring-shaped pressure pickup 
located near the inlet is protected against water intake by a 
small projecting lip. The outside diameter must of course be 
compatible with the requirements of accuracy. Such probes 
have been built and tested with excellent results (3).4 


Theoretical Studies 


We turn now to the quantitative analysis of the droplet 
motions in the neighborhood of the inlet of a test probe. 
A theoretical study of these events will indicate the orders of 
magnitude to be expected as well as the relative importance of 
the significant parameters. 


Model for Theoretical Analysis 


Exact calculation of the trajectory of a droplet is quite a 
formidable undertaking. It requires accurate knowledge of 
drop size distribution and shape; accurate knowledge of the 
viscous and compressibility effects on the gaseous velocity 
fields and on the interactions between gas and droplets. To 
make the theoretical problem manageable, we adopt rather an 
approximate but simple model which, it is believed, incor- 
porates the leading physical phenomena. The assumptions 
are as follows: 

(a) The gas flow is incompressible (Ref. 2 indicates that 
the effects of compressibility are small), non-viscous, and irro- 
tational. 

(b) The gas flow pattern is unaffected by the presence of 


the droplet. 


(c) The droplets behave like solid spheres with no mutua! 
interaction. 


‘ Numbers in parentheses indicate References at end of paper. 
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Fig.5 Working design of probe 
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Fig. 6 Fluid streamlines of two-dimensional flow, for 
ID/OD =1 
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Fig. 7 Contours of constant fluid velocity components, for axisymmetric tube with /D/OD = 1. 


Linear Scale 


Internal tube velocity equal to 5 per 


cent of free-stream velocity 


(d) The droplets are very small compared to the probe. 

(e) The drag coefficient of the droplets is a function only 
of the Reynolds Number relative to the fluid and varies ac- 
cording to (6). 

(f) The drag forces are the only significant forces acting on 
the droplet. Other forces, namely, virtual mass, gravity, and 
buoyancy can be shown to be small (1 and 5). 

(g) Droplet and gas have the same velocity far upstream. 


Internal and External Gas Velocity Fields 


The problem to be solved here deals with the mapping of 
the incompressible, irrotational, three-dimensional flow field 
uear the entrance of an open-end cylindrical tube, the axis of 
which is aligned with the flow (Fig. 1). The effect of tube-wall 
thickness was considered by studying the two extreme cases in 
which the ratio of internal to external diameter is respectively 
unity (tube with very thin wall) and zero (solid cylinder with 
very small stagnation hole). In addition, venting considera- 
tions call for the former one to have an internal velocity equal 
to 5 per cent of the undisturbed flow velocity. 

Thin-Walled Cylindrical Tube. Approximate construction 
of the flow field within and around a thin-walled tube is based 
on a modification of the potential flow solution for the cor- 
responding two-dimensional channel as obtained from (7). 
Super-position of a uniform velocity 1.05 times as large as and 
opposite in sign to the channel velocity yields the solution to 
the two dimensional problem (Fig. 6). 

Approximation to the three dimensional problem was ob- 
tained by making use of the fact that the relationship between 
the two-dimensional flow and its axi-symmetric counterpart 
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is similar to that between a circular cylinder and a sphere. 
Applying correction factors derived from these well-known 
solutions, the field of Fig. 7 for the flow about a thin-walled 
axisymmetric tube was obtained. 

Thick-Walled Cylindrical Tube (ID/OD = 0). Ina similar 
manner the velocity field for a square-ended, round tube, with 
a very small pressure hole may be derived by modifying the 
solution for the two-dimensional flow past a two-dimensional 
step of height (D/2), (6). For the case 7D/OD = 0, we are in- 
terested only in the velocity distribution on the axis of sym- 
metry. The two-dimensional solution of this is shown by a 
dashed line in Fig. 8. 
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Fig. 8 Fluid velocity distribution on center line, for /D/OD = 0 
Dashed curve is for two-dimensional flow; solid curve for three- 
dimensional flow 
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By applying the correction factors previously mentioned to 
the dashed curve, the solid curve of Fig. 8 is obtained, repre- 
senting the desired center-line velocity distribution for a semi- 
infinite cylinder with 7D/OD = 0. 


The External Motions of the Droplets 


According to the previously made simplifying assumption, 
of the several forces acting upon a droplet, only the aerody- 
namic drag is of consequence in the practically useful ranges 
of the controlling parameters. Equating the drag force to the 
mass-acceleration product, the equations of droplet motion 
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Fig. 9 Typical results of droplet trajectory calculations, for 


axi-symmetric flow, with 1, Reynolds Nubmer _ 
OD Ha 
300 and Obedience Number °°”? = 0.10 
4 pyud 


(a) Lateral deflections of trajectories (see Fig. 2 for nomencla- 
ture) from location at upstream infinity. 

(b) Lateral velocities of droplets in y-direction. 

(c) Reduction in longitudinal velocities of droplets. 
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Fig. 10 Decrease in center-line droplet velocity for axi-sym- 
metric probe with /D/OD = 0 
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in the z- and y-directions may be written respectively as 
‘a 1d? 
Vay — Vuy)? + — Vaz)?] = 
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6 p ( y ) [2a] 
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We, 
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V (Vay — Vey)? + (Vex — Vas)? 
Now, for a good probe design, the capture and momentum 
efficiencies should be nearly unity, and the droplet trajectories 
are, accordingly, nearly straight. On the right-hand sides of 
the equations of motion, therefore, we may drop the terms in- 
volving V OV »:/dy and V.,OV ,,/d, without sensible error in 
the ranges of practical interest. Making this approximation, 
and reforming the equations in dimensionless form, we obtain 
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Droplet trajectories are defined by the equation 
d(y/D) 
[4] 


d(t/D) Vuz/Vaw 


In Equations 3, the drag coefficient ¢p depends on the rela- 
tive Reynolds Number of the droplet 


ped / (Vey — + (Vor — Vaz)? 
Rey = — [5 | 
Ma 

A study of the form of Equations 3 then shows that two di- 

mensionless parameters (apart from the probe geometry) 

control the droplet motion, namely, (i) the characteristic 

Reynolds Number Rey, = paVwaod/ Ma and (ii) the obedience 

Number, 3p.D/4p.d, so named because its magnitude is a 

measure of how closely the droplet trajectories obey the gas 
motion. 

Simultaneous solution of Eq. [3-4] using the drag-coefficient 
data of (6), yields a complete description of the motions of the 
droplets. These solutions were carried out numerically with 
the method of isoclines. Extensive numerical results are pre- 
sented in (5). Typical results of practical significance are il- 
lustrated in Fig. 9 (for 7/D/OD = 1.0) and Fig. 10 (for 7D/OD 
= 0). 


Capture Efficiency 


Thin-Walled Tube. From the data of Fig. 9a it is an easy 
matter to compute the capture efficiency for the thin-walled 
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tube with JD/OD = 1.0. The fraction of uncaptured water is 
shown in Fig. lla as a function of Reynolds Number and 
Obedience Number. For (#/D) = 0, the curves show the 
fraction of the water in the capture area 7D?/4 which escapes 
entering the probe mouth. For positive values of x/D (cor- 
responding to stations within the probe), the ordinates show 
the fraction of water in the capture area which is no longer in 
the gas stream at the abscissa x/D; this fraction is accounted 
for by the amount escaping the mouth of the probe plus the 
amount impinging on the internal walls from z = 0 tow = 2. 

Thick-Walled Tube. In the case of the thick-walled tube in 
which /D/OD — 0, the capture efficiency at the stagnation 
point cannot be calculated from an analysis involving the 
central streamline only. An indeterminate condition is pres- 
ent there. Removal of the indeterminacy, yields the following 
expression for the capture efficiency 

1 2 


d(y/R) (y/R) =0 
[t is assumed here that the droplet behavior along the central 
streamline will be closely similar for the cases of 7D/OD = 0 
and JD/OD = 1.0. This assumption enables us to obtain 


the desired information from charts like Fig. 9a. The results 
for capture efficiency are shown in Fig. 11b. 


=> 
(y/R) =0 


Overpressure in Tube Mouth 


The average overpressure in the tube mouth, which is 
probably a good interpretation of the reading of a static tap 
very close to the mouth, depends on both the momentum 
efficiency and capture efficiency. Consequently, from the 
considerations leading to Figs. 9 and 10, we may expect the 
dimensionless overpressure 


(14) paVae 2(Wu/Wa) oo 


to be a function of the probe geometry, the Reynolds Number, 
and the Obedience Number. 


Thick-Walled Tube (ID/OD = 0) In this case we need 
consider only the center-line stream tube. Applying the 
momentum equation to an infinitesimal length of this stream 
tube, we have 


—dp = peVartVer + dV 


which may be rearranged in dimensionless form as 
a Ves 
—d ( = ( + 
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The first term on the right-hand side represents the contribu- 
tion from the change in momentum of the gas (assumed to be 
the same as when the droplets are not present); the second, 
that of the droplets. When Eq. 6 is integrated from x = — © 
to x = 0, the integral of the term in dV q; yields the difference 
between py and p,,, the former representing the true stagna- 
tion pressure of the gas phase alone. The integral of the term 
in dV. gives the overpressure in the tube mouth due to the 
droplets. Thus we have 


(=) 
Wa/ w 


where @ is an average value of e over the interval of integra- 
tion, from z = —@ tox = 0 
When the Obedience Number is very small (e.g., very large 
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Fig. lla Uncaptured fraction of droplets for probe with /D/OD 


= 1. Curves for = = 0 show amount of water escaping probe 
mouth. Curves for positive values of - show amount of water 


escaping probe plus amount impinging on internal probe walls 
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Fig. 11b Uncaptured fraction of droplets on center line for probe 
with ]D/OD = in plane = 0 
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Fig. 12 Overpressure in mouth of thick-walled probe, for /D/OD 
=0 
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Fig. 13 Control surface for calculation of average overpressure 
in probe mouth for 7D/OD = 1 
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Fig. 14 Average overpressure in mouth of thin-walled probe, for 
ID/OD = 1 
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droplets, or very small probe), the droplets move with un- 
changed velocity. Consequently the right-hand side of Eq. 7 
vanishes, and the probe error is zero. 

When the Obedience Number is very large (e.g., very small 
droplets, or very large probe), the droplets have the same 
motion as the gas. Consequently the mixture behaves as a 
quasi-homogeneous fluid with the effective density p.[l + 
(wW/Wa)]. In Eq. 7, (Vwz)z-0 = 0 and the dimensionless 
overpressure is equal to 22. For this case, however, the 
dimensionless over-pressure is equal to unity (Eq. 1), and 
consequently é = 0.5. The latter is the mean value between 
1.0 (atz = — ©) and 0 (atz = 0). 

For arbitrary Obedience Numbers and Reynolds Numbers, 
Eq. 7 may be used in conjunction with Fig. 11b to give the 
results shown in Fig. 12 for the dimensionless overpressure. 

Thin-Walled Tube (ID/OD = 1.0). Assuming that the gas 
motion is unaffected by the droplets, we may compute the 
average overpressure in the probe mouth by application of the 
momentum theorem to the droplets passing through the con- 
trol surface shown in Fig. 13 
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The first term on the right-hand side represents the incoming 
flux of droplet momentum through BC; the second term the 
outgoing flux through DE; and the third term the outgoing 
flux through the cylindrical surface represented by CD and 
BE. 

In evaluating the first integral, we consider the elementary 
flow through an annular-shaped control volume which far up- 
stream is bounded by y,, and y,, + dy... Then 


Similarly, the elementary flow for the second integral may be 
written as 


Substituting the latter two expressions into the momentum 
equation, rearranging in dimensionless form, and taking due 
note of the necessary changes in limits of integration when we 
change from y and z as variables to y,,, we get 

Pmeas 


Ww 
y/R= 
0 Vas z=) R 
1 
= _ Vue Yo 


for the overpressure ratio. 

If the momentum-defect curves of Fig. 9¢ are replotted 
against (y/R)?,, then the dimensionless overpressure given by 
Eq. 8 is simply twice the area under the curve formed by the 
line z = 0 joined to the line (y/R) = 1.0. Thus Eq. 8 has a 
simple physical significance: The overpressure is equal to the 
momentum loss of all the water originally in the projected 
cross section of the probe. 

By carrying out the indicated integrations in Eq. 8, the 
over-pressure curves of Fig. 14 are obtained. These curves 
also approach an asymptote of 1 as the Obedience Number 
becomes very large or the Reynolds Number very small. 

Figures 12 and 14 give working curves for design use. They 
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Fig. 15 Universal curves for determining deceleration distance 

of droplet within probe. Difference in ordinates corresponds to 

distance required to reduce droplet velocity from the value 

corresponding to initial Reynolds Number to that corresponding 
to final Reynolds Number 


show directly and quantitatively the influences of the probe 
dimensions and of the gas and particle streams on the error in 
reading at the probe mouth. A study of Figs. 1la, 12 and 14 
show that the error in probe reading may be reduced by using 
a thin-walled probe. Moreover, the error decreases and the 
capture efficiency increases when the probe diameter, the gas 
density and the gas viscosity decrease or the gas velocity, the 
droplet diameter and the droplet density increase. 


The Internal Motions of the Droplets 


Having crossed the complex flow field outside the probe 
inlet, the captured droplets enter a zone of comparatively low 
and uniform velocity. It is of considerable interest to study 
this region inside the probe for two reasons: First, some indi- 
cation may be had as to the magnitude of the error recorded 
by a tap situated too far from the entrance. Secondly, the 
particle size may be determined by a novel method involving 
the measurement of the longitudinal pressure gradient within 
the probe, together with measurement of local droplet- 
velocity and local water-air ratio. 

One-Dimensional Analysis of Internal Droplet Motion. The 
three-dimensional aspects of the internal droplet dynamics 
have already been presented in Fig. 9. For many purposes, 
however, particularly with a well-design probe of high capture 
efficiency, the gas and droplet motions may be considered one- 
dimensional. For this simplified model the equations of mo- 
tion may be written as 

— KVa,,)? wd? 
D 


Cc 
2 4 6 dx 


where KV, is the average gas velocity in the probe. 
Separating variables in this equation, introducing the two 
Reynolds Numbers 
Ve KV, d Vaw aod 
Ka Ma 


[9b] 


and integrating from z = 0 to x = 2, we obtain 


Zaz 
Pow _ | [10a] 


Pud = 33 Jz Cn(Rey)? 
where 
Z=Rey + K Rey, (2 ) [10b] 
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Fig. 16 Longitudinal pressure gradient within probe, reckoned 
at x/D = 0.5, for thin-walled probe with JD/OD = 1 


Now, for a specific value of A Rey , ( se , the integral of 


Kq. 10a is a function of Rey only. Consequently, Eq. 10a, 
when integrated, gives the distance x required for decelerating 
the droplets from the speed corresponding to (Rey), to that 
corresponding to (Rey),—:. 

The results of numerical integration, using drag data for in- 
compressible fluids, are shown in Fig. 15. For the lower limit 
of integration, Rey = 1000 was chosen, so that the curves give 
directly the distance required to reduce Rey from 1000 to any 
lower value, however, the curves are also applicable to any 
initial value of Rey less than 1000. It is merely necessary to 
take the difference in ordinates corresponding to the initial 
and final Reynolds Numbers. 

Pressure Gradient Within Probe. The longitudinal pressure 
gradient within the probe may be related to the deceleration 
of the droplets by writing the momentum equation for a 
cylindrical control volume bounded by the tube walls and of 
length dr. Thus 


d av. 


Now, rearranging and substituting for — the value obtained 
dz 


from Eq. 9a, we get 
1 dp 


3 pPaD Ve 
—K 


where cp is a function of (Rey)o, measured at stagnation 
conditions within the probe. 
Equation 12 has two important uses: 
(1) It shows how large an error in pressure may be incurred 
by locating the pressure tap too far from the probe mouth. 
To illustrate orders of magnitude, calculated results are shown 
in Fig. 16 for the probe with 7D/OD = 1, based on the values 
7 


of e and 

(2) It suggests that the drop size may be inferred from 
measurements of the pressure gradient near the mouth, the 
droplet velocity entering the probe, and the droplet flow per 
unit area entering the probe. Means of measuring these 
properties by modifications of the present method are indi- 


cated in (5). 


computed previously, and reckoned at 2/D = 0.5 
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Fig. 17. Details of experimental probe, Configuration V 


Experimental Studies 


With the object of verifying the results of analysis, an ex- 
perimental program was undertaken in the form of systematic 
testing of several probe geometries over a range of free-stream 
conditions. Of prime interest was the verification of the 
effects of probe diameter and of the ratio of internal to external 
diameter on the pressure in the plane of the probe inlet. In 
addition, measurements of internal pressure gradients were 
desired, and the effective operation of the venting system was 
to be demonstrated. 


Probe Geometries 


Table I summarizes the various geometries investigated. 
All configurations except VITA, VIIB, and VIIC were of the 
type illustrated in Figs. 17 and 18, with internal pressures re- 
corded by wall static taps. In configurations VIIA, VIIB, 


Table I Probe Geometries 


| 


Probe Outside probe 
configuration diam (in.) ID/OD 
| I 0.020 0.500 | 
II 0.035 | 0.657 
IIA 0.109 0.211 
IIB 0.245 0.0938 
Ill 0.109 0.780 | 
IIIA 0.245 0.347 | 
IV 0.177 0.726 | 
| IVA 0.245 0.524 
| Vv 0.245 0.746 
| VI 0.245 0.930 | 
| VIIA 0.350 0.838 
VIIB 0.350 0.492 
VIIC 0.350 0.217 | 


and VIIC, an effort was made to secure more readings with 
the help of a static tap installed in a needle traveling parallel 
to and along the inside wall of the probe proper. All probes 
incorporated a venting aperture. 


Test Tunnel 


Use was made of a test tunnel of transparent lucite, 21/3 in. 
in diam and 72 in. long, and located in the Gas Turbine 
Laboratory at Massachusetts Institute of Technology. The 
tunnel is subatmospheric, air flow being induced by a down- 


32 


Fig. 18 Photograph of several experimental probes 


stream steam ejector. Maximum rate of airflow is about one 
lb per sec, which condition chokes the tunnel and yields a 
Mach Number of about 0.85 in the test plane near the exit 
end. A water injection system of 0.350 lb per sec capacity, 
with parallel downstream injection is installed near the tunnel 
inlet bellmouth. The atomization is done by the air stream 
with atomization Mach Numbers up to 0.6. In every case the 
probe was placed on the axis of the tunnel in the same longi- 
tudinal plane, and was carefully aligned with the axis. 


Test Procedure and Data Reduction 


Comparative testing of the various probe configurations 
calls for subjecting all instruments to identical Mach Num- 
bers. It may be shown that this condition can be achieved 
at the test section by maintaining constant (i) the ratio of the 
static pressure at the test plane to the tunnel inlet pressure 
(atmospheric); (ii) the water-air ratio; and (iii) the inlet air 
humidity, air temperature and water temperature. Dis- 
turbances of any of these quantities, through energy effects 
associated with evaporation or condensation of water and heat 
flow between air and water produce perturbations of the flow 
properties in the test plane. 

Conditions (i) and (ii) were easily achieved in the experi- 
mental procedure. Control over humidity and temperatures, 
however, was not feasible. Instead, theoretical corrections 
were applied to test data, on the basis of one-dimensional 
equilibrium flow in the tunnel (5 and 9). Application of the 
various corrections in general has the effect of reducing the 
data to some standard condition; this condition is here taken 
as the one in which thermal energy interchanges between air 
and droplets are absent. 

Reduction to standard conditions further calls for a correc- 
tion of the error introduced by the presence of the scavenging 
vent. This requires a calibration of each instrument against 
a conventional Pitot tube in a dry stream. The smaller the 
instrument, the more essential this calibration becomes be- 
cause of the difficulty of controlling the vent area with ac- 
curacy. In general, the deviation of the probe reading from 
that of the standard Pitot tube was within 1 per cent of the 
free-stream dynamic pressure. The calibration has the 
further advantage of yielding useful information on the 
characteristics of each static pressure tap. For example, a 
tap located too near the inlet may be found to read incor- 
rectly. 

In a stream containing water droplets, a great deal of atten- 
tion must be devoted to the technique of keeping pressure 
taps free of water. A scheme used successfully involved the 
use of pressurized manometers. Pressurization followed by 
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Fig. 19 Typical pressure measurements near probe mouths with 
Pe — 0,648, = 0.19, M,, = 0.675 


Patm Va 

Inlet Inlet 

dry wet Inlet 
Configura- Atmospheric bulb bulb water 


tion pressure temp temp temp 
Fig. 19a IV 30.216”Hg 70.5F 54F 51F 

IVA 30.146 72 53 40 
Fig. 19b VIIA 30.180 75 58 56 


establishment of equilibrium through a rigidly-enforced one- 
way flow from manometer to tunnel was the solution to the 
problem of keeping the pressure-transfer passages free of 
water. 


Test Results 


Figs. 19a and 19b show typical static-pressures recorded 
within two of the experimental probes. The data of Fig. 19a 
were obtained with probe configurations IV and IVA, having 
individual pressure taps. Fig. 19b refers to probe VITA, con- 
taining a longitudinally-traversing static tap. Point M 
represents the extrapolated measured pressure in the plane of 
the inlet. Point N represents the inlet-plane pressure after 
correction to standard air flow and water flow conditions. 

Data such as those shown in Fig. 19 were obtained for each 
of the probes of Table I over a range of water-air ratios and 
air Mach numbers. In each case the stagnation pressure 
measured by the probe was interpreted to be that found by 
extrapolating the recorded static pressures to the inlet plane 
and then correcting to standard corrections (e.g., point NV in 
Fig. 19 is taken as the measured stagnation pressure). A grand 
summary of all measured stagnation pressure—taken in the 
foregoing sense—is presented in Fig. 20. Each of the four 
charts shown in Fig. 20 is for a particular combination of 
water-air ratio and Mach Number. On each chart may be 
seen the effects of probe diameter and of the ratio of internal 
to external diameter. 

In general, the data in Fig. 20 fan out in a triangular band 
from an apex on the vertical axis. Following the reasoning 
associated with Fig. 3, the position of this apex (on the vertical 
line D = 0) represents the true stagnation pressure. It should 
be noted that such an inferred value of the true stagnation 
pressure was necessary because there is no absolute method 
of determining the latter. 

Apart from the experimental results, there are shown in 
Fig. 20 the theoretically predicted error as found from Figs. 
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Symbol Config. 1.0/0.0. 


° I 0.500 Theoretical Variation for 
ao Probe with an 1.0./0.D 
s 
o 0780 icated Droplet Diamete 
IITA 0.347 
va 0.528 ———— Theoretical Variation for 
a 0.746 o Probe with an 1.0/0.0. 
< VI 0.930 of 0.0 and for the In- 
> VIIA 0.838 ; 
vig 0.492 dicated Droplet Diameter 
VIIC 0.217 
1.477 Microns+ = 1.3875 =10 Microns 
Pip 0.10 020 030 0.10 020 0.30 
Pstotic (a) (b) 
+4218 Microns 
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D, Probe Outside Diameter (inches) 


Fig. 20 Summary of corrected stagnation pressures read by the 
different probe configurations 


Chart Po/Patm (Ww/Wa) M. 
(a) 0.620 0.128 0.744 
(b) 0.648 0.204 0.673 
(c) 0.734 0.182 0.566 
(d) 0.810 0.160 0.469 


12 and 14. In order to employ these charts, it is necessary to 
know the volume-surface mean diameter of the droplets. The 
latter was estimated from Fig. 16, using measurements of the 
internal pressure gradient near the probe mouth (see Fig. 19, 
for example), together with independent measurements of 
water flow per unit area. The droplet diameters so found are 
in reasonably good agreement with those of Nukiyama and 
Tanasawa (8). The curves in Fig. 20 show the predicted 
probe errors for the range of droplet sizes in each test found as 
described, and for the two limiting cases of /D/OD = 0 and 
ID/OD = 1.0. 

The differences in stagnation pressure measured by the dif- 
ferent probe configurations are significant but are not large 
compared with the experimental accuracy. Consequently a 
considerable amount of apparent scatter in the data is intro- 
duced when these are plotted, as in Fig. 20, on a scale large 
enough to display differences. It is evident that the theoreti- 
‘al predictions are in general agreement with the test results: 

(1) The overpressure increases as the probe size increases. 

(2) The overpressure increases as the ratio 7D/OD de- 
creases. 

(3) The magnitude of the measured error is of the same 
order as that given by the theory. 

Since, in the range of variables tested, the errors in probe 
are small, it is not considered of practical importance to pre- 
dict the error with great accuracy. Moreover, the difficulties 
of precise experimental measurement preclude the possibility 
of accurately observing such small errors. 

Fig. 20 shows finally the feasiblity of constructing a stagna- 
tion pressure probe of convenient size and yet of satisfactory 
accuracy. 


Conclusions 


1. A simple, practical probe design capable of measuring 
the stagnation pressure of the gas phase in a droplet-laden 
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stream has been developed. Its principle of operation has 
been formulated and verified experimentally. 

2. A theory of probe operation is presented, adaptable to 
any stream carrying suspended bodies of density greater than 
that of the stream. Generalized curves have been worked 
out providing design data for any such application. 


Acknowledgments 


The work described here was supported partly under Con- 
tract N5ori-07878 with the Power Branch of the Office of 
Naval Research and partly by the Procter and Gamble Fel- 
lowship at MIT. It was carried out in the MIT Gas Turbine 
Laboratory with the heplful cooperation of the Laboratory 
Staff. 

It is a pleasure to acknowledge particularly the contribu- 
tions to this work of Prof. Arthur A. Fowle, Harry Foust, 
Dalton Baugh and Margaret Tefft. 


References 


1 Boucher, R. M. G., ‘Sur le fonctionnement de l’épurateur a 


microbrouillards Aerojet-Venturi,’’ Chaleur et Industrie, no. 328, 
November 1952. 

2 Brun, R. J., Serafani, J. S. and Gallagher, H. M., ‘Impinge- 
ment of Cloud Droplets on Aerodynamic Bodies as Affected by 
Compressibility of Air Flow Around the Body,’’ NACA TN 2903, 
March 1953. 

3 Carpenter, James W. and Martin, W. L., ‘Construction and 
Preliminary Test of a Traversing Stagnation Pressure Probe for a 
Large Scale Aerothermopressor,’’ Thesis for Naval E. Degree, 
May 1954, MIT, Cambridge, Mass. 

4 Cheatham, John B., Jr. ‘The Photometric Measurement of 
Particle Size,’’ Thesis for M.E. Degree, September 1954, MIT, 
Cambridge, Mass. 

5 Dussourd, Jules L., ‘‘A Theoretical and Experimental In- 
vestigation of a Deceleration Probe for Measurement of Several 
Properties of a Droplet-Laden Air Stream,’’ Sc. D. Thesis, De- 
partment of Mechanical Engineering, MIT, 1954. 

6 Eisner, Widerstandproblem,’’ Third International 
Congress of Applied Mechanics, Stockholm, 1930. 

7 Milne-Thompson, L. M., “Theoretical Hydrodynamics,”’ 
MeMillan and Co., London, 1938. 

8 Nukiyama, 8S. and Tanasawa, T., “‘An Experiment on the 
Atomization of Liquids by Means of an Air Stream,”’ Transactions 
ASME, Japan, vol. 4, Report 1-6, (1938). 

9 Shapiro, A. H., “The Dynamics and Thermodynamics of 
Compressible Fluid Flow,’’ Ronald Press, New York 1953. 


Mass-Transfer Cooling of a Laminar Boundary Layer by 
Injection of a Light-Weight Foreign Gas 


E. R. G. ECKERT,’ P. J. SCHNEIDER,’ A. A. HAYDAY? and R. M. LARSON® 


University of Minnesota, Minneapolis, Minn. 


Analytical predictions are given for the development of 
the velocity, temperature and concentration fields in a 
laminar air boundary layer on a flat plate in high-speed 
dissipative flow, the plate being considered porous and 
cooled by injection of hydrogen from its surface. The 
admixture of hydrogen, having a low density and high 
thermal capacity relative to air, is shown to greatly di- 
minish the skin friction and to markedly relieve the ad- 
verse thermal effects of intense aerodynamic heating under 
conditions of hypersonic flow. 


Nomenclature 
A = Kye — 1)M*. 
Cr = integration constants 
Cp = constant-pressure thermal capacity 
Cy = constant-volume thermal capacity 
Che = local skin-friction coefficient = (2u./p.u?.)(Ou/OYy)w 
Dy» = coefficient of ordinary diffusion 
f(n*) = modified stream function 
tw = wall injection parameter = Rez 
h = unit surface conductance 
k = mixture thermal conductivity 
M; = molecular weight of i’th component in mixture 
M. = free-stream Mach number 
Nu = = local Nusselt number = hz/k, 
= free-stream Prandtl number = 
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q = heat flux per unit area 

Re, = local free-stream Reynolds number = u.2/v, 

Sc = Schmidt number = v,,/Dy 

a = normalized mixture temperature = t/t,,, (t. = 392°R) 

t = local mixture temperature (absolute) 

ty = recovery temperature 

u = velocity component parallel to flow (x-direction) 

v = velocity component perpendicular to plate surface (y- 
direction) 

W = mass fraction of coolant in binary flow 

x = coordinate parallel to plate surface 

y = coordinate perpendicular to plate surface 

Yo = specific-heat ratio of free-stream fluid = cp./cve 

Ne = Blasius similarity coordinate = }y(u../v.2)3 

n* = similarity coordinate defined by Eq. [5] 

m = dynamic viscosity of mixture 

vs = kinematic viscosity of free-stream fluid = p./p. 

p = mass density of mixture 

= recovery factor = 2(7; — 1)/(v2 — 1)M,? 

ec = normalized mixture thermal capacity = cp/cp. 

¢a2 = normalized pure-component thermal capacity difference 

Cp. — 

= 

¢k = normalized mixture thermal conductivity = k/k, 

= normalized mixture viscosity = 

Yp = normalized mixture density = p/p. 

= stream function = Vy 

Subscripts 

1 = foreign component in binary mixture (hydrogen) 

2 = primary component in binary mixture (air) 

r = recovery 

w = wall(y = 0) 

o = free stream (y > ~) 


Primes (’) denote differentiation with respect to n* 
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I. Introduction 


HE USE of mass-transfer cooling is being investigated as 
a means of protecting the surfaces of high-speed aircraft, 
missiles and satellites exposed to intense aerodynamic heat- 
ing. The need for cooling such surfaces becomes critical under 
one of three flight situations: (a) Rapid acceleration at con- 
stant or varying altitude, (b) sustained flight at high Mach 
numbers, and (c) re-entry flight into the earth’s atmosphere. 
The essential mechanical feature of mass-transfer cooling is 
represented by injection of a foreign gas into the boundary- 
layer fluid. This feature distinguishes it from conventional 
cooling methods which depend on either a high thermal- 
capacity heat sink within the structure or artificial cooling at 
the inner surface of the skin. If the injection gas is identical 
with the boundary-layer fluid, then the method is known as 
transpiration cooling, while if the injection and mainstream 
fluids are dissimilar the method is called mass-transfer cooling. 
The foreign injection fluid may be a gas, or it may be a liquid 
which covers the porous surface as a film and evaporates into 
the boundary layer. An alternate technique is represented by 
sublimation cooling in which the walls of the vehicle are con- 
structed of a solid material whose surface continuously 
evaporates into the high-temperature boundary-layer air. In 
all such processes the heat and momentum transfer are ac- 
companied by a simultaneous diffusion of the injection or 
venerated gas away from the surface. This diffusion process 
has a marked effect on the boundary-layer development, par- 
ticularly when the properties of the diffusing foreign gas are 
widely different from the properties of the main stream. 
Results of an exact calculation were reported in (1) for the 
velocity and concentration fields in a laminar isothermal 
houndary layer on a flat plate with locally uniform pressure 
and with injection of a foreign gas from the surface, assuming 
an uncontaminated coolant reservoir in which no air penetrates 
the surface into the interior of the plate. In that analysis the 
fluid properties were allowed to vary with local mixture con- 
centration in the binary isothermal boundary layer, and 
numerical calculations were performed for injection of hydro- 
gen into air. Hydrogen was selected as a coolant because its 
properties are widely different from those of air. The results 
of that study therefore indicated the strong effect of property 
variations on the mass-transfer process. The purpose of this 
present paper is to report on similar calculations for a high- 
velocity boundary layer on a cooled plate. All fluid proper- 
ties in the non-isothermal case are again considered variable, 
the dependence now being on both mixture concentration and 
mixture temperature. 


II. Differential Equations and Boundary 
Conditions 


The differential equations for a non-isothermal binary 
boundary layer with variable fluid properties and in the 
absence of a pressure a are (2) 


ou 
ow ow re) ow 
Diffusio1 »— — — [| —) = 0........ 3 
sion pu + pr Oy 12 ~) [3] 
Energy 
at +. at od (: 
ou ot OW 


® Numbers in parentheses indicate References at end of paper. 
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The effect of thermal diffusion (mass diffusion in the presence 
of a temperature gradient) has been neglected in Eq. [3-4]. 
The last term in the energy equation is absent in the con- 
ventional single-component boundary layer; it represents the 
diffusion of enthalpy in a binary mixture. 

The continuity equation is satisfied by introducing a 
stream function ¥. Eq. [2-4] are then reduced to a corre- 
sponding set of ordinary simultaneous differential equation 
through a similarity transformation * defined by 


0 Pu 


and a modified stream function f(n*). The transformation 
n* is related to the familiar Blasius similarity coordinate n, 
= (y/2)(u../v.x)' by the relation 


n* 
0 


which is valid for any fixed x. The transformed results 


corresponding to Eq. [2-4] are 


Momentum 

Diffusion 

Energy 


+ = 0. .[9] 


In these equations the primes (’) denote differentiation with 
respect to 7*; all other symbols are defined in the Nomencla- 
ture. The mixture properties Sc, Gu, Ge ANd de- 
pend, in general, on both the local mixture concentration W 
and local mixture temperature 7 (Article II). 

Equations [7-9] must satisfy a total of seven boundary con- 
ditions, all of which are the same as in [1] except for two ad- 
ditional conditions on the energy equation 


=0 =2 
= 0,W n*— @...... 


The second condition at »* = 0 implies that there is no mass 
flow of air through the wall (1). Although a diffusion of 
boundary-layer air into the plate occurs because of the finite 
concentration gradient at the surface, it is assumed that the 
diffusion velocity of this air is small compared to the convec- 
tive velocity of the coolant in the porous plate, and that ac- 
cordingly the mass of air that does enter the plate is negligible. 
Thus (1 — = where vo = 
—[( Gp) Transformation of this to gives 
the second condition in Eq. [10], with C, defined by Eq. [15]. 
Either the surface injection rate f,, or the surface mass con- 
centration W,, may be chosen arbitrarily. 


III. Property Values 


Suitable methods for calculating the concentration and 
temperature dependence of the various transport properties 
appearing in Eq. [7-9] have been summarized in (8). 

1. Mixture Schmidt number, Sc = v./D2; assumed de- 
pendent on 7’ only (4). Values of Se calculated by Hirsch- 
felder’s method (3) are shown in Fig. 1 for t,, = 392 R, v. = 
0.340 ft?/hr, and 1 Z T Z 10. 

2. Pure-component thermal capacity, and selected 
from (5) for molecular hydrogen and dry air. Calculated 
values of ¢.,, = (cp, Cp)/€po are shown in Fig. 1 with 
Cpa = 0.239 Btu/Ib-R. 

3. Mixture density ratio, ¢,; calculated assuming per- 
fect-gas behavior as 
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Fig. 1 Variation of Schmidt number Sc and pure-component 

thermal capacity difference ¢,,, with mixture temperature 7’ = 
for hydrogen-air mixture 


Gp = + (Mo/M, WIT = 1/11 + 13.5W)T 


4. Mixture viscosity ratio, gy; calculated by Wilke’s 
method (3) as 


~ 1 + 0.129611 — W)/W 1 + 3.9730W/]1 — W) 
the pure-component viscosities 4; and yw: being selected from 
NBS-NACA tabulations (6), Tables 7.39 and 2.39. Calcu- 
lated values of g, = uw/m. are shown in Fig. 2 with uw. = 
0.0345 Ib/hr-ft. 

5. Mixture conductivity ratio, ¢,; calculated by Enskog’s 
method (3) as 


Pu 


+ k2/a2 + 0.2890 VW 
(1 — 0.0209 /ara2)k 

where a; = 1 + 0.1886(1 — W)/W and az = 1 + 14.2280W/- 
(1 — W), the pure-component conductivities k; and k, being 
calculated by Hirschfelder’s method (4). Values of g, = 
k/k.. are shown in Fig. 3 with k,, = 112.6  10~4 Btu/hr-ft-R. 
6. Mixture thermal capacity ratio, ¢,; calculated as- 

suming perfect-gas behavior for each component as 


Coe Coo 


IV. Method of Solution 


Equations [7-9] together with the boundary conditions in 
Eq. [10] can be transformed into the following system of 
simultaneous integral equations 


Gc = 


Momentum 


n* 
* fan* 


Diffusion 


* (euSe/ go) fdn* 


Energy 


n* 
n* -f X(n*)dn* n* 
T=T,+ JO dn* — f, (Pu/ 


n* n* 
= A(n*)dn* (°n* + A(n*)dn* * * ‘ 
f, f, 0 dyn*dy*.. . .[13] 


H,~ Air 
4 
3 | 
2 
/ 
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Fig. 2 Variation of mixture viscosity », with mass concentration 
W and mixture temperature 7’ for hydrogen-air mixture 


where 
-f fdn* 
C= 2/ e Jo [14] 
n* 

and 


n* 
-f, A(n*) dn* 


* 


* A(n*)dn* 
g(n* Je dn*dn* 


C3 = ..[16] 
f A(n*)dn* 
f, Jo dn* 


In these integral equations, the \(7*) and g(n*) are defined as 
A(n*) = + 
and 
g(n*) = 


wherein Pr., = 0.73 and the A is constant. Of the six im- 
posed boundary conditions, three appear as constants C,, Co, 
and C3, and the remainder are surface values f,, wy, and 
T,. The explicit relationship f,, = C./(1 — W..) is obtained 
directly from Eq. [2, 3, and 12]. If the boundary condition 
is one of an adiabatic wall (7’),, = 0, then Eq. [11-12] re- 
main unchanged and Eq. [13] changes to 


n* 
n* -f Mn*)dn* 
T=T,.- f, (gu/gnle JO 


n* 
* + AC *) d * 


where 


n* 
Te = 1+ (en/erre Sc 


n* 
An*)d 
g(n*)e Si dntdn*... (18! 


The integral Equations [11-13] or [17] were solved by suc- 
cessive approximations. A reasonable guess at the functions 
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Fig.3 Variation of mixture conductivity ¢, with mass concentra- 
tion W and mixture ‘emperature 7’ for hydrogen-air mixture 


f’(n*), W(n*), and T(n*) together with a selection of W,, and 
T,, and an estimated f,, constituted the initial approximation. 
It was found that rapid convergence of successive iterations 
could be achieved by first iterating on f’(n*) alone, using Eq. 
[11] and holding W(n*) and 7'(n*) fixed. This procedure 
vields a set of approximations which represent consistant 
guesses. These generated guesses were then used as first ap- 
proximations for a complete iterative cycle in which ordered 
changes in all functions were made simultaneously. The cycle 
which yielded the most rapid convergence proceeded in the 
order (W, 7,7 2,7, 

Machine Calculations. The iterative method described 
above formed the basis of solution on a Univac Scientific 
Computer. Operations were performed on groups of numbers 
stored in the computer, each group representing a column of 
values for the index cAn*. All numbers were in scientific nota- 
tion with eight significant digits and an exponent being stored. 
A floating-point interpretative system (FLIP) developed by 
Convair was used for the basic arithmetic, and a pseudo-in- 
terpretative vector-operation library was developed to handle 
operations on the columns of numbers. These systems greatly 
simplified the coding. It should be remarked in this connec- 
tion that early attempts using the conventional Blasius trans- 
formation were difficult to code because a perturbation method 
was required to force convergence. Approximately 65 sec 
were used to set up the program for a particular case, and 45 
sec were required for each iteration and monitoring. Each 
solution was carried through an average of 10 iterations. 

A tabular interval of An* = 1 was used in all solutions. 
However, repeat solutions were carried out in several cases 
with An* = 14 and 4 but no significant changes were noted 
with the finer subdivisions. Integration of the tabular entries 
was performed by the trapezoidal rule. The iterative scheme 
used was such that round-off errors were damped out; thus 
the sources of error are due only to inaccuracies in property 
values, the integration formula, and incomplete convergence. 
An integration of the difference between two consecutive 
solutions provided a convergence criterion; if this result was 
less than 10~4 for each of the functions f’, W, and 7’, then the 
solutions were terminated. With this criterion the last two 
approximations in each solution agreed point-by-point to four 
significant figures. 

7 The particular transformation Eq. [5] was chosen after at- 
tempts to solve the problem with the usual Blasius transforma- 
tion Eq. [6] failed owing to lack of convergence. 
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V. Results 


Typical results of the calculation are shown in Fig. 4 for in- 
jection of hydrogen into air, all results being plotted against 
the Blasius coordinate nx = (y/2)Wu./vox. The upper 
set of curves represent temperature profiles for M.. = 12, T., 
= t,/t. = 2 and 6, and W,, = 0.0, 0.1, 0.2, 0.4, 0.6 and 0.8. 
The Van Driest curve (7) for = 6 and W,, = 0 (zero injec- 
tion) differs from the present results for W, = 0 because of 
different assumptions on the temperature dependency of pure 
air. The temperature profiles for finite injection rates indicate 
that the peak boundary-layer temperatures are substantially 
reduced by injection of a light-weight gas of high thermal 
capacity. The surface temperature gradients are seen to be 
greatly reduced at this high Mach number even with a com- 
paratively low surface mass concentration of W,. = 0.1. This 
corresponds to a high reduction in the rate at which heat must 
be removed from the plate to maintain the specified surface 
temperature. The middle set of curves in Fig. + show the 
shapes of the mass concentration profiles for VW. = 12 and 
T = 6. The velocity profiles shown in the lower set of curves 
in Fig. 4 were determined from the calculated mass-velocity 
profiles as u/u.. = (4%4¢,)df/dn,. 

Results corresponding to those in Fig. 4, but for the case of 
zero dissipation (designated M., = 0), are given in Fig. 5 for 
T, = 2and 6. Note that the concentration profiles are not 
markedly altered by neglect of the frictional heating. The 
effect of wall temperature on the temperature profiles is shown 
in Fig. 6 for VM. = O and 8. In contrast to the results of 
Smith (8), no reversals in the heat flow (surface temperature 
gradients) were noted in the combination of Mach numbers, 
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Fig. 4 Temperature, mass concentration and velocity profiles in 
a laminar hydrogen-air boundary layer as a function of wall con- 
centration W,, for a free-stream Mach number of 12 and a wall- 
to-free-stream temperature ratio of t/t. = 6 (t.=392 R) 
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Fig.5 Temperature, mass concentration and velocity profiles for 
nondissipative flow 


wall temperatures, and surface concentrations considered. 
Results for the adiabatic-wall boundary condition are il- 
lustrated in Fig. 7 tor M,, = 12. Note the large increase in 
ny. The transport properties at temperature ratios as high 
as T = t/t. = 20 are quite insecure, particularly the mixture 
conductivity. In practice, the independent parameter may be 
considered as fy = — 26, (WV. Re,) since this fixes the 
“coolant’’ flow requirements. The connection between 
and W,,, is given in Fig. 8 for five cases. Evidently the effect 
of Mach number for a given wall temperature is not very 
pronounced. Lower injection rates are required at higher Mach 
numbers to achieve a given surface concentration for the same 
wall temperature. For constant Mach number the injection 
rate required for a given surface concentration decreases as the 
surface temperature increases, the adiabatic wall representing 
a limiting case. 

The local skin-friction coefficients are calculated from the 
definition = c;,(p.U?../2) as Re, = (Gy) un 
(Ce) . These results are shown in Fig. 9 for four 

B w 

Mach numbers and three wall-to-free-stream temperature 
ratios. Again, the effect of Mach number is small for 7, > 1. 
The reduction in skin friction with hydrogen injection is evi- 
dently very pronounced. This decrease is due to the reduc- 
tion in surface velocity gradient resulting from the change in 
profile shape (Figs. 4, 5, and 7) and the increased boundary- 
layer thickness accompanying injection and the diffusion 
process in the binary flow. The skin-friction coefficient for 
zero injection (W. = 0) is c;, = 0.468/ V/ Re, for M o = 12 
and T,, = 6 (7). Hydrogen injection under the same condi- 
tions and at the highest wall concentration in Fig. 9 is seen to 
reduce this value by approximately 90 per cent. 
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Fig. 6 Temperature profiles for various wall-to-free-stream 
temperature ratios 


The local Nusselt numbers required for the calculation of 
heat transfer gq” = h(t; — tw) are derived from the surface 
temperature gradients as Nu/V Re, = (—¢x)w(dT/dnp) w/- 
2(T,. — T,). For M, = 0 the recovery temperature is ¢... 
The recovery temperature 7, = (,/t,, for M., = 12 is shown 
in Fig. 10 as a function of the injection parameter f,. The 
recovery factor may be calculated from these results as ¢ = 
2(T, — 1)/(y~ — 1)M?,. Also shown in Fig. 10 are the local 
Nusselt numbers for three cases. An examination of the re- 
sults for M., = 12 shows that a hydrogen mass concentration 
of 14 at the surface reduces the local Nusselt number by 
nearly 70 per cent of its value for pure air at the surface. The 
hydrogen injection is expected to influence the heat transfer 
in two ways. Although the higher thermal conductivity of 
hydrogen relative to air might be expected to increase the 
heat transfer to the surface at low injection rates, an increase 
in boundary-layer thickness caused by injection and the 
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Fig. 7 Temperature, mass concentration and velocity profiles 
for adiabatic-wall case 


higher thermal capacity of hydrogen tend to markedly de- 
crease the heat transfer. The results in Fig. 10 suggest that 
the last two effects coupled with the lowered density at the 
surface with hydrogen injection strongly overshadow the 
conductivity effect. 

The use of helium injection, although more attractive from 
the point of view of application, is known to be more effective 
than pure air-to-air injection but somewhat less effective than 
hydrogen-to-air injection. This is a reflection of its higher 
density and lower thermal capacity relative to hydrogen.* 
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Fig. 9 Variation of local skin-friction coefficient with surface 
mass concentration 
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The Temperature History in a Thick Skin Subjected to 
Laminar Heating During Entry into the Atmosphere 


GEORGE W. SUTTON! 


General Electric Company, Philadelphia, Pa. 


During high speed entry into the earth’s atmosphere, a 
vehicle can be afforded thermal protection for the short 
period of entry heating by a thick outer skin, sometimes 
called a ‘“theat sink.”” The temperature distribution in 
such a heat sink has been found by integrating the product 
of the laminar aerodynamic heating rate and the appro- 
priate Green’s function for a finite-thickness wall over the 
generalized trajectory for a vehicle entering the earth’s 
atmosphere at high speeds. The result is given in closed 
form, thus solving completely the one dimensional heat 
conduction problem for laminar heating. The maximum 
surface temperature that occurs during the generalized 
entry trajectory for any combination of wall thickness and 
thermal properties is obtained. This maximum tempera- 
ture is given by a single dimensionless curve, from which 
the performance of any material can be found, provided 
that the average thermal properties may be used. As an 
example of the use of the solution, the performance of 
copper, graphite, molybdenum and tungsten are com- 
pared. 


Nomenclature 
a = aconstant, taken as 13/4 
A = reference area of vehicle 
B = aconstant, taken as 20,800 Btu/ft-sec °F 
Cy = hypersonic drag coefficient 
Cy, = specific heat of heat sink material 
C = constant 
G = Green’s function 
g = 32.2 ft/sec? 
H = 21,400 ft 
h = stagnation enthalpy of air, Btu/Ib 
h, = stagnation enthalpy of air at edge of boundary layer, 
Btu/lb 
k, = thermal conductivity of heat-sink material, Btu/ft?-sec 
L = heat sink thickness ; 
q = heat transfer rate, Btu/ft?-sec 
Qr = total heat pulse per unit area, Btu/ft? 
r = dummy of integration 
R = radius of curvature of vehicle at stagnation point, ft 
V = velocity, ft/sec 
V. = sea level satellite velocity, ft/sec 
V. = entry velocity, ft/sec 
7 = temperature, deg F 
t = time, sec 
w = weight of heat sink per square foot of exposed area 
W = weight of entry vehicle, lb 
x = distance from heat sink outer surface, ft 
y = altitude, ft 
Z = Y — Ym ft 
a, = thermal diffusivity of heat sink material, ft®/sec 
B = HpgCpA/2W sin 
y = rVaH/LVV, sin 6 
«, = emissivity of heat sink surfaces 
6 = entry flight path angle, measured from horizontal 
o = Stephan-Boltzmann constant 


= (T = pr sin 0/qnWVH 
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p = density of air, slug/ft* 
py» = sea level density of air, slug/ft* 
pe = 0.00363 slug/ft® 
p, = weight density of heat sink material, lb/ft* 
7 = dimensionless time, (V, sin 6)t/H 
Subscripts 
h = heat sink material 
m = maximum 
Introduction 


HEN a vehicle enters the earth’s atmosphere at high 
speed, it is necessary to protect the vehicle structure from 
the intense aerodynamic heating that exists for a short period 
of time. One possible protection scheme is the use of a thick 
outer skin on that portion of the vehicle that is subjected to 
the intense heating, thus allowing the heat energy to be ab- 
sorbed by the temperature rise of the skin. This thick skin, 
sometimes called a “heat sink,” is thermally insulated from 
the vehicle structure. 
In principle, the maximum energy that such a heat-sink can 
absorb is just 


Tm 
Qr = f LpiC p(T 
To 

per unit area, if the problem of heat conduction is ignored. 
However, this cannot be achieved in practice, because the 
heat energy takes a finite time to be conducted from the heated 
face of the heat sink to the insulated back face; hence during 
the heating cycle, the temperature of the former will be greater 
than that of the latter. This problem was originally con- 
sidered by Wagner (1),? but was limited to two simple cases: 
Constant velocity with a finite wall, and an approximate tra- 
jectory with a semi-infinite wall thicknesses. Mason and 
Gazley (2) have also considered this problem. 

Herein, a method will be developed for determining the 
temperature profile in a heat sink of arbitrary thickness, and 
the maximum surface temperature that will be reached during 
a complete entry trajectory. Numerical results are presented 
in the form of a single curve which relates the maximum wall 
temperature to the maximum heating rate, and which also 
includes the effect of trajectory, material properties and wal! 
thickness. The heat sink capabilities of several materials are 
calculated from this curve. 

The following assumptions were made: The material prop- 
erties are represented by their average values; the wall thick- 
ness is small in comparison to the radius of curvature of the 
outer contour of the heat sink; and the meridional variations 
of heat sink thickness and heating rate are small, so that the 
heat flow is essentially one-dimensional from the heated front 
face to the insulated backface. 


Entry Trajectory 


During the entry into the earth’s atmosphere of a high spee:| 
vehicle, the deceleration is large in comparison to gravity. 


2 Numbers in parentheses indicate References at end of paper. 
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Hence gravity can be neglected. With this assumption, the 
flight path angle during entry will remain constant while the 
actual trajectory is slightly curved. Also, the calculated 
values of time will have an error which increases with de- 
creasing altitude. However, at the time of maximum wall 
temperature, which is the region of interest, these errors are 
small and are within the accuracy of analysis. The equation 
of motion of the vehicle along the flight path is then 
W dV 


Since the flight path angle is taken constant 


dV dV dy —sinOd 
dt dy dt 2 dy 


The density can be approximated by an exponential, that is 


For p. = 0.00363 slugs/ft® and H = 21,400 ft, the maximum 
error in p is six per cent compared to the standard atmosphere, 
but this is well within the experimental error of density de- 
termination at high altitudes, for altitudes up to 200,000 ft. 
Since at hypersonic velocities Cp is a constant, Equation 
[1] can be integrated to obtain the vehicle velocity V as a 
function of altitude by use of Equations [2 and 3] as follows 


[4] 
where 
[5] 
oWend | 
To obtain time as a function of altitude, note that V = 
—(dy/dt)/sin 6 and integrate [4] to obtain 
H —y Be —2u/H 
t= —y/H 23° 
3p —3y/H y 
3°3! ym 


where t = 0 is chosen to correspond to that altitude y,, at 
which the convective heat transfer isa maximum. Equations 
|4] and [6] are derived in (1); however, later in this analysis, 
it will be convenient to neglect terms in the series higher than 
Bexp(—y/H). At the altitudes of greatest heating, the error 
so introduced is negligible; at lower altitudes, where the heat- 
ing rate is 10 per cent of maximum, the error in time is less 
than 20 per cent. The net effect of this assumption is negligi- 
ble. 


Heat Transfer Rate 


In the vicinity of the stagnation point of a blunt-nosed hy- 
personic vehicle the boundary layer will be laminar, in which 
case the convective heat transfer rate at the stagnation point 
may be approximated by (3) 


Ry hy 
R po V, 


Equation [7], when compared to the results of a similarity 
solution of the laminar boundary layer equations, is accurate 
to within five per cent for the following values of the con- 
stants: B = 20,800 btu/ft?sec, Ro = 1 ft, and a = 13/4, ac- 
cording to (3). At points on the body other than at the stag- 
nation point, the heating rate at large Mach numbers is given 
by 


= ‘ [8] 
Methods for determining C, which is primarily a function 
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of geometry, are contained in (4). However, in the following 
analysis, the heating cycle to a point on the body will be 
normalized to the maximum heating rate at that point and, 
hence Equation [8] will have no effect upon the results. 
Thermal! radiation from the shock layer to the vehicle will be 
neglected; this is permissible if the vehicle skin is plated with 
a material of low absorptivity, or if the convective heat 
transfer is much larger than the shock layer radiation. 

In the relations to be derived, it is assumed that h../h. ~ 0; 
the results, however, are applicable for cases where €, and 
h./h, are not negligible, by calculating q,, as follows 


Ry p ‘( he 


This approximation is valid when the correction to qm is small; 
this is true because the heating cycle including radiation from 
the surface, when normalized against ¢,, as given by Equation 
[9], is virtually the same as the heating cycle given by Equa- 
tion [7], when h,./h, ~ 0. Hence q» is obtained by iteration 
on the value of 7',, for which an expression is derived later. 


Heating Cycle 


When Equations [3 and 4] are substituted into [7], the 
heating rate as a function of altitude is obtained 


e=B Pe (7) e-(y/2H+ape 
R po Wa 


By differentiating [10] with respect to y, the altitude, y,,, at 
which the heating rate is a maximum is obtained 


Ym = Hin(2aB8)............. 


which corresponds to a velocity of 


If a = 13/4, then V = 0.858V,... Thus the maximum heating 
rate oecurs before the vehicle has slowed down appreciably. 
The maximum heating rate is obtained from the substitution 
of Equation [11 and 5] into [10] 


(® W Sin 

Qn = B 
po ageH 

To simplify the equations, normalize [10] by dividing it by 


[13] and by measuring altitude from the altitude of maximum 
heating; that is, let 


The following expression for the normalized heat transfer is 
then obtained 


/H—1+e7 2/4) 


Hence the normalized heat transfer rates in terms of Z are 
identical and are independent of the value of a. This same 
result has been obtained (5) for a = 3.0. It is instructive to 
determine whether a similar relation exists for time. The sub- 
stitution of Equations [11 and 14] into [6] vields 


V. sin Z 1—e% H 
H H 2a 
_.. . [16 
2-2'(2a)? 3-3!(2a)3 


Thus, time scales proportionally to the vertical component of 
the entry velocity, where for convenience a dimensionless time 
is defined as 
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Fig. 1 Convection heating rate during entry 


NORMALIZED HEATING RATE 


° 


Equations [15] and [16] are a parametric representation of 
q/qm vs. r, Which is independent of missile weight, drag, or 
trajectory. This is shown in Fig. 1, fora = 13/4. 

The ratio of the total heat transfer to the maximum heat 
transfer rate can be obtained from an integration of Equation 


[15] 


Note that dt/dZ = —1/V sin @, and substitute Equations [4 
and 5] into [16] to obtain 
% 
Qr (2) -erf[B(a — 1)]%4..... [19] 


qm V. sin 6 


Since 6 > 1, erf [B(a — 1)]% =~ 1. Asan example, a vehicle 
with an entry velocity of 36,000 fps, and flight path angle of 
80 deg has a ratio of Qr/qm of 3.60. The following explicit 
expression for Q,r is obtained by substitution of Equation [13] 
into [19] 


B 
V.\a —1 sin@ pp CpA gR 


Maximum Temperature for Back Radiation 


Before proceeding to the calculation of the heat conduction 
into the skin of the vehicle, one may inquire as to whether 
it is possible to radiate all of the convective heating from 
the vehicle skin. As a first approximation, the heat capacity 
of the skin is neglected. Then q, as given by Equation [13] 
may be equated to oe,7’.*. This. procedure neglects hu/he; 
thus the calculated values of 7’, will be somewhat too high. 
The resulting dependence of 7’, on W/CpAR and entry angle 
is shown in Fig. 2. It may be seen that a small entry angle 
which decreases the convective heating, and a refractory ma- 
terial which increases the radiation, are highly desirable. Note 
that relatively heavy entry bodies with large entry bodies can 
back-radiate the convective heat transfer if the outer shell is 
made of a refractory material with very high temperature sta- 
bility. 

For a given body and trajectory one may obtain a more 
precise value of 7’. by correcting the convective heating with 


the factor 1 — as, see Equation [9]. The value of h, may be 


obtained from normal shock calculations (6) for the value of 
V, and y» of that trajectory; the value of h. may be ob- 
tained from the value of the stagnation pressure behind the 
shock and, as a first approximation, the value of 7, as ob- 
tained from Fig. 2. The corrected value of gm is then equated 
to and solved for T'.. 

Back radiation appears to be practical only for light bodies 
with shallow entry angles. In the remainder of this paper, 
back radiation will be neglected, but the thermal capacity of 
the skin will be included; this corresponds to large entry 


42 


198 T T T T 
he 
\ OK B000°R 
102 1 
8 + 
+ 
4 
6000°R 
10! 
6 
\ 
5000°R 
2 
10° 
\ 
2 
3000°R 
\ 
4 
2 
10. +20 40 50 60 70 60 


30 
©, ENTRY ANGLE, OEG. 


Fig. 2 Surface temperature for back radiation 


0.10 


0 04 08 12 
SIN@ 
e 


Fig. 3 Effect of trajectory and heat sink properties on heat sink 
maximum temperature 
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angles, high entry velocities, and temperature rises which are 
much smaller than that required for back radiation. As an 
introduction to the method, two limiting cases are derived: 
The temperature rise of a very thin wall, and the temperature 
rise of a semi-infinite wall. 


Maximum Temperature of a Thin Wall 


The calculation for the temperature of a thin walled heat 
sink is simple, for as L > 0, T(z, t) > T.., hence 


Substitution of Equation [19] into [20] vields 
mH 2 
Lp. \a 1 


It is convenient to define a new independent variable, which 
reflects the wall thickness, as 


arg [22] 


and a new dependent variable which is a measure of the tem- 
perature rise of the wall as 


T(Z) — Ty 
qm 
Then equation [21] becomes 
2ae % 
(a —1 


for a thin wall. Equation [24] is shown in Fig. 3, as “hk > 


A thin wall usually does not have sufficient heat capacity to 
absorb the aerodynamic heating without surface melting. To 
prevent surface melting, the heat sink must be sufficiently 
thick, in which ease the thin-wall approximation is in- 
adequate. The limiting case occurs when the back face ex- 
periences zero temperature increase; this corresponds to 
semi-infinite thickness, which is considered next. 


Maximum Surface Temperature of a Semi- 
Infinite Wall 


The heat sink that can absorb the aerodynamic heating with 
the least surface temperature rise is sufficiently thick so that 
the back insulated side experiences little or no temperature 
rise. Hence, for the purpose of mathematical analysis, it is 
permissible to let L ~ ©; then calculated the temperature 
response of this wall and also the value of x for which the tem- 
perature rise is negligible when the surface temperature 
passes through its maximum. This value of z is called the 
semi-infinite thickness. 

The surface temperature is calculated as follows: consider 
a unit amount of heat transferred to the wall at 2 = 0 at time 
t;; the temperature distribution at some later time ¢ is given 
by the following Green’s function (7) 


4) = — 5 
V ky th) 

Since the heat conduction equation is linear, the temperature 
Tat x = 0 due to a heating rate q(t) is given by 


] 


The factor ¢ — t, is obtained from Silden [6 and 11] by 
taking only the first two terms as follows 
V. sin 0 E H 2a (e (27) 
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= ¥ — Ym, and ¢ is the altitude corresponding to 4. 


where Z 


By use of the relationship dy = —V sin 6 and Equation [4], 


dt 
dt, becomes 
jae 

and from Equation [15] 
Q(t) = [29] 

To simplify notation, let 


Substitution of Equations [27-30] into Equation [26] and the 
definition of [23] yield 


/H) 


e 


where r = ¢/H is the dummy of integration. 

The independent variable y does not appear in Equation 
[31] because y ~ 0 as L ~ ©, The integral in Equation 
[31] has been evaluated by means of a digital computer. 
The maximum value @¢,, is 2.06 for a = 13/4, and occurs when 
Z/H = —0.80. This value of y = 0 gives the maximum 
amount of ene rgy that a given waahieti al can absorb for a given 
temperature rise. Equation [381] is shown in Fig. 3 as “IL —> 

Since it is usually desirable to have as thin a heat sink as 
possible in order to conserve weight, the wall thickness may 
be less than the semi-infinite thickness. The limitation on L 
depends on the maximum allowable temperature rise. In 
the next section, this dependence is derived. 
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Fig.4 Heat sink capacities of several materials 
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Maximum Temperature of a Heat Sink of 
Arbitrary Thickness 


Although it is possible to correct Equation [31] for the case 
of a finite wall thickness by considering “reflections” of the 
heat at the insulated back face, the problem can be solved 
directly by using the Green’s function for a wall of thickness 
L. It is easily shown that the Green’s function is given by 


G(z,t —t) = - J cos sailed | 
piCprL L 


n=1 


The temperature rise is obtained in a similar manner to 
Equation [31]. By the use of the definitions of @ and y and 
Equations [22, 23], this becomes 


P(x, 7, Z) 
T(z, Z) — J 2ae 
H ma — 1) 


2a 


0 


n=1 


1 -r 
cos (=) e dr. . [33] 


Because of the complex way in which y and Z appear in 
Equation [33], there exists no generalized temperature curve 
which isolates the effect of y. Hence $(Z) was calculated at 
x = 0 with a digital computer for various values of y, with a 
= 13/4; and the maximum values @¢,,(7) were obtained. 

The results are shown in Fig. 3. It is seen that, if y < 0.75, 
the results correspond to an “infinite”? heat sink; while if 
7 > 2.0, the heat conduction through the wall can be neglected 
and only the heat capacity of the heat sink need be considered. 

From Fig. 3, the maximum heat sink temperature for any 
thickness of any material can be found, provided that the 
thermal properties of that material can be represented by 
temperature-averaged values. First the value of y is calcu- 
lated and from Fig. 3 the value of ¢, is obtained. Then, by 
the use of the value of q¢,, for that trajectory, T, — To is ob- 
tained by means of Equation [23]. On the other hand, if the 
maximum temperature rise is specified, the required heat sink 
thickness L can be obtained. A practical upper limit to the 
thickness corresponds to y = 0.75. 

To illustrate more fully the use of Fig. 3, the heat absorp- 
tion capacities of four different materials have been obtained 
and are shown in Fig. 4. The weight per square foot 
has been obtained by multiplying L by the weight den- 
sity of that material. These calculations neglect the effect 
of surface combustion; this may be permissible if the con- 
vective heating rate is much larger than the rate of energy re- 
lease due to combustion. The thermal properties of the four 
materials are shown in Table I. 

It may be seen from Fig. 4 that tungsten is capable of ab- 
sorbing the greatest amount of heat. For maximum heating 
rates which are somewhat less, a shell fabricated of graphite 
will weigh the least. Copper is limited to a maximum heating 
rate of 1818 [V, sin 6/V,]% Btu/ft?-sec. The calculated ca- 
pacity of copper agrees to within 5 per cent of a Dusinberre 
numerical calculation of the temperature rise, using an actual 
trajectory. 

The altitudes corresponding to the maximum heat sink 
surface temperatures are shown in Fig. 5. Since Z = y — ym, 
the maximum temperature always occurs after the maximum 
heat transfer. 

For an actual trajectory, 7, — To and q» are interde- 
pendent but the following iteration scheme may be used to 
find the temperature increase: The zero approximation of q,, 
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Fig. 5 Altitude of maximum surface temperature 


is calculated, assuming 7',, — Ty) = 0. With q,,.(0), the zero 
approximation for 7,,,(0) is obtained from Fig. 3. Since the 
T »(Z) curve is flat in the vicinity of Z = 0, assume that the 
surface temperature at Z = 0 is 7,,(0). Then q(1) is ealeu- 
lated from Equation [9], and the first approximation T’,,(1) is 
again obtained from Fig. 2. This process converges rapidly 
because of the small dependence of 4», on T',. 


Temperature Profiles 


The temperature profiles at various times in a wall for which 
y = 1.245 are shown in Fig. 6. These were obtained from 
Equation [33]. They illustrate the fact that the back face of 


Table I Selected Physical Properties 

| Temperature 

at which k 

evaluated Btu/ft p 
Material deg F sec F lb/ft? Btu/Ib 
Copper? 1112 0.0567 545 0.1057 
Graphite? 1472 0.0106 105 0.344 
Tungsten’ 2192 0.0162 1204 0.0377 
Molybdenum! 2192 0.0189 638 0.075 


'W.H. McAdams. ‘Heat Transmission,’’ 3rd Edition, 
McGraw-Hill, 1954. 

* Courtesy V. C. Hamister, National Carbon Company. 

8]. C. T. and “Handbook of Chemistry and Physies,”’ 
33rd Edition, p. 1875. 

4 “Metals Handbook,’’ 1948. 
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the heat sink reaches temperatures almost as high as the 
heated side. Note that the value of ¢ is 2.30, which is only 10 
per cent greater than the infinite-wall value of 2.06 and yet the 
infinite-wall approximation cannot be used to predict the 
temperature profiles correctly. 


723.55 


T=0.95 


Fig. 6 Heat sink temperature profiles at various times for 
y = 1.245 
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Conclusions 


By use of the Green’s function for heat conduction, the 
temperature distribution in a heat sink has been found in 
closed form, thus solving completely the one-dimensional 
heat conduction problem for a laminar boundary layer. The 
values obtained are five per cent lower than a Dusinberre 
calculation for an exact trajectory. 

The maximum surface temperature has been calculated in 
dimensionless form for any trajectory. A calculation of a 
temperature profile history shows the inadequacy of the semi- 
infinite slab solution. The heat sink capacities, which are 
easily found from the analysis, have been calculated for 
several different materials. 

The accuracy of the solution depends on the accuracy of the 
expression for the heat transfer rate. The solution is ap- 
plicable to a laminar boundary layer and one-dimensional heat 
conduction. It is possible to extend the method for the case 
of the turbulent boundary layer and for axially-symmetrie 
heat conduction. 
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Factors Which Influence the Suitability of Liguid 


Propellants as Rocket Motor Regenerative Coolants 


D. R. BARTZ! 


Jet Propulsion Laboratory, California Institute of Technology, Pasadena, Calif. 


A brief review of methods of utilizing liquid propellants 
as rocket-motor coolants is given, and criteria for accept- 
able cooling are established. Methods of analytical and 
experimental determination of the pertinent heat-transfer 
characteristics of propellants are described. The measured 
maximum rates at which liquid ammonia can satisfacto- 
rily accept heat are presented to show typical trends with 
variations of fluid pressure, temperature and velocity. 
Measured maximum rates for several propellants are com- 
pared with predictions of heat fluxes from the combustion 
gases of these propellants when producing the same thrust 
in a hypothetical rocket motor. From this comparison, an 
insight into the importance of various propellant proper- 
ties and operating conditions is gained. 


Nomenclature 
c* = characteristic velocity, ft/sec 
Dep = deposition 
FC = forced convection without boiling 
LS = specific impulse, sec 
L* = characteristic length, in. 
NB = nucleate boiling 
p = pressure 
q = heat flux, Btu/in.? sec 
dul = heat flux at the upper limit of nucleate boiling, Btu/ 
in.? sec 
r = mixture ratio, oxidizer to fuel by weight 
t = height of coolant-passage annulus, in. 
T = temperature 
V = coolant velocity, ft/sec 
W = coolant flow rate, lb/sec 
Ap = coolant-passage pressure drop 
Ap,-0 = coolant-passage pressure drop without heat transfer 
Hw/w, = ratio of viscosity at the wall to that in coolant-passage 
free stream 

Subscripts 

ce = chamber 

L = bulk 

out = coolant-passage outlet 

sat = saturation 

w = wall 

* = at the nozzle throat 

Introduction 


A MAJOR problem confronting the designers of long- 
duration liquid-propellant rocket systems is designing 
motor walls compatible with the extreme temperatures result- 
ing from the combustion. Three methods for achieving the 
desired resistance to the high temperatures are commonly 
considered. They are refractory liners, regenerative cooling 
and film or transpiration cooling. Of the three, regenerative 
cooling is most commonly utilized in propulsion systems cur- 
rently in production or development. Regenerative cooling is 
accomplished by flowing one or both propellants through 

Presented at the ARS Spring Meeting, Washington, D. C., 


April 4-6, 1957. 
' Research Group Supervisor, Power Plant Research Section. 


46 


motor-wall coolant passages on their way from the propellant 
tanks to the injector. The missile performance penalty of re- 
generative cooling is the additional turbo-pump weight an:| 
power consumption or additional gas pressurization require 
to overcome the coolant passage pressure drop. In order to 
select the method best suited for a contemplated design, this 
penalty must be compared with the additional dead weight 
and fabricational difficulties associated with refractory liners 
and with the penalty of additional propellant consumption 
associated with film or transpiration cooling. The limitations 
of each of these methods are such that for extremely high- 
performance motors, optimal combinations of all three 
methods may be necessary. This paper presents and dis- 
cusses the factors which determine the suitability of a given 
propellant as the coolant for a completely regeneratively 
cooled motor. 

The criterion for adequate cooling is simply that the 
coolant flow must maintain the walls of the motor at tem- 
peratures below which failure might occur due to melting 
and/or stress. For metals commonly used for rocket motor 
walls, such as stainless steel, nickel, inconel, or mild steel, the 
limiting combustion side surface temperature is in the region 
from 1500 to 2000 F. Thus, with combustion temperatures 
from 3000 to 8000 F, resultant differences between the free- 
stream gas temperature and the wall temperature are from 
1500 to 6000 F. The consequence of such high driving po- 
tentials are heat fluxes to the wall, which, in units of Btu/in.* 
sec, range from about 1 to 20. Such fluxes are orders of mag- 
nitude higher than those encountered in conventional heat 
exchange equipment and even significantly higher than en- 
countered in most nuclear reactors. 

These high levels of heat flux in conjunction with conditions 
of coolant velocity, temperature, and pressure are such that, 
over most of the cooled surface, heat is transferred by boiling 
when the pressure of the system is below the critical pressure 
of the coolant. Some coolants, however, are necessarily 
operated at pressures above their critical pressure; hence, 
both the boiling and supercritical heat-transfer characteristics 
of a wide variety of liquids are of interest in connection with 
rocket-motor regenerative-cooling applications. A qualitative 
description of the characteristics of the various pertinent re- 
gimes of heat transfer and the methods utilized for theoretical 
and experimental evaluation of these characteristics will serve 
to define terms and to introduce the problems involved. Such 
descriptions are best illustrated by a plot of heat flux as a 
function of wall temperature for constant pressure bulk tem- 
perature, and velocity, Fig. 1. 

Considering first a liquid at a pressure below its critical 
pressure, and a bulk temperature well below its saturation 
temperature, the line segment A-B represents the range of heat 
fluxes which can be sustained without boiling. The wall tem- 
perature at B exceeds the saturation temperature by any- 
where from 10 to 100 F, which constitutes superheating of the 
liquid. At heat fluxes slightly higher than B, a few bubbles 
are formed at favored nucleation spots on the hot wall, and 
they grow rapidly out into the cold-liquid stream. They 
continue to grow until condensation at the cold-liquid surface 
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begins to exceed the rate of vaporization at the base of the 
bubble, whereupon the bubbles begin to collapse. The inertia 
of the fluid following the collapsing bubble is responsible for 
completely collapsing the bubble. This process, which occurs 
at exceedingly high frequencies, is called nucleate or surface 
boiling. As still higher heat fluxes are encountered in the re- 
gion B-C, an increasing number of nucleation sites are 
brought into play, thus increasing the bubble population. 
This is accomplished with only slight increases in wall tem- 
perature. The resulting increase in heat-transfer coefficient 
represented by the near vertical part of the curve (B-C) is be- 
lieved by many, although not unanimously (1),? to be the 
result of increasing turbulent mixing or stirring created by the 
growth and collapse of the bubbles. Unfortunately, this 
favorable situation has a limit. At this limit C, a slight in- 
crease in the heat flux leads to such a dense bubble population 
that the bubbles coalesce into a vapor film with an attendant 
large decrease in heat transfer coefficient. The resulting shift 
in wall temperature, C-D, is for many liquids so high that 
failure of commonly employed metal surfaces occurs. Hence, 
this transition from nucleate to film boiling has been used 
synonymously with the term burnout. However, for many 
other liquids, notably light hydrocarbons and alcohols, the 
temperature at D is not much above 1000 F, and failure does 
not occur. Thus, the term burnout becomes misleading. To 
avoid this difficulty, point C shall be referred to as the upper 
limit of nucleate boiling with the symbol q,,. At heat fluxes 
above C, even for liquids which do not cause failure, the wall 
temperature increases rapidly with increases in heat flux 
(D-E) so that for practical purposes the value of the heat flux 
at the upper limit of nucleate boiling must be used as the de- 
sign limit for regenerative cooling systems. 

The picture for supercritical pressure liquids is much 
simpler. Since no boiling can occur, the wall temperature in- 
creases monotonically with increasing heat flux (F-G), the 
heat-transfer coefficient remaining essentially constant until 
the wall temperature reaches the critical temperature. At 
higher heat fluxes (G-H), the heat is transferred through a 
gaseous layer adjacent to the hot wall, resulting in somewhat 
lower heat-transfer coefficients. Failure temperatures are 
usually reached at much lower heat fluxes when a liquid is 
operated above its critical pressure than when operated below 
its critical pressure. 


Determination of Heat-Transfer Characteristics 


Analytical Methods 


The relationship between heat flux and wall temperature, 
i.e., the heat-transfer coefficients, in the nonboiling, subcritical 
temperature regions of both subcritical and supercritical pres- 
sure liquids (A-B and F-G of Fig. 1) can be predicted with 
sufficient accuracy for design purposes by the Sieder-Tate 
equation (2). The data for a wide variety of liquid propel- 
lants have been well-correlated by this equation. It has suc- 
cessfully predicted the effects of pressure, velocity and bulk 
temperature. The principal difficulty in applying this equa- 
tion to new propellants is the need of values of the transport 
properties, viscosity and thermal conductivity. They gen- 
erally cannot be predicted from theory. While measurement 
of viscosity is a rather simple laboratory procedure, measure- 
ment of thermal conductivity of liquids is considerably more 
difficult. 

The wall temperature at which boiling is initiated (B of Fig. 
1) can usually be predicted to within about 50 F if the relation 
between saturation temperature and pressure is known for the 
liquid. The amount of superheat for most liquids is from 10 
to 100 F, the higher values corresponding to higher liquid 
velocities, although no general correlation is known. For de- 
sign purposes, it is sufficiently accurate to assume the wall 


? Numbers in parentheses indicate References at end of paper. 
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P= Peritical = 
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x CONSTANT VELOCITY, ABOUT 30 ft/sec 
z | a CONSTANT BULK TEMPERATURE, ABOUT 100° F 
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Fig. 1 Heat flux vs. wall temperature of typical propellant in 
various heat transfer regimes 


temperature at the inception of nucleate boiling is 50 F above 
saturation temperature. It should be noted that bulk tem- 
perature has only a second-order effect on superheat. Satura- 
tion temperature curves for most pure liquids are given in the 
literature. Low-pressure data are best extrapolated by plot- 
ting log p vs. reciprocal of absolute temperature. For many 
multicomponent liquid mixtures, especially those which de- 
compose at elevated temperatures, the saturation tempera- 
tures (or more properly the bubble-point temperatures) are 
not known and, hence, must be measured. 

The slope of the heat flux vs. wall-temperature curve in 
nucleate boiling (B-C of Fig. 1) is not correlated for various 
liquids except that it is generally steeper than the third power 
of the difference between wall temperature and saturation 
temperature. Fortunately, this defines the wall temperatures 
within about 100 F over the range of heat fluxes covered by 
nucleate boiling of most liquids. It should be noted that the 
position of the boiling line is independent of either velocity 
or bulk temperature except through second-order effects on 
superheat which are probably less than 20 or 30 F in most 
instances. 

The determination of the characteristic curve up into nu- 
cleate boiling is thus possible from the considerations men- 
tioned, but this is not of much utility in predicting the 
limitations of the liquid as a regenerative coolant. This is be- 
cause the resulting wall temperatures in nucleate boiling are, 
for most liquids, well below limiting temperatures for suc- 
cessful cooling. The real limitation from a practical view- 
point is the heat flux at the upper limit of nucleate boiling q,; 
since at heat fluxes above this level, the wall temperatures are 
generally beyond the limit for successful cooling. Unfor- 
tunately, this limiting heat flux cannot be predicted as yet 
from knowledge of a liquid’s physical or transport properties. 
This is still the unfulfilled goal of a great deal of boiling-heat- 
transfer research throughout the world. Definite progress has 
been made toward this goal in recent years through con- 
siderations of the dynamics of bubble growth and collapse 
(1, 3, 4). 

The shift in wall temperature (C-D of Fig. 1), the heat- 
transfer coefficients in film boiling (D-E of Fig. 1), and the 
heat-transfer coefficients in the region near the critical tem- 
perature (@ of Fig. 1) are also unpredictable from fluid proper- 
ties with reasonable accuracy. Prediction of the first two 
characteristics is not particularly essential if ¢,; is considered 
to be the limitation. However, predictions of the supercritical 
heat-transfer coefficients are essential since they determine 
the design limitation directly for some coolants. Such pre- 
dictions are particularly difficult near the critical temperature 
because of unusual variations of fluid properties with tem- 
perature. 


' 
eat 
ons 
lat, 
i 
ing 
ure 
rily 
ies 
= 
i 
ith 
cal 
ch 
m- 
ral 
at 
he 
eS 
id 
N 
47 


Experimental Methods 


number of experimental programs have been 


various organizations on a variety of liquid propellants. An 
index of references to these data, including ranges of variables 
investigated and types of data reported, is given in Table 1. 


In order to determine the heat-transfer characteristics of 
propellants that cannot be predicted by analytical means, a 


Most of these data were obtained by flowing the propellant 
through electrically heated tubes or annuli at velocities, pres- 


conducted by 


Table 1 Index to measurements of heat transfer to propellants with and without boiling under forced velocity, 


subcooled conditions 


Pressure Bulk tem- 


Velocity 
range, 
ft/sec 

5-100 
15-43 
30 
30-90 
5-35 
10-130 
5-60 
10-50 


15-60 
(3-40)* 


30 


30 


30-40 


3-40 


(0.6-1.9)* 


6-90 


Data* 


FC, NB, 
qui, Dep. 
FC,NB,qui, 
AP, Dep. 
FC, ‘NB, 
quit, Dep. 
FC, NB, qui 


FC, NB, 
Dep., AP 
FC, NB, 
Dep., AP 
FC, NB, 
qul 

FC, NB, 
qul 

FC, NB, qui 
FC, sub- 
critical, 
super- 

critical 


FC, NB, qui 


FC, NB, 
qul, AP 


FC, 
qul 


NB, 


FC, NB, qui 


FC, NB, qui 


FC, 
Dep. 

FC, NB, 
AP 

FC, NB, 
qut, Dep., 
AP 

FC, NB, qui 


NB, 


FC, super- 
critical 


FC, NB, qui 


Propellant - Organiza-_ Ref- range, perature 
Designation Nominal composition Authors tion rence psi range, °F 
Oxidizers: 
Type III-A 83% HNOs, 14% NO», 2.5% Dean Bell Air- 5 100-950  —45-170 
RFNA H.O, 0.6% HF craft 
RFNA 921/.% HNOs, 6'/2% { 50-500 90-304 
SFNA 83% HNOs, 14% NOs, 2.5% } Bart, Noel 60-500 80-235 
0.6% HF | | 
WFNA 97% HNOs, 1% NO», 2% H20} | 300-900 96-247 
Type III 82% HNOs;, 15% NOs, 3%) | 90-615 25-95 
RFNA H.O, 0.59% HF Wolfe, Gray, Purdue 15-17 
WFNA 98% HNOs, 0.5% NO», 1.5%) Reese, Gra- | 65-710 47-160 
ham 
WFNA 97.5% HNO;, 0.5% NO», Ashley Bell Air- 18 65-400 60-90 
2.0% craft 
RFNA 921/.% HNO;, 6'/2% NOs, Hatcher, Bartz JPL 19 50-500 70-90 
1% HO 
CIF; Commercial grade Noel JPL 20 150-860 0-200 
Baner, Ball Michigan 
rine 
Fuels: 
Liquid Am- 99% NH3;, 1% H2O Noel JPL 23-28 100-1600 50-150 
monia 
Corporal fuel 46% An, 46% FA, 7% N2Hs, Noel JPL 29 = 200-400 
1% HO 
80% An, 20% FA ) 
50% An, 50% FA 
33% An, 67% FA 
47.5% An, 47.5% FA, 5% ) Noel JPL 6 400 125 
40% An, 40% FA, 20% H2N, 
Diacetone Commercial grade Noel JPL 6 400 125 
Alcohol 
Isopropyl Commercial grade Noel JPL 6 400 125 
Alcohol 
DETA (Di- 90.4% diethylenetriamine Noel JPL 11,30 400 10-350 
ethylenetri- 9.6% amino-ethy]-piperazine 
amine) 
Aniline * 96% aniline, 4% HO Kreith, Sum- JPL 31 35-365 90-135 
merfield 
n-Butyl 99.5% n-butyl alcohol, 0.5% Kreith, Sum- JPL 32 50-200 85-100 
Alcohol H,O merfield 
JP3 Mil-F-5627 Hatcher JPL 19,33 50-500 80-150 
JP4 Mil-F-5624A Beighley,Dean Bell Air- 34 30-500 75 
craft 
[Eee Powell JPL 21 750 —390-700 
Hydrogen 
Hydrazine __............ Banchero, Univ. of 22 
Baner, Ball Michigan 
Hydrazine 97% NeHy, 3% HO Noel JPL 35 300-1200 70-320 
* Units: Ib/in.? sec; FC—nonboiling forced convection, NB—forced velocity nucleate boiling; 


boiling; AP = pressure drop; Dep. = deposition or coking; An = aniline; FA = furfuryl alcohol. 
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sures, and bulk temperatures common to rocket-motor cooling 
passages. Generation of heat by use of the tube wall as a re- 
sistance element was found to be the most practical means 
of achieving and controlling the level of heat fluxes require. 
Measurements required for determination of the desired heat- 
transfer characteristics are power consumption for comput:- 


qur—upper limit of nucleate 
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tion of heat flux, flow rate, inlet and outlet pressures and 
temperatures, and tube-wall temperatures. Because the hy- 
draulic diameters of most rocket motor coolant passages, even 
up to very large scale, are between } and 3 in, tests made 
with tubes and annuli in this range of sizes suffer little in 
scaling. Fortunately, from an experimental standpoint, it is 
not necessary to scale channel length for boiling heat-transfer 
measurements since uniform heat-transfer conditions are 
established very close to the entrance of a forced-velocity 
boiling system. Tubes and annuli of 10 to 20 diam in length 
are commonly used for such measurements. A necessary but 
complicating part of some of these experimental investiga- 
tions has been measurements of such adverse effects as de- 
position of solids on the heated wall, decomposition of the 
liquid, and corrosion. 


Typical Results of Experimental Investigations 


Since there is rather complete experimental confirmation of 
the methods described for predicting the nonboiling curve, the 
inception of nucleate boiling, and the wall temperatures in 
nucleate boiling, the most significant results of experimental 
investigations are values of qu, as a function of the three 
s\stem variables: Pressure, bulk temperature and velocity. 
Although, as mentioned, no general correlation of this quan- 
tity has been achieved, certain trends have been fairly well 
established. These trends are well represented by data re- 
cently obtained by Noel for anhydrous liquid ammonia 
(23-25). 

1. Effect of pressure on q,,: It has been found that maxi- 
mum values'of q,, are obtained at pressures which are about 
0.4 to 0.6 of the liquid critical pressure for given constant bulk 
temperatures and velocities. The curve usually slopes quite 
gently on both sides of the peak. The maximum value for q,; 
for ammonia at 60 F was obtained at a pressure of about 850 
psi or about 52 per cent of critical pressure (Fig. 2). 

2. Effect of bulk temperature on q,,: It has been found that 
values of q,, decrease with increasing bulk temperature, reach- 
ing a near-zero value at saturation temperature. (The only 
exceptions to this trend have been observed at bulk tempera- 
tures less than about 40 F, where it is believed that viscous 
effects begin to control the boiling mechanism; (5, 26-29). 
For some liquids, the decrease is nearly linear with bulk tem- 
perature; for others, the values of q,, are greater than for a 
linear decrease. The data for ammonia (Fig. 3) exhibit this 
latter trend. Several correlations, notably for water (36, 37) 
have successfully made use of a correlation variable called 
subcooling, which is the difference between saturation tem- 
perature and bulk temperature. For such correlations, ¢,,, de- 
creases linearly or less than linearly with decreasing sub- 
cooling. 

3. Effect of velocity on qu: The q,, data for ammonia plot- 
ted for constant pressure and various bulk temperatures in 
Fig. 4 show trends with increasing velocity, which are typical 
of several propellants investigated. The values of ¢,, increase 
linearly with velocity from zero-velocity values which are 
dependent on pressure and bulk temperature. The slope of 
the curve is found to be dependent on subcooling with lower 
slopes at lower subcooling. 

4. Effect of deposition on q,,: Deposition of carbon by JP3 
(33) and deposition of metallic salts by fuming nitric acids 
(5-14) have been observed to cause rather rapid increases in 
wall temperatures both with and without nucleate boiling. 
The effect of such deposition on q,, from these tests is not 
clear, but the general conclusion drawn from the data of 
(6-14) is that it did not particularly affect the values of q,, 
except in the case of SFNA. 


Dependence of q,; on System Operating 
Conditions 
Having previously discussed the reasons for considering 4,,; 
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Fig. 2 Heat flux at upper limit of nucleate boiling as a function 
of pressure for ammonia at bulk temperature 60 F +5 
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Fig. 3 Heat flux at the upper limit of nucleate boiling as a func- 
tion of bulk temperature for ammonia at pressure 500 psi veloc- 
ities of 30 and 60 ft/sec 
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Fig. 4 Heat flux at upper limit of nucleate boiling as a function 
of velocity for ammonia at pressure 500 psi 
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as the practical limiting heat flux and having shown the man- 
ner in which g,; depends on pressure, bulk temperature, and 
velocity, it is now of interest to show the effect of the rocket- 
motor operating conditions on the regenerative cooling prob- 
lem. 

The pressure at the outlet of the coolant system is estab- 
lished by the combustion pressure and the injector pressure 
drop. It is not generally economical to operate the coolant 
passage at any higher pressure than necessary. It was noted 
that pressure affected the value of g,,; only slightly over the 
pressure range from about 0.3 to 0.7 of critical pressure. On 
the other hand, the net heat flux from the combustion gas to 
the motor wall is expected to increase with pressure to about 
the 0.8 power (38) the effect of pressure on nucleate boiling 
wall temperatures being of second order in determining the 
heat-transfer driving potential, (7, — 

The coolant bulk temperature at any point in the coolant 
passage is simply the sum of the inlet bulk temperature and 
the temperature rise to that point. The temperature rise is 
determined by the heat transfer to the wall, the coolant flow 
rate, and coolant specific heat. Since the higher values of q,,, 
are associated with the lower values of bulk temperature, it is 
desirable to utilize the largest coolant flow rate possible. For 
a fixed total propellant consumption, operation at high mix- 
ture ratios (0/F) favors cooling with the oxidizer, and opera- 
tion at low mixture ratios favors cooling with the fuel. Since 
both propellant consumption and heat flux increases nearly 
linearly with chamber pressure, bulk temperature rise is not 
particularly influenced by changes in chamber pressure. 
Hence, the bulk temperature rise is principally dependent on 
selection of motor configuration and propellant combination. 

With pressure and bulk temperature rise essentially deter- 
mined by the selection of motor configuration, propellant 
combination, and combustion pressure, only velocity, of the 
three parameters affecting q,,, is left for manipulation to gain 
cooling advantage. Fortunately, ¢,, can be readily increased 
by increasing velocity. However, since total coolant flow 
rate is fixed, velocity can only be increased by decreasing the 
net flow area of the cooling passage. With axial flow cooling 
in thin annular passages, it is found that velocity is inversely 
proportional to the coolant passage height ¢ and that the 
coolant-passage pressure drop increases with the third power 


Table 2 Effect of coolant velocity on value of g,,; and on 
turbopump power requirements 


First Second 


Item design design 
Coolant velocity at throat (ft/sec) 50 100 
Coolant bulk temperature at throat 
(°F)? 100 100 
Chamber pressure (psia) 300 300 
Injector pressure drop (psia) 75 75 
Coolant passage pressure drop (psi)® — 75 600 
Turbopump outlet pressure (psi) 450 975 
Pressure at throat in coolant passage 
(psia) 423 675 
Value of gy; at throat (Btu/in? sec )° 4.2 ‘ef 


Total propellant consumption re- 

quired for fuel turbopump (%)@ 2.0 4.25 
Increase in value of qui (%) a5 83 
Increase in total propellant con- 

sumption required for fuel turbo- 

pump (%) 2.25 


* Calculation for NH; cooled motor. 

Coolant passage pressure drop proportional to 

© Interpolated from Figures 2 and 4. 

“For constant mass flow rate turbopump power require- 
ment assumed proportional to head; assuming constant effici- 
ency, propellant consumption is proportional to power. 
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Fig.5 Ratio of pressure drop to adiabatic pressure drop vs. heat 
flux for RFNA and for Corporal fuel 


of velocity for turbulent flow. A sample calculation for a 
300-psi chamber pressure, turbo-pump-fed, ammonia-cool:d 
rocket motor (summarized in Table 2) shows the cooling ad- 
vantage and performance penalty resulting from doubling the 
coolant velocities. While the calculation is, of necessity, not 
general, the results are considered to be typical. They show 
that the value of q,, can nearly be doubled for as little as about 
2 per cent increase in propellant consumption (for turbo- 
pump) and perhaps a moderate increase in system dead weight. 
Hence, such changes should receive full consideration during 
preliminary design of a regeneratively cooled rocket propul- 
sion system. 


Dependence of Pressure Drop on Heat Transfer 


Considerations of coolant-passage pressure drop such us 
discussed in the section on Dependence of q,,, on System Oper- 
ating Conditions require knowledge of the manner in which 
the pressure drop varies with type and magnitude of heat 
transfer. A convenient method of showing this variation is 
through the use of the ratio of pressure drop per unit length 
with heat transfer to that without heat transfer for equal bulk 
temperature, and mass flow rate. This ratio is plotted vs. 
heat flux in Fig. 5 from experimental data for two propellants, 
RFNA and Corporal fuel (a mixture of aniline, furfuryl alco- 
hol, and hydrazine). The results for RFNA from (5) showed 
that Ap/Ap,.-o decreased from 15 per cent over the heat-flux 
range from zero to the inception of nucleate boiling (point B of 
Fig. 5). This decrease is somewhat less than would be pre- 
dicted by the Sieder-Tate friction factor correction (jw/u,)°'* 
which is usually recommended. Data presented in (16) for 
WENA were correlated with this viscosity ratio to the 0.12 
power. 

Pressure-drop data measured in the nonboiling region for 
n-butyl alcohol (32) and the Corporal fuel (26-29) decreased 
with increasing heat flux to a much greater extent than would 
be predicted by the factor (u,,/u,)°:'4, the resultant exponents 
being 0.3 and 0.2, respectively. Curve A’B’C’ (Fig. 5) drawn 
from data for the latter fuel at 100 F bulk temperature shows 
experimental values of the pressure-drop ratio as low as about 
0.4 at the inception of nucleate boiling. 

With heat fluxes increasing into the nucleate boiling range, 
the values of the ratio level off, then increase. Maximum 
values of the pressure-drop ratio at the upper limit of nucleate 
boiling both above and below unity have been observed, de- 
pending upon particular conditions of bulk temperature and 
pressure. The amount of pressure-drop data measured during 
the programs referenced was insufficient to achieve a correla- 
tion. Since the total pressure drop is an integration of friction 
on the cool wall under adiabatic conditions and on the hot wall 
both with and without boiling the general conclusion drawn 
from the data was that a calculation based on the adiabatic 
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conditions would be high by less than about 10 per cent for the 
acid cooling and 25 per cent for the Corporal fuel cooling. 

Unfortunately, very little data on pressure drop have been 
reported from which to establish correlations or draw general 
conclusions. This is partially due to the fact that pressure- 
drop measurements must be made over short sections since at 
nucleate boiling heat fluxes the bulk temperature rise per unit 
length is quite large. Such pressure-drop measurements are 
difficult to make with acceptable accuracy. Additional 
measurements are definitely needed to settle the questions 
concerning coolant-passage pressure drops. 


Operational Comparison of Cooling Capabilities 
of Several Propellants 


Having discussed various aspects of the utilization of pro- 
pellants as regenerative coolants, it is still not possible to cal- 
culate a parameter which will characterize the potentialities 
of a particular propellant as a coolant relative to other pro- 
pellants. This situation is a result of (1) the apparently com- 
plex relation between q,, and liquid properties and (2) the 
variables which are closely related to the selection of the 
combustion system, motor configuration, and combustion 
pressure. Although the final evaluation of the merits of a 
particular coolant must await a detailed design analysis, some 
measure of comparison between propellants can be achieved 
by considering several propellants as potential coolants for a 
hypothetical motor of given configuration and computing the 
margin between the local values of q,; of the coolant and the 
heat flux to the wall in a critical region such as the throat. 

Such a comparison was made using the propellants listed in 
the first column of Table 3 as the coolants. The propellants 
listed in the second column were considered to react with the 
coolant-propellant at peak-performance mixture ratios with 
resulting values of c* which were 95 per cent of theoretical 
values, based on equilibrium composition expansion. The 
propellant consumption was such as to produce 50,000 lb of 
thrust at chamber pressure near 300 psia with each combina- 
tion burning in the same motor. 

The pertinent configurational parameters of the motor 
selected were L* 40 in., throat area 117 in.?, 2:1 contraction 
area ratio, 7:1 expansion area ratio, 30-deg contraction half 
angle, 15-deg expansion half angle, and throat radius of curva- 


ture equal to throat diameter. The motor was assumed to 
have an annular coolant passage of height 4 which was con- 
stant along the whole length of the motor. The value of ¢ 
was allowed to vary from coolant to coolant so as to result in 
a 75-psi coolant-passage pressure drop for the appropriate 
flow rate of each particular coolant. Since the manner in 
which the friction factor for each coolant varies with heat flux 
is unknown, the adiabatic value was used for this comparison, 
using values of viscosity and density which were the average 
between inlet and outlet bulk temperatures. 

An inlet bulk temperature of 100 F was assumed for each of 
the coolants except liquid ammonia, which it was assumed 
could be reduced to 32 F by venting the tank prior to flight. 
The bulk-temperature rise of the coolant from the nozzle exit 
end to the throat and over-all bulk temperature rise were 
‘calculated from the coolant flow rate, specific heat and heat 
flux to the wall, computed with the following assumptions. 
Chamber gas temperatures 90 per cent of theoretical equi- 
librium combustion chamber temperature were assumed, the 
fraction being the result of assuming c* to be 95 per cent of its 
theoretical value. Wall temperature was 50 F (superheat) 
above coolant saturation temperature at 410 psi, the average 
coolant-passage pressure. This wall temperature results from 
assuming the combustion side wall to have negligible thermal 
resistance, and assuming the mechanism of heat transfer along 
the whole passage to be nucleate boiling. The local heat flux 
was computed by multiplying the difference between chamber 
gas temperature and wall temperature by heat-transfer co- 
efficients computed by the methods of (38). For reasons dis- 
cussed in (38, 39), the values of the computed heat flux are 
expected to be generally too high in the chamber, and 
throughout the nozzle as well, for carbonaceous systems which 
deposit carbon on the gas side walls during combustion. On 
the other hand, computed heat flux values are expected to be 
too low for systems burning either with appreciable dissocia- 
tion or with unstable combustion or both. These complicating 
factors cannot be accounted for analytically at present, but it 
is believed that values predicted on the basis of pure convec- 
tion are still useful for the comparison of this report; however, 
in drawing specific comparative conclusions between systems 
calculated, qualitative allowances should be made for these 
uncertainties. 

With these assumptions, it was possible to compute the 


Table 3 Comparison between calculated convective throat heat flux and throat values of g.: for ten propellant systems 


1 z 3 4 5 6 7 8 qui* porn 
value 
W V*,. Ty, Ref- at375 Trout, 
Coolant Reactant r ft/see °R see Ib/sec in. ft/see ° F sec erence psi °F 
a RFNA UDMH 2.5 5083 4796 5.85 158.0 0.109 60.8 179 8.8 6-14 355¢ 9-242 
b Corporalfuel SFNA 3.0 4763 4938 5.07 59.3. 0.0716 49.3) 243) 8.20 26-29 580F 
ec DETA SFNA 3.0 4930 4808 5.02 S73 0.0724 58.6 214 7.20* 26-29 740 306 
d 50 An-50FA SFNA 3.0 4740 4955 4.90 59.6 0.0719 50.0 236 5.40* 6 680t 346 
e JP-3 N20, 3.0 5069 5206 4.91 55.8 0.0712 63.4 238 5.25* 33 670¢ 349 
f NH; RFNA 2.2 4997 4224 5.58 70.7 0.0614 72.5 92 530 2325 140 140 
g Isopropanol SFNA 4757 55.3 0.0744 630 207 4.50* 6 550 293 
h NH; Oz 1.41 5499 4882 6.52 85.3 0.0688 78.0 92 5.80 23-25 140 135 | 
i JP-3 Os 2.25 5593 5577 6.00 62:5 67,1 276 5.16" 33 6707 
j CIF; NH, 2.50 5546 5881 5.83 145.0 0.0912 57.8 213 4.25 20 258 304 


* Extrapolated from lower velocity and lower bulk temperature values. 
+ Estimated from nucleate boiling wall temperature measurements. 


Notes: 1. 


Approximately peak performance (Isp @ 300 psi) mixture ratio, oxidizer/fuel. 2. 95 per cent of theoretical value 


based on equilibrium expansion. 3. 90 per cent theoretical value. 4. Computed from g* = hg*(7', — Tw): hg+ computed from Ref. 25 Cp 
based on constant composition value from thermochemical performance calculation, 7'~ taken to be 50 F higher than saturation tempera- 
ture of the coolant at throat pressure, 410 psia. 5. Value based on total utilization of propellant for regenerative cooling. 6. Computed 
coolant passage height for 75 psi coolant passage pressure drop, using isothermal friction factor values, and average values of density and 


viscosity between inlet and outlet temperatures. 
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7. Coolant velocity of throat. 8. Coolant bulk temperature at throat, based on 100 
F inlet bulk temperature for all coolants except NH; where 32 F was used. 
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CORPORAL FUEL (SFNA) 
LS 


3° 


HEAT FLUX (Btusin® sec) 
Fig. 6 Comparison between calculated throat convective heat 
flux and throat values of qu of ten propellant systems (see Table 
III for performance and mixture ratio) 


values of pressure, velocity and bulk temperature of each 
coolant at the throat. Values of q,,, at these conditions, (q,,:*), 
were determined from experimental data and compared with 
computed values of throat heat flux, g*. These results are pre- 
sented in Table 3 and Fig. 6. 

It is most significant to note that the computed bulk tem- 
perature rise with the CIF; system, is such that the saturation 
temperature is exceeded by the bulk temperature within the 
coolant passage, and saturation temperature is nearly reached 
with the NH; systems. Successful cooling even in the cham- 
ber under such conditions is extremely unlikely since values of 
dur drop well below 1 Btu/in.? see when bulk temperature 
reaches saturation temperature. This isa good point at which 
to start a proposed regenerative cooling design. 

Of the systems calculated, it appears that the RFNA, the 
Corporal fuel and DETA coolants offer the best margins for 
cooling, the values of qg,,* being comfortably above throat heat 
flux values. It must be remembered that the expected throat 
heat flux for the JP3-O, system would probably be as much 
as 50 per cent lower than predicted because of expected effects 
of carbon deposition. Hence, in practice, this system would 
also have a satisfactory but very unpredictable margin for 
cooling. Systems d, e, and f were predicted to be marginal 
within the uncertainties in predictions of heat flux to the wall. 
For most fuel-cooled systems, shifts toward lower mixture 
ratios, and for most oxidizer cooled systems, shifts toward 
higher mixture ratios, could be made to gain cooling advantage 
with only negligible penalties in performance. 

As indicated by the calculation in Table 2, even a coolant 
with apparently poor cooling capabilities such as isopropanol 
could be employed successfully if significantly greater coolant 
passage pressure drop is allowed in order to increase coolant 
velocity. 


Conclusions 


From such calculations and comparisons as presented in 
this report, some general conclusions can be drawn relative to 
liquid physical and chemical properties and to motor-operat- 
ing conditions which are favorable for regenerative cooling. 

In order to achieve successful regenerative cooling in a 
rocket motor with propellants below their critical temperature 
and pressure, sufficient coolant flow must be provided to 
achieve local values of the heat flux at the upper limit of 
nucleate boiling g,, which, everywhere throughout the motor, 
exceed corresponding local heat fluxes from the combustion 
gases to the wall. 

The total heat rate that can be absorbed by a coolant at 
bulk temperatures below saturation temperature is equal to 
the product of flow rate, specific heat and difference between 
saturation temperature and coolant inlet temperature. This 
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total heat rate must be safely above the total expected heat- 
rejection rate from the combustion gases to the motor wall 
since values of q,; decreases sharply as bulk temperatures ap- 
proach saturation temperature, becoming too low for success- 
ful cooling even back near the injector. For this reason a 
high specific heat, and/or a saturation temperature at least 
about 300 F is desirable. However, higher saturation tem- 
peratures than about 700 F allow wall temperatures to climb 
too high before being limited by nucleate boiling. Even this 
temperature would prohibit the use of aluminum as a wall 
material and would significantly reduce the allowable stresses 
of steel materials. 

A liquid critical pressure such that coolant passage pres- 
sures are from about 0.3 to 0.7 of critical pressure is desirable. 
Liquids operating in this pressure range usually exhibit a 
maxima in values of q,, for given velocity and bulk temper:- 
ture. Fluids which must be used at pressures above their 
critical pressure are usually poorer coolants than those 
operated below critical pressure, especially if the bulk tem- 
perature goes from subcritical to supercritical somewhere in 
the coolant passage. 

A low viscosity is a favorable liquid physical property, since 
it permits higher coolant velocities for a given pressure drop 
at a given mass flow rate. Furthermore, increasing viscosity 
is believed to be the explanation for the two cases mentioned 
where ¢,,, Was observed to decrease with decreasing bulk tem- 
perature, contrary to the normal trend. 

Liquids which decompose either at limited rates, or ex- 
plosively, which deposit salts or carbon on the hot wall, or 
which corrode the wall are generally inferior to liquids which 
are essentially chemically inert under coolant-passage 
conditions. 

The cooling of a motor with the oxidizer is best accom- 
plished by operating the motor at high-oxidizer-to-fuel mix- 
ture ratios. Conversely, cooling with the fuel is best ac- 
complished at low mixture ratios. For many systems, sig- 
nificant changes in mixture ratio from the peak performance 
value can be made to gain cooling advantage with but little 
sacrifice in performance. 

The only design parameter that is available to the designer 
to improve cooling conditions once the propellants, mixture 
ratio, combustion pressure, and motor configuration are 
selected is the coolant velocity which may be advantageously 
increased by decreasing coolant passage flow area with an at- 
tendant increase in coolant passage pressure drop. For some 
systems cooling with both propellants may be required. 


References 


1 Bankoff, S. G., Colahan, W. J. and Bartz, D. R., Summary 
of Conference on Bubble Dynamics and Boiling Heat Transfer 
Held at the Jet Propulsion Laboratory, June 14-15, 1956, 
Memorandum 20-137. Pasadena: Jet Propulsion Laboratory, 
Dec. 10, 1956. 

2 Sieder, E. N. and Tate, G. E., ‘‘Heat Transfer and Pres- 
sure Drop of Liquid in Tubes,’ Industrial and Engineering 
Chemistry, vol. 28: 1936, pp. 1429-1436. 

3 Ellion, M. E., ‘A Study of the Mechanism of Boiling Heat 
Transfer,’ Memorandum 20-88. Pasadena: Jet Propulsion 
Laboratory, March 1, 1954. 

4 Forster, H. K. and Zuber, N., ‘Dynamics of Vapor Bub- 
bles and Boiling Heat Transfer,’’ Journal American Institute of 
Chemical Engineers, vol. 1, 1954, p. 531. 

5 Dean, L. E., ‘Heat Transfer Characteristics of RFNA,” 
Report 117-982-004. Buffalo: Bell Aircraft Corporation, July 
1956. 

6 Combined Bimonthly Summary No. 36 (April 20 to June 
20, 1953). Pasadena: Jet Propulsion Laboratory, July 20, 195: 
(Confidential). 

7 Combined Bimonthly Summary No. 37 (June 20 to Aug. 
20, 1953). Pasadena: Jet Propulsion Laboratory, Sept. 20, 195: 
(Confidential). 


JET PROPULSION 


OPERATIONAL COMPARISON FOR: 
50,000 THRUST, 2:1 CONTRACTION, 7:1 EXPANSION, 40° L", 
COOLANT (REACTANT) ABOUT 300 psic, 95% OF C BASED ON EQUILIBRIUM, 2 
FOR MAX. Ig, 75 psi COOLANT PASSAGE PRESSURE DROP 
q 
F 
i 
2 
A 
ge 
| cr 
to 
Je 
1¢ 
th 
de 
of 
st 
1g 


neat- 
wall 
ap- 
CeSs- 
on a 
least 
tem- 
limb 

this 
wall 


Tes- 
ible. 
it a 
heir 
hose 
emi- 
ein 


ince 
sity 
ned 
em- 


ex- 
age 


ym- 
1iX- 
ac- 
sig- 
nce 
ttle 


ner 
ure 
are 
sly 
at- 
me 


8 Combined Bimonthly Summary No. 38 (Aug. 20 to Oct. 
20, 1953). Pasadena: Jet Propulsion Laboratory, Nov. 20, 1953 
(Confidential). 

9 Combined Bimonthly Summary No. 39 (Oct. 20 to Dec. 
20, 1953). Pasadena: Jet Propulsion Laboratory, Jan. 20, 1954 
(Confidential). 

10 Combined Bimonthly Summary No. 40 (Dec. 20, 1953 to 
Feb. 20, 1954). Pasadena: Jet Propulsion Laboratory, March 20, 
1954 (Confidential). 

11 Combined Bimonthly Summary No. 41 (Feb. 20 to April 
20, 1954). Pasadena: Jet Propulsion Laboratory, May 20, 1954 
(Confidential). 

12. Combined Bimonthly Summary No. 42 (April 20 to June 
20, 1954). Pasadena: Jet Propulsion Laboratory, July 20, 1954 
(Confidential). 

13 Combined Bimonthly Summary No. 43 (June 20 to Aug. 
20, 1954). Pasadena: Jet Propulsion Laboratory, September 20, 
1954 (Confidential). 

14 Combined Bimonthly Summary No. 50 (Oct. 1 to Dee. 1, 
1955). Pasadena: Jet Propulsion Laboratory, December 15, 
1955 (Confidential). 

15 Wolf, H., Gray, F. L. and Reese, B. A., “Heat Transfer 
and Friction Characteristics of Red and White Fuming Nitric 
Acid,”’ Jet Propusion, vol. 26, Nov. 1956, p. 979. 

16 Reese, B. A. and Graham, R. W., “Experimenta! Investi- 
gation of Heat Transfer and Fluid Friction Characteristics of 
White Fuming Nitric Acid,’’ Technical Note 3181. Washington: 
NACA, May 1954. 

17 Reese, B. A., and Graham, R. W., ‘Heat Transfer and 
Friction Pressure Drop Characteristics of White Fuming Nitric 
Acid,”’ Jev Proputsion, vol. 24, no. 4, July-Aug., 1954, p. 228. 

18 Ashley, E., “Heat Transfer Measurements for White 
Fuming Nitric Acid, Report No. 56-982-016. Buffalo: Bell Air- 
craft Corporation, Feb. 1953. 

19 Hatcher, J. B. and Bartz, D. R., “High Flux Heat Transfer 
to JP-3 and RFNA,” External Publication No. 119. Pasadena: 
Jet Propulsion Laboratory, November 1951. 

20 Combined Bimonthly Summary No. 56 (Oct. 7 to Dee. 1, 
1956). Pasadena: Jet Propulsion Laboratory, Dec. 15, 1956 
(Confidential). 

21 Powell, W. B., “Heat Transfer to Fluids in the Region of 
the Critical Temperature,’ Progress Report No. 20-285. Pasa- 
dena: Jet Propulsion Laboratory, April 1956. 

22 Banchero, Baner and Ball, ‘““Heat Transfer Characteristics 
of Boiling Oxygen, Fluorine, and Hydrazine,’”’ Report 51-11. 
Ann Arbor: University of Michigan Engineering Research In- 
stitute. 

23 Combined Bimonthly Summary No. 52 (Feb. 1 to April 1, 
1956). Pasadena: Jet Propulsion Laboratory, April 15, 1956 
(Confidential). 

24 Combined Bimonthly Summary No. 53 (April 1 to June 1, 
1956). Pasadena: Jet Propulsion Laboratory, June 15, 1956 
(Confidential). 


JANUARY 1958 


25 Combined Bimonthly Summary No. 54 (June to Aug. 1, 
1956). Pasadena: Jet Propulsion Laboratory, Aug. 15, 1956 
(Confidential). 

26 Combined Bimonthly Summary No. 44 (Aug. 20 to Oct. 20, 
1954). Pasadena: Jet Propulsion Laboratory, Nov. 20, 1954 
(Confidential). 

27 Combined Bimonthly Summary No. 45 (Oct. 20 to Dee. 
20, 1954). Pasadena: Jet Propulsion Laboratory, Jan. 20, 1955 
(Confidential). 

28 Combined Bimonthly Summary No. 46 (Dec. 20, 1954 to 
March 31, 1955). Pasadena: Jet Propulsion Laboratory, April 1, 
1955 (Confidential). 

29 Combined Bimonthly Summary No. 47 (April 1 to June 1, 
1955). Pasadena: Jet Propulsion Laboratory, June 15, 1955 
(Confidential). 

30 Combined Bimonthly Summary No. 51 (Dee. 1, 1955 to 
Feb. 1, 1956). Pasadena: Jet Propulsion Laboratory, Feb. 15, 
1956 (Confidential ). 

31 Kreith, F. and Summerfield, M., “‘Heat Transfer From an 
Electrically Heated Tube to Aniline at High Heat Flux,’’ Progress 
Report No. 4-88. Pasadena: Jet Propulsion Laboratory, Feb. 
1949. 

32 Kreith, F. and Summerfield, M., “Investigation of Heat 
Transfer at High Heat Flux: Experimental Study of Heat Trans- 
fer and Friction Drop with n-Butyl Alcohol in an Electrically 
Heated Tube.’ Progress Report No. 4-95. Pasadena: Jet Pro- 
pulsion Laboratory, May 1949. 

33 Hatcher, J. B., “High Flux Heat Transfer and Coke Dep- 
osition of JP-3 Fuel Mixture,’ Progress Report No. 20-157. 
Pasadena: Jet Propulsion Laboratory. 

34 Beighley, C. M. and Dean, L. FE. “Study of Heat Transfer 
to JP-4 Jet Fuel,” Jer Propunsion, vol. 24, No. 3, May-June 
1954, p. 180. 

35 Combined Bimonthly Summary No. 57 (Dee. 1, 1956 to 
Feb. 1, 1957). Pasadena: Jet Propulsion Laboratory, Feb. 15, 
1957 (Confidential). 

36 Gunther, F. C., “Photographic Study of Surface-Boiling 
Heat Transfer to Water with Forced Convection,’ Transactions 
of the American Society of Mechanical Engineers, vol. 73, 1951, 
pp. 115-123. 

37 Rohsenow, W. M. and Clark, J. A., Technical Report No. 
4, DIC Project No. 6627. Cambridge: Massachusetts Institute 
of Technology, 1952. 

38 Bartz, D. R., “A Simple Equation for Rapid Estimation of 
Rocket Nozzle Convective Heat Transfer Coefficients,’ JET 
PROPULSION, vol. 27, Jan. 1957. 

39 Bartz, D. R., “An Approximate Solution of Compressible 
Turbulent Boundary Layer Development and Convective Heat 
Transfer in Convergent-Divergent Nozzles,’ Progress Report No. 
20-234. Pasadena: Jet Propulsion Laboratory, July 1954 (also 
Transactions of the American Society of Mechanical Engineers, vol. 
77, no. 8, Nov. 1955, p. 1235). 


3 
| 
fer 
D6, a 
ry, 
es 
ri 
on 
ib- 
| 
hy 
. 
ig. 
5 
53 


Technical Notes 


Criteria for Orbital Entry 


LOUIS G. VARGO! 


The Ramo-Wooldridge Corporation 


HE PROBLEM of entry into a circular orbit has been in- 

vestigated using the concept of characteristic velocity 
(1).2. We may regard use of the characteristic velocity as a 
definition of an optimization criterion. With this criterion, an 
important result shown in Fig. 3 of (1) is that the character- 
istic velocity reaches a maximum value at a finite distance 
from the departure point when the Hohmann ascent method 
(2) isemployed. This note gives a corresponding result for a 
different criterion, viz., the radial distance at which the pay- 
load reaches a minimum for a fixed initial weight. 

Consider the final stage(s) of total initial weight wo boosted 
to an initial tangential velocity 1 at a radial distance 1. 
Equations [1 and 2] give the initial and final velocity incre- 
ments required to achieve circular velocity at a radial distance 


1 20. 
Av, = = ——].......... 2 
[2] 


Here o = r,/ro and y is the constant of the gravitational field 
under consideration. The characteristic velocity, Avy + 
Av,, may be differentiated to obtain the cubic equation 


— — 90 — 1 =0.............B] 


which has the positive root ¢ ~ 15.6. This is the aforemen- 
tioned result of (1). Let w denote the payload weight. The 
relationship between the over-all payload ratio w/w) and o 
takes two forms depending on whether or not weight is jetti- 
soned after Avy is added. Equations [4 and 5] give this rela- 
tionship for the two cases. 


w 
jettisoning 
Wo 


Co Cs 
— — &,) 


In Eq. [4], & is the empty weight factor relative to wy — w; 
in [5], &o and k, denote the empty weight factors relative to 
the gross weights of the two steps which total w — w. c) 
and ¢, are the exhaust velocities of the propulsion systems pro- 
ducing the two velocity increments. 
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From an inspection of [4] we see that the payload ratio is 
monotonically related to the expression Ay/t + Av,/c.. 
The minima of w/wo occur at points (omin, Cs/¢o) which 
satisfy 


— E -3| o? — 
Cs 
+2) 


Fig. 1 shows the essentially linear dependence of omin on 
c./C for a range of exhaust velocity ratios of interest. Thus, 
for the case of no jettisoning, the two criteria give identical 
results only if c, = ¢o. 

If weight is jettisoned, Eq. [5] can be differentiated to ob- 
tain 

exp [— afi(a) — abfo(o) + d] — k, exp [— afi(o) + d] 


exp [— — abfo) + d] — ko exp[— abfe(o)] 
bf(o).... [7] 


Kprtor’s Note: This section of Jer PROPULSION is open to short manuscripts describing new developments or offering comments on 
apers previously published. Such manuscripts are published without editorial review, usually within two months of the date of receipt. 
equirements as to style are the same as for regular contributions (see masthead page of this issue). 
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As an example of the use of Eq. [7] consider the numerical 
casea = d = 3,b =1andk) =k, =k. Fig. 1 showsa plot 
of omin as a function of k. 

The above results point out the necessity for separate 
analyses of orbital entry problems when both vehicle and tra- 


jectory criteria are considered. 
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Minimality for Problems in Vertical and 
Horizontal Rocket Flight 


STANLEY ROSS! 


Lockheed Missile Systems Division, Palo Alto, Calif. 


Nomenclature 
t = time 
rc = instantaneous range 
y = instantaneous height 
= instantaneous velocity 
» = instantaneous mass 
4 = gravitational acceleration, assumed constant 
c = relative exhaust velocity, assumed constant 
m, = payload 
£ = gx/c? nondimensional range 
» = gy/c? nondimensional height 
vy = v/ec nondimensional velocity 
u = m/m, nondimensional mass 
= me” 
t = poe where W> is proportional to Cp, the drag coefficient 
) = modified drag function, assumed dependent on »v and y 
only 
\ = Lagrange multiplier function 


Introduction 


OCKET thrust programs for stationary values of fuel 
consumption have been investigated by several authors. 
In particular, vertical motion of the sounding rocket has been 
studied (1-3),? as has the case of purely horizontal flight (4). 
These treatments have employed some combination of maxi- 
mum-thrust ares (boosts), varying-thrust arcs (sustain) and 
zero-thrust arcs (coasting) to establish an over-all optimum. 
By a simple artifact, a condition may be found for each case 
by whose satisfaction that entire trajectory furnishes not 
only a stationary path, but actually a true local minimum. 


The Sounding Rocket 


We define a new dependent variable ¢ = we’, and seek 
those functions g(n) and v(n) which minimize ¢(0) = g, 
under the constraint? 


Case I: 


+" + D0) = 0... 
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such that »(0) = 0, v(m) = 0, ¢(m) = 1. (Differentiation 
is with respect to 7, here the independent variable.) This 
may be associated with Mayer’s problem (5), for which the 
following extremals arise 


(Fy =F, 


(F = AW). From Equation [2], \ = Agexp(gt/c). Since 
Equations [2 and 3] are homogeneous in X and X’, the multi- 
plier may be assigned an arbitrary nonzero‘ value at a par- 
ticular point of the interval (0, m). We choose A(0) = A» = 1, 
from which it follows that \ # 0. Therefore, on the path in- 
dicated by [2 and 3] 


g=vD,.. [+] 


When D = Qve’~*”, i.c., quadratic drag law, condition [4] re- 
duces to 


p= + in (O, m). [5] 


as has been noted previously (3). The second variation of 
~y can be shown® to be of the form 


1 1 
+ 2F + ... . [6] 
0 


Now, 
dF = + F6¢ + Fyov' + = 0 
and by [1, 3, 2] respectively 


Fy =0; F,=0: =F 


Thus (F = 0, in (0, 7). At» = 0, é¢ = 0 for a station- 
ary value of g. And since Py = X + 0, we conclude that 
6¢ = 0 identically in (0, 4:). Equation [6] becomes 


= f, F,,6v°dn 


A necessary and sufficient condition for a true local minimum 
is then that F,, > 0, along the trajacetory. In particular, 
whenever D is of the form Qv’e’—* then 


2 
Fy, =X + Qv?-2[(v + p)? ple” One 


which is positive for all values of p, except perhaps for 0 < p 
<1. Specifically, it is positive for both linear and quadratic 
drag laws. 

Equation [4] determines values of the two end-point veloci- 
ties, which will be different from zero. Terminal boosts of in- 
finite thrust are therefore introduced in order that the re- 
quired velocity discontinuities may be achieved at these points. 
And since ¢ is continuous across such instantaneous boosts, 
¢o is still minimal over the entire path, including these boosts. 
However, if the thrust is to be bounded by a zero minimum 
value (i.e., coasting), it may be necessary to join this ‘““bound- 
ary-curve” program to the path defined by Eq. [3 and 4] in 
such a way that the desired end conditions are attained.‘ 
The over-all trajectory in such cases is then a constrained 
minimum,’ subject to the above thrust limitation. 


‘If \ = 0 at any point where [2 and 3] hold then it must be 
zero identically. Then it can be shown that 6F = 0, without 5¢o 
necessarily being zero, which is not a legitimate solution. 

5 Cf. reference (5), page 23, Eq. 7. 

6 The Weierstrass-Erdmann First Corner Condition must be 
satisfied at all junctions between two dissimilar extremals. Here 
the condition is that \ be continuous across each junction. 

7 Cf. reference (6), page 42. 
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Case Il: Hibbs’ Problem of Optimum Horizontal Flight® 
In this case gp is to be minimized subject to 
W =o! + Kwe’ + = 0......... [7] 


where K, and K; are constants, and the last term on the left 
represents the effect of lift-induced drag. (Differentiation is 
with respect to &, here the independent variable.) Analysis 
following that of the previous case shows that 

(a) By the Euler-Lagrange equations 
1+» 
is the burning program to be followed for a stationary go 


(b) 6g =0, in (£o, &) 
(c) The second variation reduces to the form 


f 


where /’,, > 0 is again a necessary and sufficient condition 
for a true local minimum. 

(d) Here F,, > 0 always, and the burning program indi- 
cated in (a) is to be employed, unless physical restrictions 
dictate the insertion of coasting arcs. Insertion of terminal 
boosts is also indicated, and the solution again provides either 
a free, or thrust-restricted, local minimum. 


p= Ky 


5 Cf. reference (4). 
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On Turbulent Heat Transfer Through 
a Highly Cooled Partially Dissociated 
Boundary Layer' 


RONALD F. PROBSTEIN,? MAC C. ADAMS,? PETER H. 
ROSE* 


Aveo Research Laboratory, Everett, Mass. 


T IS the purpose of this note to present some analytical 
results regarding the properties of a compressible dissociat- 
ing turbulent boundary layer in thermodynamic equilibrium. 
For brevity, detailed derivations are omitted and symbol 
definitions are given only when it is not clear what is meant. 
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Fig. 1 Fluid property distributions through a highly cooled 


partially dissociated turbulent boundary layer. 


Of interest are the mean continuity, momentum and energ\ 
equations for the plane or axisymmetric turbulent boundary 
layer in thermodynamic equilibrium. These can be derived 
from the general equations for a dissociating boundary layer 
(1)* by writing the fluid properties as the sum of a time averaged 
quantity plus a fluctuating quantity (e.g.,u = @+ u’). I 
in addition we assume the Prandtl] number Pr = é,u/k = |, 
and the Lewis number, Le = ¢pD/k = 1 (where D is the 
difiusion coefficient) these equations can be expressed in the 
usual form (2) 


(piir’)x + (pir’)y + (p’v'r’)y = 0....... 


Pitts + (pd + p'v')ily = + (Wily + Mrthy)y..... [2] 


pall, + (pi + p’v’)H, = A,+— 
Cp Cp 


Here, r is the cross-sectional radius of the body with 7 = 0 for 
plane flow and j = 1 for axi-symmetric flow, FT = hiacat + 
Naissoe + 7/2 = h + a@/2 is the total enthalpy, and wr = 
—pu'v'/i, is the eddy viscosity. In order to derive Equations 
[1-3] the additional definitions and assumptions not given in 
(2) are: (a) The eddy conductivity is defined as ky = 

,and (b) the turbulent Lewis number Ler = 

u/ Cp 

=1, with D> the turbulent diffusion coefficient. Thus analo- 


gous to the fact that in a dissociating gasthe heat transported by 
ordinary conduction and diffusion is - h, only for Le = 1, 
see (1), so too is the rate of energy transported by turbulent 
conduction and diffusion given by —j v’h’ = tt h, only when 


Ley = 1. As with the turbulent Prandtl] number, Ley is 
probably nearer to unity than the laminar Lewis number. 

If in addition to Le = Pr = 1, we have Prr = épur/kr = 1, 
with 
en with 
equations are similar. For either an insulated body or a con- 
stant surface enthalpy we immediately obtain the familiar 
Crocco integral between enthalpy and velocity 


= 0, both the mean momentum and energy 


H = he + (He he) [4] 


where the bars have now been dropped, and the subscripts 
and e represent conditions at the wall and edge of the bound- 


5 Numbers in parentheses indicate References at end of paper. 
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ary layer respectively. From this relation we note that for a 
highly cooled boundary layer the enthalpy and temperature 
distributions will not have a maximum within the layer as in 
the usual supersonic case with moderate temperature dif- 
ferences. This follows by differentiating Eq. [4] with respect 
to y and noting that since u, > 0, then h, > Oif (H. — hw)/u. 
> u,. For large wall cooling where h,/H.— 0, this implies 
u?/H. < 1 for hy > 0. This condition will in general 
be satisfied on a blunt body in hypersonic flight where 
u-’/H.< 1. Interpreted in terms of loca] external Mach num- 
ber u.2/He = 1 corresponds to M, = 2.24 for y = 1.4 and 
M, = 3.16 for y = 1.2. In general other fluid properties such 
as density p, temperature 7’, etc., will follow the same trend 
as the enthalpy and will have no maximum. Some typical 
distributions of static enthalpy h, temperature, density, and 
degree of dissociation a, obtained using Eq. [4] are shown in 
Fig. 1 for the case of a highly cooled turbulent boundary layer 
in thermodynamic equilibrium with an assumed velocity pro- 
file of the form u/u, = (y/6)'/7.. The nature of the distribu- 
tion is not critical to the velocity profile so long as it is of the 
turbulent type, i.e., a full profile. 

From Eq. [4] Reynold’s analogy between heat transfer and 
skin friction follows; that is, the Stanton number defined with 
stagnation enthalpy instead of recovery enthalpy is 


_ Nu q 7 Cy 
Re.Pr peue(H. — hw) Pee” 2 


St 


where Nu is the Nusselt number and Re, the Reynolds num- 
ber. For Prandtl numbers different from unity it is usual to 
apply a correction factor (3-4), which turns out to be ap- 
proximately the same as for laminar flow, namely St = o 
Pr-2/3, 

Analogous to this laminar Prandtl number correction we 
may expect that a laminar Lewis number correction for values 
different from unity should also be applied for a dissociated 
flow. Such a correction is to be expected since, as indicated in 
Fig. 1, atom concentrations persist to the extremities of the 
laminar sublayer. The form of this correction can be seen by 
considering the heat flux crossing any plane which is 


where h° is the dissociation energy per unit mass of atomic 
products, D, the atomic diffusion coefficient, and c, the atom 
mass fraction. From Eq. [6] we find that the heat transfer 
ratio is 


Vie = 1 {hideat + h°caly 


Now for the one-dimensional case of no convection in which 
we have linear distributions between a higher temperature 
surface and a cold wall Eq. [7] reduces to (dz # 1)/(due = 1) 
= 1 + (Le — 1)Ap./H., where hp, is the energy in dissociation 
at the high temperature surface, and h, = H, the total en- 
thalpy at this surface. On the other hand for a flat plate 
highly cooled laminar boundary layer this correction enters 
approximately as 1 + (Le?/? — 1)hp./H., see (5) and (1), with 
the change in the exponent resulting from the convection 
present in the boundary layer. Without going into the de- 
tailed reasoning it can be shown that it is logical to assume for 
a highly cooled turbulent boundary layer where the profiles 
in the laminar sublayer are nearly linear that 

hoe 2 


die #1 
Le® — 1) — ,-<B<1........ 8 


with recourse to experiment required for the determination of 
the exact value of B. 

With the previous considerations in mind it is now possible 
to formulate for a blunt body of revolution, for the practical 
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hypersonic case of high wall cooling, a semi-empirical heat 
transfer rate theory dissociated turbulent boundary layers 
in thermodynamic equilibrium. The simplification follows 
from the fact that for such a case it is not difficult to show that 
the boundary layer behaves very much as an incompressible 
turbulent boundary layer (see Fig. 1), and the pressure 
gradient affects the boundary layer growth, but not markedly 
the shapes of the profiles. In particular it can be shown that 
so far as the skin friction coefficient is concerned the boundary 
layer behaves as if it were an incompressible constant property 
layer with a pressure gradient, so that it is justified to use the 
Ve 
ud 
velocity profile represented by a power law. Justification of 
this is based on the reasoning that: 1. The major changes in 
fluid properties occur through the laminar sublayer, which 
transmits the turbulent shear to the wall; 2. the shear is 
approximately constant through the sublayer; and 3. the 
boundary layer growth is practically unaffected by the sub- 
layer, since only a small fraction of the total momentum is in- 
cluded there. 

The boundary layer growth can now be determined utilizing 
the integrated momentum equation 


24 2 due 
Pe dz 


Blasius incompressible relation C,/2 = 0.0225 ef and a 


dr Ue az r dz 


To integrate the above equation we require expressions for the 
momentum thickness 0, and the form factor Hr = 6*/6. 
Using Eq. [4] and a 1/7 power law these relations are found 
for a highly cooled blunt body to be; 6 = (7/72)6(1.2), with 
1.2 representing a mean value, and Hy < 2. Thus with 6 = 
6(6) known and Hy ~ 0 we find on integrating Eq. [9] from 
x = 0 where wet set 6 = 0 that 


0.289 fy 
[10] 
Jo 


Ue 
Utilizing Eq. [10], the Blasius skin friction relation, Eq. 
[5], and the Lewis and Prandtl] number corrections, we can 
now write the final heat transfer equation for a highly cooled 
turbulent boundary layer subject to a pressure gradient as 
Nu 

(Re;)°-8 
where we remind the reader that Nu and Rex, are based on 


properties at the edge of the boundary layer. Here F(z) isa 
form factor for the pressure gradient effect and is given by 


1 
F(z)" = in [12] 
le 0 le 


Clearly F(x) ~ 1 asx — © and 4 — 0, as for example far 


65/4 = 


= 0.029Pr'/3[F(z) + (Le® — 1) 


downstream on a hemisphere cylinder. Furthermore it can 
be shown that F(z) is always greater than 1 for a positive 
pressure gradient, and will always have a maximum in the 
region of maximum mass flow, and hence local Mach number 
of one 
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Heat Transfer to a General Three- 
Dimensional Stagnation Point 


ELI RESHOTKO! 


California Institute of Technology, Pasadena, Calif. 


The equations for the compressible laminar boundary 
layer at the stagnation point of a body having unequal 
principal curvatures are presented. A simple representa- 
tion for the heat transfer to such a body is obtained for the 
case of the very cold wall. Estimates based on this repre- 
sentation are made for the heat transfer at other tempera- 
ture levels. 


Nomenclature 


specific heat at constant pressure 

stream functions for similar solutions 

enthalpy 

ratio of principal velocity gradients (Eq. [11]) 

thermal conductivity 

Prandtl number 

pressure 

heat transfer rate to wall 

gas constant 

principal radii of curvature 

temperature 

transformed velocity in X and Z directions, respec- 
tively 

velocities in x, y, and z directions, respectively 

transformed coordinates along body surface 

coordinates along body surface 

transformed normal coordinate 

normal coordinate 

boundary layer similarity variable 

H dH, 

w 


density 

absolute viscosity 

kinematic viscosity 
components of stream function 


5 
= 


x 


Subscripts 


local condition outside boundary layer (external) 
stagnation point value 

wall value 

free-stream stagnation value 

free-stream static value 

A coordinate symbol used as a subscript generally denotes partial 
differentiation with respect to that coordinate. 


sp 
w 
0 


eo 


Boundary Layer Equations 
HE COORDINATE system chosen is shown in Fig. 1. 
In the vicinity of the stagnation point, these coordinates 
are Cartesian. The pertinent compressible boundary layer 
equations are 


Continuity 
(pu)e + (pv), + (pw), = 0.............. {1] 
Received November 12, 1957. 
1 Guggenheim Aeronautical Laboratory. On leave of absence 
from Lewis Flight Propulsion Laboratory, NACA. 
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Stagnation Point 
ee 
Fig. 1 Coordinate system for general three-dimensional 
stagnation point 
Momentum 
due 
pu Uz + pruy + pwu, = pet | + (muy)y.....[2a] 
/sp 
dw, 2 
+ pyvw, + pww, = + (pwy)y. . . [2c] 
dZ } sp 
Energy 
pub. + pvb, + = (4 [3] 
r 
State 


The boundary conditions for these equations are 


Aty =0 


u=v=w=60=0 


=%w= = 1 


The viscosity is assumed to be a linear function of the tem- 
perature according to the relation 


The velocities in Equations [1 and 2] can be replaced 
through the definition of the two component stream functions 


(Y, ¢) 
= 


Equations [6] automatically satisfy the continuity equation. 
Upon making the transformation of coordinates 


ion 
= d. Y= =e Z= d. 
I=) op o\ 


and introducing the following similarity transformation 


dU, aU, 


[8] 
Vo 


Equations [1-3] reduce to the ordinary differential equations 


Hy» 
= 1 = — 1)(1 — 6).... [9a] 


Hw 
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+ Pr(f + Kg)0’ =0.............. [9e] 
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. [2a] 
. [2b] 


. [2c] 


[3] 


[4] 
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with boundary conditions 


At =0 f=f'=g= 


The parameter K in Equations [9] is related to the principal 
velocity gradients 
(dW./dZ)sy _ 
(dU./dX ) sp 


(dw./dz)sp 
(du./dz)sp 


K= 


It is readily seen from Equations [9 and 11] that the case 
K = 0 corresponds to the two-dimensional stagnation point 
while = 1 represents the axially-symmetric stagnation point. 
If the X coordinate is identified with the larger principal 
curvature then all stagnation points have 0 < K < 1. 

For incompressible flow, H../Ho is taken as unity in Equa- 
tions [9a and 9b]. The resulting equations were originally 
obtained and numerically solved by Howarth (1).2 In 
principle, [9] with boundary conditions [10] can be solved for 
many values of K, H,,/Ho, and Prandtl number. It was sug- 
gested by Professor L. Lees, however, that perhaps a simple 
solution could be obtained for the very cold wall where H./Ho 
<1. This is in fact the case. 


Solution for Very Cold Wall 
For H,,/Hy = 0, Equations [9] become 


+ (f + Kg)g’”’ = K(g” — @)....... [12b] 
0” + Prif + Kg)0’ =0........... [12¢] 


It has already been shown in (2) that at this temperature 
level the right hand sides of [12a] and [12b] have a small effect 
on heat transfer and for simplicity can be neglected. Equa- 
tions [12] then become 


+ = ©...:.. 
+ + = 0......... 


Since the functions f and g satisfy the same boundary 
conditions (see Eq. [10]), then for Equations [13], f = g, and 
Equations (13) become 


6’ + Pr(l + K)f@’ = 


Now let F = fv + K and = 7V1+ K. Equations 
[14] then reduce to the Blasius equation and its accompanying 
energy equation 


+ FF) = of 
+ PrF@, = 


which have the well-known wall values (F),). = 0.470 and, 
for Pr = 1, (0,)~ = 0.470. Returning then to the functions 
of Equations [14] 


fol’ = 0.470071 + K 
= 0.470071 + K 


Caleulation of Heat Transfer 


The local rate of heat transfer per unit area to a stagnation 
point is expressed 


q= te) (Ao Hw) [17a] 
Cp Mw\dx sp 


*Numbers in parentheses indicate References at end of paper. 
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or 


q= (He — Hi)...... [17b] 


dx 


where the value of @.,’ is that for Prandtl number | and may be 
read from Fig. 2. The lowest curve in Fig. 2 is that of Eq. 
[16]. The other curves are estimates based on [16] of the 
variation of 6.’ for H./Ho = 0, 1, 2 where the end-points are 
obtained from the exact solutions of (3-5) for two-dimensional 
and axially-symmetric stagnation points. The variation of 
heat transfer at the stagnation point with Prandtl number is 
included in Equations [17] as Pr’. 

The velocity gradient (du./dx)sp is possibly evaluated from 
modified Newtonian flow using the expression 


due Pip — Po 18 


where R, is the smaller principal radius of curvature. The 
ratio of velocity gradients K is related to the principal radii of 
curvature by the expression 


IR 
Ve. 


The effects of dissociation for Lewis number of one can be 
roughly accounted for by using external fluid properties in 
place of the wall fluid properties in Eq. [17b] as was done in 
(2). A better approximation to the effects of dissociation 
would be obtained through multiplying the right side of 
[17b] by (Reference 6). 
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On the Rational Choice of Flame- 
Holder Shape 


ABBOTT A. PUTNAM! 


Battelle Memorial Institute, Columbus, Ohio 


N DESIGNING a flame-holder system, the first require- 
ment is that the flame-holder system be capable of main- 
taining a stable flame over a given range of mixtures and flow 
velocities. Next, this flame must propagate across the mix- 
ture from the holders so that all the fuel is consumed. Finally, 
combustion should take place with minimum pressure losses 
and in a system of minimum length. 

For simplicity, let us consider a combustor in which the most 
stringent requirement is related to flow velocity and to mix- 
ture ratio, and the length of the combustor is fixed. Because a 
flame can spread only so fast in a stream, the fixed length of 
combustor fixes the area served by each component of a grid 
flame holder. In an open system, the larger the flame holder 
the greater the velocity of final instability (1).2 However, 
when the flame holder is in a restricted region, the flow ve- 
locity past the flame holder increases as its size increases, for 
a fixed average flow velocity over the entire cross section. 
Thus, a point will be reached where increasing the size is not 
beneficial. It follows that for each shape of flame holder, and 
each mixture ratio, a curve may be drawn giving the maximum 
average flow velocity that may be obtained in a fixed size of 
duct or associated area. 

Fig. 1 shows plots of such curves for two different configura- 
tions of flame holder, designated by F and U, as a function of 
reciprocal area of the duct. Since all operating conditions of 
a given shape must fall below the corresponding curve on this 
plot, these curves may be termed envelope curves. The 
relative positions of the curves are related to the angle between 
the initial flame front and the flow direction (2). 

If the design condition is such that only one shape of holder 
will satisfy the flame-holding requirements, then there is no 
further problem. The shape has been fixed. However, sup- 
pose that several holders will satisfy the requirements. Then, 
the advantages of one shape relative to another must be 
weighed. The question of flame-holder drag must be con- 
sidered. 

Fig. 2 gives an indication of the variation of drag for a given 
shape of flame holder. The dashed curve, corresponding to 
either envelope curve of Fig. 1, may be considered to be com- 
posed of infinitesimal segments of curves of constant drag 
force. The left-hand parts of two such curves of constant 
drag force are shown as solid lines; the right-hand parts are 
not indicated because a greater flame-holding ability can be 
obtained with smaller drag using a smaller holder in the right- 
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Velocity 


Fig. 1 


Velocity 


Velocity 


Shape F 


Size of flame holder decreases 


Curves of maximum flame-holding ability for two shapes 


Reciprocal Area 


of flame holder 


A-23592 


Reciprocal Area 


Fig. 2 Curves of constant drag force relative to envelop curve 
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Use U 


Use F 


Fig. 3 


Reciprocal Area 
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Curves of constant drag holding for two shapes 
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hand region. It is seen from this figure that savings in drag 
force can be obtained by not overdesigning a holding system. 

The next question is whether it is better (a) to use the shape 
of holder, the envelope curve of which just meets the design 
requirements, or (b) to use a shape, the envelope of which is 
too high, but use a smaller size holder with a lower drag. 
This question is answered in (3) for the case 1/S = 0, that is, 


and duct size. 


for a holder in an infinitely wide stream. For minimum drag, 


one should use the shape associated with the lowest envelope 


curve that will meet the requirements. 


Fig. 3 shows the extrapolation of this information to finite 
values of duct area. A drag force is chosen, and then the 
curves of constant drag force are drawn. 
shape U should be used to obtain the maximum flame-hold 
ability for a given amount of drag, or minimum drag for a 
given amount of flame-holding ability. For 1/S— ©, orS— 


For 1/S — 0, 


0, shape F should be used The change-over point is a function 
of the design velocity. 

These considerations demonstrate that, to make a scientific 
decision as to the optimum shape of flame holder to choose, the 
flame-holding ability and hot drag of various shapes of flame 
holders must be known as a function of both flame-holder size 
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Book Reviews 


Physical Measurements in Gas Dynamics 
and Combustion, edited by R. W. La- 
denburg (Part 1) and B. Lewis, R. N. 
Pease and H. 8. Taylor (Part 2); Vol- 
ume IX, High Speed Aerodynamics 
and Combustion, Princeton University 
Press, Princeton, N. J., 1954. $12.50. 


Reviewed by F. H. Wricut 
Jet Propulsion Laboratory 


“Physical Measurements in Gas Dy- 
namics and Combustion’’ contains a wealth 
of information about experimental tech- 
niques and supplies sufficient background 
material to permit intelligent application 
of the techniques. A refreshing feature of 
the book is the choice of authors. The 
twenty-two authors of the several sections 
have developed the methods they write 
about, or have had long experience in their 
practical application. As a result, they 
write with vigor and authority and in a 
way that makes interesting reading. 

This book is much more than a hand- 
book for the experimentalist—it is also a 
book for the general reader. It succeeds 
in reaching both audiences. However, it 
is not completely successful in coping with 
the duality of presentation required. The 
task of satisfying both audiences is an al- 
most impossible one. For the experi- 
mentalist, complete information on each 
method is required, either within an arti- 
cle, or in easily available references. Eval- 
uation of techniques, analysis of accura- 
cies, delineation of fields of application, 
critical discussion of references and com- 
parison of a method with other ways of 
performing the same task are also needed. 
Further, the experimentalist wishes the 
author to share his wealth of experience 
with the reader, suggesting experimental 
tricks and pointing out pitfalls. 

On the other hand, for the general 
reader many of the practical details are 
unnecessary. The general reader needs a 
book that describes the experimental 
methods in a simple fashion that will lead 
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Ali Bulent Cambel, Northwestern University, Associate Editor 


him to an understanding of the object of 
experimentation, and that will help him 
to interpret experimental results. He re- 
quires discussion of the quantities to be 
measured, of the possible lines of attack on 
the problems and of the principle experi- 
mental techniques and attainable accu- 
racies. 

However, an abundance of detail is 
provided for the worker who plans to use 
the methods, and there are presented with- 
out boring the general reader. But the 
experimentalist, particularly the novice, 
will have difficulty in applying some of the 
methods. Expanded analysis and cata- 
loguing of the references, perhaps designat- 
ing the basic and important references by 
asterisks, would have helped to supple- 
ment the necessarily concise text. Also, 
more discussion of the details and tricks 
not readily found in references would have 
been greatly appreciated by the experi- 
menter. 

Little space is devoted to the orientation 
of the general reader. There is no intro- 
ductory chapter to acquaint the reader 
with the important properties to be meas- 
ured in high speed and combustion flows 
and with the circumstances of measure- 
ment. An introductory chapter might 
also have outlined the methods of measure- 
ment and compared them. It might have 
indicated possible accuracies and ranges of 
applicability. In addition to this general 
introduction, introductory chapters for 
the several divisions of the book would 
have been valuable. The first chapter of 
the book, is, indeed, an introduction to 
optical methods, supplying the necessary 
theoretical background for interferometric 
schlieren and shadow methods. A similar 


but less theoretical chapter for each of the 
other divisions would have been appropri- 
ate. All divisions would have benefited 
from discussion of the quantities to be 
measured, conditions of measurement and 
comparison of possible techniques. An 
introduction would also have provided the 


opportunity for relating one division to 
other portions of the book, showing rela- 
tions between measured quantities and 
comparing widely differing techniques. 

Both the experimentalist and the general 
reader should be pleased with the coverage 
provided by this book. Personal pref- 
erence might, in particular instances, dic- 
tate a different choice of subject matter 
but there can be no great quarrel about 
the selection that has been made. The 
present reviewer, for example, would like 
to have seen the section on pressure meas- 
urements expanded. He would also like 
to have seen more discussion of flow visu- 
alization techniques, which frequently 
provide the avenues of approach to under- 
standing the phenomena in new and com- 
plex fields. Space in the book might also 
have been found for methods of measure- 
ment of fluid properties such as humidity, 
composition, ete. 

These criticisms are minor and the au- 
thors are to be congratulated on the articles 
they have produced. All of the chapters 
areAcientifically sound and concisely pro- } ( 
vide the necessary information/ Several 
chapters are superb. This book isfrecom- 
mended both as an introduction to the 
experimental methods of high speed gas 
dynamics and combustion and as a refer- 
ence for an experimentalist working in 
these fields. / 


Thermodynamics and Statistical Me- 
chanics, by A. H. Wilson, Cambridge 
University Press, New York, 495 pp. 
$9.50. 


Reviewed by E. F. Obert 
Northwestern University 


There are so many books on thermo- 
dynamics and statistical mechanics that 
each new book should have a definite ob- 
jective and a definite place in the litera- 
ture. This new book by A. H. Wilson ful- 
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fills these prerequisites since it is directiy 
aimed at the theoretical physicist and pre- 
sents a detailed but compact survey of the 
subjects of thermodynamics and statis- 
tical mechanics and their interrelation- 
ships. The engineer in research or in 
teaching will find the book to be valuable 
since the physical significance of the 
various steps in analysis is emphasized 
while the mathematical detail is mini- 
mized. Thus the presentation stresses 
the salient points of thermodynamics as 
illustrated by the statistical mechanics, 
without the intervention. of detail which 
elementary texts must include for full 
understanding. 

The first three chapters would seem to 
belie the objectives since classical thermo- 
dynamics is covered through a consider- 
able range of elementary as well as ad- 
vanced topics. The concepts of tempera- 
ture, the first and second laws and equilib- 
rium are discussed briefly in the usual 
manner, with greater detail paid to the 
stability requirements of homogeneous 
substances. But starting with Chapter 
Four entitled, ‘“The Axiomatic Foundation 
of Thermodynamics,’ the diversity of 
coverage and the objectives of the author 
come alive. Here temperature and the 
first and second laws are adequately re- 
viewed in the axiomatic manner of Cara- 
theodory following the simplified analysis 
of Born. In Chapter Five “Statistical 
Mechanics,” the fundamental concepts of 


statistical mechanics, are reviewed briefly 
and then applied to obtain the concepts, 
first, of temperature and then of entropy, 
and other thermodynamic functions. The 
chapter is completed by sections on the 
perfect gas, the transition to classical 
mechanics, and the quantal statistical 
mechanics. 

Chapter Six, “Some Applications of 
Statistical Mechanics,’ reviews the spe- 
cific heats of gases, crystalline solids, Fermi- 
Dirac and Einstein-Bose systems as pre- 
requisite material for Chapter Seven, ‘“The 
Third Law of Thermodynamics.”’ The 
discussion of the third law is excellent (and 
illustrates why the author emphasized 
certain topics in previous chapters). 
Chapter Eight, ‘Imperfect Gases,” is also 
an excellent and complete job and dis- 
cusses the behavior of gases from the em- 
pirical through the thermodynamic as- 
pects with added discussion from the view- 
point of statistical mechanics. 

The diversity of coverage is shown by 
the topics in the remainder of the book: 
Equilibria, the properties of helium, elec- 
tric and magnetic phenomena, polar sub- 
stances, superconductivity, gas mixtures, 
chemical reactions, solutions and solids. 

Whether or not the book will appeal to 
the theoretical physicist must be answered 
by a member of that group; this reviewer 
is of the opinion that the book will be a 
valuable addition to the library of the en- 
gineer engaged in teaching or research 


since it is one of the relatively few books in 
thermodynamics which stresses the phys- 
ics, rather than the chemistry of the sub- 
ject, and supplies immediately the added 
information offered by the statistical 
mechanics. 


Thermal Stresses, by B. E. Gatewood, 
McGraw-Hill Book Company, Inc., 
New York, 1957, 247 pp. $7.50 


Reviewed by P. N. STEVENS 
Convair-Fort Worth 


This book is based on a course which the 
author gave to graduate aeronautical e1- 
gineering students at the Air Force Insti- 
tute of Technology. The author presents 
a comprehensive approach which is not 
restricted only to the classical elastic 
thermal stress problem. All phases of 
the structural design problem are investi- 
gated and includes areas such as the tem- 
perature distribution, the elastic and in- 
elastic applied and thermal stresses, allow- 
able stresses, and various effects of ele- 
vated temperatures. 

In Chapter 1, the author gives the ther- 
mal-stress equations, starting with simple 
cases and building up to the general three- 
dimensional equations. In the rest of the 
book, reference is made to Chapter 1 for 
the thermal-stress formulas. Chapter 2 
gives the temperature equations for heat 
conduction, heat radiation, heat convec- 
tion and aerodynamic heating, and ap- 
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plies them to obtain the heat-balance 
equations for an idealized structural ele- 
ment on a typical aircraft or missile struc- 
ture. Reference is made in the later 
chapters to the temperature equations in 
Chapter 2. 


In Chapter 3 the temperatures and 
thermal stresses for an idealized structural 
element consisting of a skin element at one 
uniform temperature joined to a stringer 
element at a different uniform tempera- 
ture are determined. In this case time is 
the only independent variable for the 
idealized element. With certain assump- 
tions concerning the aerodynamic heating 
of the skin element, the heat radiation to 
the stringer, and the thermal resistance of 
the joint between the skin and the stringer, 
analytical solutions for the temperature 
and thermal stress are obtained. The 
important parameters involving geome- 
try, material properties and time are de- 
termined, and curves of maximum stress 
against these parameters are constructed. 


In Chapter 4 the procedures of Chapter 
3 are extended to the case of one-dimen- 
sional temperature and _ thermal-stress 
variation. With time and the thickness 
dimension as independent variables, ana- 
lytical solutions are obtained for a thick 
plate or a beam or a multiweb box beam 
of an aircraft structure. Chapter 5 gives 
some of the simpler two-dimensional 
thermal-stress problems of particular in- 
terest, but the general two-dimensional 
problem is considered in Chapter 9. 


Chapter 6 is devoted to material prop- 
erties, allowable stresses under various 
loading conditions, and efficiencies of ma- 
terials at elevated temperatures. Some 
discussion of thermal fatigue and thermal 
shock of materials also is given. In 
Chapter 7 the stress-analysis problem of 
combining the applied stresses, thermal 
stresses and allowable stresses is consid- 
ered. Particular emphasis is put on the 
inelastic thermal stresses, the inelastic 
combination of applied and thermal 
stresses and the time at which the applied 
and thermal stresses are combined. In 
Chapter 8 the deflection and stiffness of the 
structure under elevated temperatures and 
thermal stresses are considered. Simple 
procedures are given to determine the ap- 
proximate deflection and stiffness of the 
structure and the resultant effects upon the 
thermal stresses and upon flutter. 


In Chapter 9 some cases of the general 
two-dimensional thermal-stress problem 
are considered. In particular the energy 
and complex-variable methods are dis- 
cussed and applied to several problems. 


It is the reviewer’s opinion that this 
book fulfills in part the general need for a 
comprehensive coverage of the general 
general stress problem. However, the 
basic information and problem-solution 


procedures presented/lean heavily toward ) 


the aircraft structures field and are not 
directly relatable to problem areas in the 
nuclear and other fields. / The basic dif- 
ference is in the heat generation source 
terms and in the materials of construc- 
tion. It is felt that additional emphasis 
should have been given to these problem 
areas because of the strong influence im- 
posed by the thermal stress problem on all 
nuclear reactor design considerations. 
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cast it for you. 

Actually, no such casting exists. 

It could. 

We have developed unique foundry 
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thousands of pounds. These castings 
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thick which must be porosity free. 

We are also casting solid and hol- 
low mandrels used to form solid pro- 
pellant. 

The outstanding contribution of our 
large aluminum castings is drastically 
reduced first and ultimate costs. 
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now. 
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New Patents 


driven helicopter (2,811,324). 
Ralph P. Alex, Fairfield, Conn., assignor 
to United Aircraft Corp. 

Gas turbine engine mounted on a gear 
box aft of the rotor shaft, its jet discharge 
directed rearwardly. The power shaft is 
directed forwardly toward the rotor shaft, 
Variable area nozzle (2,811,831). Fred- 
erick L. Geary and Wm. Taylor Jr., 
South Glastonbury, Conn., assignors to 
United Aircraft Corp. 

Flaps pivotally mounted on a ring 
around the outlet. Bell crank levers 
mounted on the outlet duct, one for each 
flap, are connected to actuate adjacent 
flaps for varying the outlet area. 
Combustor for turbine of aircraft starter 
(2,811,832). Cameron D. re and 
Robert Wiedis, Maywood, N. J., assignors 
to Bendix Aviation Corp. 

A compact combustor to provide gases 
to the turbine of a jet engine starter. A 
turbulator plate provides an annular in- 
clined passage for air moving in pocket- 
like louvres projecting into a bow]-shaped 
combustion chamber. 

Turbine cooling (2,811,833). Wilgus S. 
Broffitt, Indianapolis, Ind., assignor to 
General Motors oan 

Conduits conduct cooling air from outer 

air jackets to a chamber at the front of the 


turbine rotor. A ring adjacent the cham- 
ber has passages rotatable in and out cf 
registery with the conduits to control the 
delivery of air from the conduits to the 
chamber 
Turbine blade cooling systems (2,812,157). 
Wm. A. Turunen, Patrick W. O’Connell 
and De Owen Nichols Jr., Ann Arbor, 
Mich., assignors to General Motors Corp. 
Rotor blades each with a coolant pas- 
sage in its central portion sealed at its 
outer end. Radially extending passages 
in the turbine wheel communicate with 
passages in the blades. A duct from a 
condenser chamber conveys coolant from 
the turbine wheel passage to the chamber, 


et engine or gas turbine with electric 
ignition (2,811,676). W. B. Smits, Voor- 
x 4 Netherlands, assignor to Smitsvonk 

Low voltage surface discharge spark 
plug for each burner, connected in series. 
Condensers, sequentially connected with a 
source of electrical energy, fire atomized 
fuel associated with each plug. 


Process for preparing a liquid oxygen 
explosive (2,812,246). Samuel W. Mar- 
tin, Oak Park, IIl., assignor to Great 
Lakes Carbon Corp. 

Method comprises preparing a hard, 


Eprtors Nore: Patents listed above were selected from the Official Gazette of the U.S. 


Patent Office. 


Printed copies of Patents may be obtained from the Commissioner of 


Patents, Washington 25, D. C., at a cost of 25 cents each; design patents, 10 cents. 


When you join Telecomputing’s 
Engineering Services Division, 
you will be given full scope to 
allow you to grow. ..your talents 
will be used to the fullest... 
recognition and rewards will be 
yours as a matter of course. 
Engineering services is a mem- 
ber of an integrated five-com- 
pany missile systems corporation 
which designs and manufactures 
its own data-processing equip- 
ment. Our efforts are directed to- 
ward the ‘‘systems concept.’’ We 
are engaged in the reduction of 
large amounts of flight test data 
being generated by the daily mis- 
sile firings on the integrated Hol- 
loman-White Sands range. 


DATA PROCESSING SPECIALISTS! 


Get in now —at the beginning 
of the new era in missiles! 


Send resume to Director of Technical Personnel 


TELECOMPUTING CORPORATION 
Engineering Services Division 
BOX 447 » HOLLOMAN AIR FORCE BASE + NEW MEXICO 


ATTRACTIVE SALARIES 
PROFIT SHARING 


RELOCATION PAY 
ACCREDITED EDUCATION 
GROUP INSURANCE 


A NEW LIFE IN NEW MEXICO'S 
FABULOUS “LAND OF ENCHANTMENT" 


MOUNTAIN SKIING AND DESERT 
RESORTS WITHIN 30 MINUTES! 


A WONDERFUL PLACE 
TO MAKE YOUR HOME — 
GRAND COUNTRY TO RAISE KIDS! 


George F. McLaughlin, Contributor 


granular nonactivated absorbent carbon 
by the thermal expansion of finely divided 
bituminous coal suspended in liquid 
oxygen. The oxygen contains gas main- 
tained at a temperature about 1200 to 
1700 C, and an upheat rate higher than 
1000 C per sec. 


ack (2,812,- 
Graham B. 


Delayed opening parachute 
148). Gordon L. Fogal an 
Brown, Alexandria, Va. 

Mechanism for releasing a parachute 
from a pack in response to barometric 
pressure corresponding to a_ preselected 
altitude. Firing of an explosive charge, 

es by switches, causes separation 

he chute from the pack. 
Pulse jet engines with rearwardly opening 
air inlet (2,812,635). J. Le Foll and J. H. 
Bertin, Neuilly-sur-Seine, France, as- 
signors to SNECMA Co. 

Air inlet connected to one end of the 
combustion chamber and having an orifice 
facing toward the rear. An exhaust pipe 
is fed solely from the chamber. A pipe 
bend of 180 deg connects the beginning of 
the exhaust pipe with the other end of the 
chamber to form a continuous gas con- 
veying duct. 

Process and device for deflecting jets 
@ 812,636). M. Kadosch, F. G. Paris, 

H. Bertin and R. H. Marchal, Paris, 
—— assignors to SNECMA Co. 

Outwardly curved convex extension of a 

wall of a propulsive nozzle in which fluid 
is blown in a tangent direction. A valve 
controls the supply of pressure fluid to 
the blowing means. 
Fuel air ratio regulation for combustion 
systems (2,812,637). Homer M. Fox, 
Bartlesville, Okla., assignor to Phillips 
Petroleum Co. 

System of recovering low pressure and 

high pressure air streams from a combus- 
tion apparatus. Means for flowing any 
stoichiometric excess of air from the high 
pressure air stream to the low pressure air 
stream. 
Fuel system (2,812,715). Arnold H. 
Redding and Donald F. Winters, Folsom, 
Pa., assignors to Westinghouse Electric 
Corp. 

Fuel tank pump for providing a pre- 
determined head pressure in the supply 
conduit. A valve in the recirculation 
conduit is responsive to a fuel demand 
signal to control recirculation of fuel 
pumped in excess of the quantity de- 
manded by the engine. 

Gas turbine engines (2,812,897). L. Mc 
T. Cameron, Bristol, England, assignor 
to The Bristol Aeroplane Co. 

Ring of hollow inlet guide vanes in the 
air intake. Air from the compressor is 
led to the interior of the vanes, and con- 
duits lead heat into exchange relationship 
with the turbine. 
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+BY THE  JRTCAL 


DOUBLE | 
DIAGNOSIS 
TO PREVENT 


“HEART DISEASE” 


Two of the most important factors that affect jet engine life, efficiency, and safe 
operation are Exhaust Gas Temperature (EGT) and Engine Speed (RPM). Excess 
heat will reduce “bucket” life as much as 50% and low EGT materially reduces 
efficiency and thrust. Any of such conditions will make operation of the aircraft 
both costly and dangerous. The JETCAL Analyzer predetermines accuracy of the 
EGT and (interrelatedly) Tachometer systems and isolates errors if they exist. 


The JETCAL ANALYZES JET ENGINES 10 WAYS: 


1) The Jetcat Analyzer functionally tests 
EGT thermocouple circuit of a jet aircraft or 
pilotless aircraft missile for error without 
running the engine or disconnecting any wir- 
ing. GUARANTEED Accuracy is +4°C, at en- 
gine test temperature. 

2) Checks individual thermocouples ‘fon the 
bench” before placement in parallel harness. 
3) Checks thermocouples within the harness 
for continuity. 

4) Checks thermocouples and paralleling 
harness for accuracy. 

5) Checks resistance of the Exhaust Gas 
Temperature system. 

6) Checks insulation of the EGT circuit for 
shorts to ground and for shorts between leads. 
7) Checks EGT Indicators (in or out of the 
aircraft). 

8) Checks EGT system with engine removed 


from aircraft (in production line or overhaul 
shop). 


9) Reads jet engine speed while the engine is 
running with a guaranteed accuracy of +0.1% 
in the range of 0-110% RPM. Additionally, 
the TaKCAL circuit can be used to trouble 
shoot and isolate errors in the aircraft tachom- 
eter system. 


10) Jetcat Analyzer enables engine adjust- 
ment to proper relationship between engine 
temperature and engine M for maximum 
thrust and efficiency during engine run (Tab- 
bing or Micing). 

ALSO functionally checks aircraft Over- 
Heat Detectors and Wing Anti-Ice Systems 
(thermal switch and continuous wire) by 
using TEMPCAL Probes.. Rapid heat rise .. . 
3 minutes to 800°F! Fast cycling time of 
thermal switches . . . 4 to 5 complete cycles 
per minute for bench checking in production. 


BeH INSTRUMENT Co., INC. 
3479 West Vickery Blvd. * Fort Worth 7, Texas 


Sales-Engineering Offices: 


Tests EGT System Accuracy to 
+4°C at Test Temperature 


(functionally, without running the engine) 


Tests RPM Accuracy to 10 RPM 
in 10,000 RPM (=0.1%) 


The JETCAL is in worldwide use .. . by the 
U.S. Navy and Air Force as well as by major 
aircraft and engine manufacturers. Write, 
wire or phone for complete information. 


VALLEY STREAM, L.I., N. Y.: 108 So. Franklin, LO 1-9220 @ DAYTON, ©.: 209 Commercial Bidg.,BA 4-4563 ® COMPTON, CAL.: 105 N. Bradfield St., NE 6-8970 
ENGLAND: Bryans Aeroquipment Ltd. (Licensee), 15, Willow Lane, Mitcham, Surrey 
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Nuclear Propulsion 


Tabulation, Bibliography and Structure 
of Binary Intermetallic Compounds, I. 
Compounds of Lithium, Sodium, Potas- 
sium, and Rubidium, by S. G. Epstein, 
D. M. Bailey, R. L. Smythe, G. R. Kilp 
and J. F. Smith, Atomic Energy Comm., 
1SC-795, Sept. 1956, 31 pp. 

Thermal and Electrical Properties of 
Graphite Irradiated at Temperatures from 
100 to 425 K, by G. E. Deegan, Atomic 
Energy Comm., NAA-SR-1716, Dec. 1956, 
77 pp. 

The Determination of Power and Tem- 
pa Overshoot in the Reactor Flash 

oblem, by K. T. Spinny, Atomic Energy 
Comm., RNL-421, Nov. 1956, 16 pp. 

Studies of Flow Distribution in the Core 
of a Quarter Scale Flow Model of the PWR 
Reactor, by Herbert R. Hazard and John 
M. Allen, Atomic Energy Comm., BMI- 
1141, Oct. 1956, 46 pp. 


Design Evaluation of BER (Boiling 


EpiTor’s Note: Contributions from E. R. 
G. Eckert, J. P. Hartnett, T. F. Irvine Jr. and 
P. J. Schneider of the Heat Transfer Labora- 
tory, University of Minnesota, are gratefully 
acknowledged. 


Technical Literature Digest 


M. H. Smith, Associate Editor, and M. H. Fisher, Contributor 
The James Forrestal Research Center, Princeton University 


Experimental Reactor) in Regard to In- 
ternal Explosions, by F. B. Porzel, Atomic 
Energy Comm., ANL-5651, Jan. 1957, 107 
pp. 

A Review of the Reaction Kinetics of 
Deuterium and Tritium Compounds, V: 
Association—Addition and Elimination 
Reactions, by Lawrence M. Brown, Atomic 
Energy Commission, NBS-4674, June 1956, 
25 pp. 

A Review of the Reaction Kinetics of 
Deuterium and Triterium Compounds, VI: 
Solvolysis Reactions, by Lawrence M. 
Brown, Atomic Energy Comm., NBS-4712, 
July 1956, 18 pp. 

Fusion for Flight, Missiles and Rockets, 
vol. 2, June 1957, p. 65. 

Atomic Power Packages for Missile 
Launching Sites, Missiles and Rockets, vol. 
2, June 1957, pp. 66-67. 

Air Force Studies Ion Power, by Henry 
T. Simmons, Missiles and Rockets, vol. 2, 
June 1957, pp. 76. 

M/r Bibliography, Atomic Power for 
Aircraft and Missiles, Missiles and 
ae vol. 2, June 1957, pp. 136, 138, 
140. 


Ion Propulsion and Why, by Ernest 
Stuhlinger, Missiles and Rockets, vol. 2, 
June 1957, pp. 82, 84-85. 


The Ion Rocket, by John E. Naugle, 
Missiles and Rockets, vol. 2, June 1957, 
pp. 87-88. 

Fuels for the Atomic Rocket, by Alfred 
J. Zaehringer, Missiles and Rockets, vol. 2, 
June 1957, pp. 93-94. 

Nuclear Powered Rockets, by Rudolf H. 
Reichel, Missiles and Rockets, vol. 2, June 
1957, pp. 96-97. 

Analysis of the Interaction of Electro- 
magnetic Radiation with a Plasma in a 
Magnetic Field, by Laurence 8. Hall, An- 
drew L. Gardner and Osmund T. Fund- 
ingsland, Atomic Energy Comm., UCRL- 
4744, Sept. 1956, 34 pp. 

Present Status of the Problem of Rocket 
Propulsion by Nuclear Energy, by |’. 
Winterberg, Fusées et Recherche Aeron., 
vol. 2, no. 1, March 1957, pp. 9-22 (in 
French). 

Some ——_ of Nuclear Energy Pro- 
ulsion, by Thomas Reis and Paul Ageron, 

usées et Recherche Aeron, vol. 2, No. |, 
March 1957, pp. 23-33 (in French). 


Pinch Effect, by Louis C. Burkhardt, 
Robert E. Dunaway, Joseh W. Mather, 
James A. Phillips, George A. Sawyer, 
Thomas F. Stratton, Emory J. Stovall Jr. 
and James L. Tuck, J. Appl. Phys., vol. 
28, May 1957, pp. 519-521. 


Self-Aligning Bearings 


CHARACTERISTICS. 


High Ambient 
Temperature 
Magnetic Transducers 


translate mechanical motion 
into AC voltage 


ANALYSIS 
Stainless Steel 
Ball and Race 
Steel Ball and Race 


Bronze Race and 
Chrome Steel Ball 


Thousands in use. Backed by years of service life. Wide variety 
of Plain Types in bore sizes 3/16” to 6” Dia. Rod end types in 
similar size range with externally or internally threaded shanks. 
Our Engineers welcome an opportunity of studying individual 
requirements and prescribing a type or types which will serve 
under your demanding conditions. Southwest can design special 
types to fit individual specifications. As a result of thorough 
study of different operating conditions, various steel alloys 
have been used to meet specific needs. Write for revised Engi- 
neering Manual describing complete line. Dept. JP-58. 


SOUTHWEST PRODUCTS CO. 


1705 SO. MOUNTAIN AVE., MONROVIA, CALIFORNIA 


RECOMMENDED USE 


For types ating under high temper- 
ature (800.1200 Gores 


Chrome Alloy For types operati der high radial 
{ Foe type loads (3000-893,000 ibs... 


{ For types operating under normal loads 
with minimum friction requirements. 


Model 3030 HTAN $26.50 list 


Operate in temperatures 
up to 500° F. 


Internal Teflon construction and 
Teflon insulated wire. 


Available in miniature size: 
Model 3015 HT $16.50 list 
(1-9/32" long — 1/3 oz.) 
Write for new bulletin CMP 656 
and quantity prices! 


ELECTRO PRODUCTS LABORATORIES 
4501-J Ravenswood Ave., Chicago 40, Ill. 
Canada: Atlas Radio Ltd., Toronto 
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Guidance Systems and 
Components 


Functional Operation of Guidance and 
Control Components in Kay-1 (LARK) 
During Noise and Vibration Test, by E. D. 
Bush, Naval Air Missile Test Center, Tech. 
Rep. 13, April 1947, 44 pp. (Declassified 
from Confidential by authority by NAM- 
TC letter A9/Ser. 362, 7/20/56.) 

German Data on SG-66 Stabilizer, by 
E. L. Burnside, North American Aviation, 
Inc., Rep. AL202, June 1947, 55 pp. (De- 
classified from Confidential by authority 
of North American Aviation, Inc., letter 
54D-6233, citing WADC letter WCSG, 
8/6/54.) 

Performance of the German SG-66 
Stabilized Platform under Laboratory 
Conditions, by R. M. Benson, North 
American Aviation, Inc., Rep. AL357, 
June 1948, 24 pp. (Declassified from 
Confidential by authority of North 
American Aviation, Inc., letter 54D-6233, 
citing WADC letter WCSG, 8/6/54.) 

Autopilot System of the NATIV (RTV- 
A-3) Missile, by K. R. Jackson and W. D. 
Lyon, North American Aviation, Inc., Rep. 
AL894, Nov. 1949, 85 pp. (Declassified 
from Confidential by authority of North 
American Aviation, Inc., letter 54D-6233, 
citing WADC letter WCSG, 8/6/54.) 

Doppler Navigation, by William J. Tull, 
Aeron. Engng. Rev., vol. 16, May 1957, 
pp. 58-61. 

Radio Tracking, Orbit and Communica- 
tion for the Earth Satellite, by John PI. 
Hagen, Aeron. Engng. Rev., vol. 16, May 
1957, pp. 62-66. 

Self-Contained Navigation Systems, by 
Robert W. West, Aeron. Engng. Rev., vol. 
16, May 1957, pp. 96-98. 

Missile Guidance Techniques and 
Guided Missile Weapons System Con- 
siderations, by William C. Hodgson and 
John C. Ryon, Fusées et Recherche Aeron., 
vol. 2, no. 1, March 1957, pp. 53-63 (in 
English). 


Out-of-Print Copies 
Of JP Now Available 


Under an agreement with the 
AMERICAN Rocket Society, Kraus 
Reprint Corp., 16, E. 46th St., New 
York 17, N. Y., has made available 
complete back files of Jer Proput- 
sion, dating back to the first year of 
publication. 

Vols. 1-22 (1930-1952) are offered 
at $148.50, bound. In addition, Vols. 
1-2, 3-7, 8-10 and 11-13 are avail- 
able at $14 each, unbound; Vol. 14 
and 15 at $9 per volume, unbound; 
Vol. 16, 17, 18, 19, 20 and 21 at $14 
per volume, unbound; and Vol. 22 at 
$12, unbound. Single units or vol- 
umes can be supplied bound at $1.50 
additionally per unit or volume. Vols. 
23-26 are scheduled for reprint early 
this year. 

Kraus is also offering back copies of 
the Journal of the British Interplane- 
tary Society. Vols. 1-14 (1934-1955) 
are offered at $180, bound; Vols. 1-5 
in one volume (sold as a unit only) at 
$20, unbound; and Vols. 6-14 at $20 
per volume unbound. 


JANUARY 1958 


STEPPER 


COUNT DIAL INPUT COUNT SELECTOR FOLLOWER 
SELECTOR PINION T PINION Ls-2 
COUNT SELECTOR 
GEAR 
72:1 72:1 
INDEX CAM 
DRIVE GEAR- \ 
SM-130 
STEPPER 
MOTOR 
INDEX\, “5-3 P 
CAM 
<cam 
FOLLOWER 
CLUTCH 
CLUTCH SOLENOID 
H-{A] GROUND -28Vv. DC. 
(B} POWER +28V. DC. 
EXTERNAL START 
D) EXTERNAL RESET 
(E) PULSE INPUT 
OPERATION KNOB 
AN3I02E-14S-5P 
STOP & RESET—-=START CONNECTOR 


CSTART, STOP, & RESET) 


This Automatic Pulse Timer mounts in a 
standard 314” mounting. The initial usage 
of the Automatic Pulse Timer was for a difficult instrumentation prob- 
lem encountered on test aircraft—timing the pulses from a fuel flow 
transducer and thus determining specific fuel consumption. It success- 
fully replaced a complex and unreliable method. 

The Automatic Pulse Timer incorporates an uni-directional Stepper 
Motor along with complimentary gears, cams, solenoids, switches, an 
indicator light and—for an accurate independent time base—a stop 
watch. It is designed to visually record the lapsed time of an occurance 
of a specific number of electrical impulses. The Pulse Timer can count 
pre-selected quantity of 2 to 60 pulses, having a uniform or variable 
rate up to 25 pulses per second. 

In this application the combined accuracy of the fuel flow transmitter 
and the automatic pulse timer is better than 1%, and of this the timer 
contributes essentially no error. When the broad input requirements 
are available, the unit can be used for timing pulses regardless of the 
source from which they may originate. 


DETAILED OPERATIONAL SEQUENCE IS AVAILABLE UPON REQUEST. 


STEPPER MOTORS corporation 


; Subsidiary of California Eastern Aviation, Inc. 


7441 West Wilson Avenue @ Chicago 31, Illinois 
@ WEST COAST FACILITY 11879 W. FLORENCE AVE. . CULVER CITY, CALIF. 
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SENIOR 
RESEARCH 
SCIENTIST 


INTERIOR BALLISTICS 


A challenging position now exists 
at Aerojet-General Corporation 
for a man of broad, fundamental 
background capable of original 
experimental and theoretical 
studies. The complex interrelations 
of mathematics, physics, and 
chemistry as applied to the inte- 
rior ballistics of rockets offer chal- 
lenging problems to the creative 
research man. 


The application of unsteady-state 
theory to ignition of propellant 
and to the study of temperature- 
resistant materials is typical of 
the problems encountered. Com- 
bustion of solid propellants, the 
rheological behavior of solid pro- 
pellants in complex grain configu- 
ration, and the gas flow and heat 
transfer in nozzles are related 
important fields for research. 


Particularly important is the ability 
to formulate experiments for the 
verification of theoretical calcula- 
tions, devise the necessary equip- 
ment, conduct the experiments, 
and analyze the data using the 
most advanced techniques, 
including digital and ana- 
logue computers. 


For further details on this in- 
teresting position in an ideal 
Western climate, 


Dr. James H. Wiegand, 
write: Head, Interior Ballistics 
Research Dept. 


AEROJET-GENERAL 


CORPORATION 
P.0. Box 1947 


SACRAMENTO, CALIFORNIA 


Recent Advances in Cermets 
(Continued from page 15) 


tions in addition to their high melting 
and oxidation resisting properties. 


Bearing Materials 


It might be worthwhile mentioning 
that a use of true cermets at relatively 
low temperatures has lately become 
interesting: The possibility of using 
these materials, with or without binder, 
as non-lubricated bearings is being in- 
vestigated at temperatures between 800 
and 1500 F (5). At these temperatures, 
the ordinary steels or common bearing 
materials either become too soft or 
their frictional properties are no longer 
favorable. It is believed that cermets 
might have a very interesting and im- 
portant role to play in such applications. 
Not only does the high melting point 
indicate that they are unlikely to seize 
or gall at these temperatures, but the 
number of possible combinations of un- 
like materials is very large, especially if 
it is considered that both binder and 
hard phase can be varied. It would be 
desirable if not only the sleeve of the 
bearing could be produced from cermets 
but also the shaft; the latter should, 
however, be made of a different mate- 
rial. Problems of joining these mate- 
rials to the standard steels would still 
have to be solved. The low friction 
properties against steel of WC-TiC cut- 
ting materials (which are a kind of a 
cermet) are well known, and it is 
believed that the further development 
and testing of hard metals against dif- 
ferent hard metals will finally result in 
obtaining materials for high tempera- 
ture bearing applications. Very little 
work has so far been done on the deter- 
mination of frictional properties in the 
range of 1000 F, but the interest has 
been aroused, and it is believed that this 
field might bring about the first large- 
scale application for some of these 
materials. 


Research and Development Work 


Research and development work dur- 
ing the last few years has gone in two 
directions: The first one is an attempt 
to better understand the reasons why 
cermets behave as they do and to find, 
if possible, means to improve the impact 
resistance and ductility of some com- 
positions. The other development work 
has investigated new production meth- 
ods for cermet materials with the inten- 
tion to increase the sizes of the parts 
which can be manufactured and to make 
fabrication of sheet-like parts possible. 


Microstructure 


The fundamental research work has 
concentrated on the study of the micro- 


structure of a number of compositions 
and of the role the surface energy be- 
tween metal binder and hard material 
plays in ductility and high temperature 
properties. This has been done both in 
systems involving true oxide ceramics 
and others using metallic refractories, 
The wetting properties between ferrous 
metals and copper on one side and 
ceramics and hard metals on the other 
have been carefully studied (6), and it 
was found that a zero wetting angle is 
at least desirable, if not necessary, for a 
good cermet material. The original 
idea that the strength at high tempera- 
ture is due to a continuous skeleton of 
ceramic material is no longer considered 
correct (7), at least for the largest num- 
ber of cases. On the contrary, it was 
found that isolated particles uniformly 
dispersed in a matrix result in much 
more desirable properties (8). Impact 
behavior is improved without sacrifice 
of strength. Norton’s idea (9) that the 
strength of a cermet is essentially due to 
the peculiar state of stress in the binder 
material seems to be substantiated. 
The dependence of strength and other 
properties on the distribution of these 
isolated particles in the matrix, as well 
as on the size of the particles them- 
selves, has also helped to bring about a 
better understanding of these materials 
(10). Infiltration of a carbide skeleton 
with metal binder material has been 
used as a means to obtain a desirable 
microstructure. After infiltration and 
soaking of the material at high tempera- 
ture, the original skeleton is not pre- 
served, but the carbide particles become 
separated from each other and well 
distributed in the matrix (11). How- 
ever, it must be said that all these 
attempts (in spite of the progress they 
might have brought towards a better 
understanding) have not been able to 
produce a material which would be 
usable without fundamentally changing 
the design presently in use for engines or 
other applications. It is still too diffi- 
cult to work with essentially brittle ma- 
terials. A lot has to be learned about 
the interactions between binder and 
hard material, about the nature of the 
bond between the two constituents, and 
about the effect the microstructure has 
on the general properties. We are far 
from the time where we can prescribe a 
certain microstructure in order to im- 
part to a material the properties desir- 
able for a certain application. We know, 
however, that the degree to which the 
fundamental properties of both con- 
stituents become apparent in the be- 
havior of the final product depends to a 
large extent on the microstructure 
present. It seems now, even more so 
than some years ago, that the kind and 
nature of refractory hard metal or 
ceramic used in the cermet material is 
of lesser importance than the distribu- 
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tion of the particles within the matrix. 
How far grain size plays a role in the 
pure hard metals or ceramics is still not 
too certain. 

The materials most frequently used 
for all these investigations are the stand- 
ard cermets. No essentially different 
material has been developed or investi- 
gated during the last few years. This 
is due to the reasoning mentioned 
above. Titanium carbide-nickel cer- 
mets have been most widely studied and 
some work has been done on chromium 
carbide, titanium boride, zirconium bo- 
ride and some similar materials. In the 
case of mixtures of true ceramics with 
metals, alumina-base cermets have been 
studied (12) and a magnesium oxide- 
nickel material (13) has been developed. 

A peculiar kind of cermet was recently 
produced for applications in the nuclear 
energy field, especially for control rods. 
The material is used rather extensively 
and is called “Boral’’; it consists of 
boron carbide dispersed in an aluminum 
matrix. It can be rolled out into sheets. 
The same idea is applied to a stainless 
steel material containing Samarium or 
Europium or Gadolinium oxide (14). 


An excellent, uniform distribution of | § 


the neutron absorbing material within 
the steel matrix can be obtained. On 
the other hand, the use of pure borides 
and silicides of rare earth metals such as 
Gadolinium, Dysprosium or Europium 
was recently made possible (15). 


New Production Methods 


The development of new manufactur- 
ing procedures has borrowed a lot of ex- 
perience from the ceramic industry. 
Slip casting (16) has found applications 
in the production of larger parts, cruci- 
bles or containers from cermet mate- 
rials. The standard methods are 
applied with only slight modifications. 

To produce even larger pieces, or to 
impart the desirable cermet properties 
to very large areas, coating processes 
have been developed. The well-known 
method for spraying oxide materials can 
also be applied, with some modifications, 
to cermets. It is now possible to coat 
steels, graphite, refractories and other 
materials with borides, silicides and car- 
bides which melt below 5000 F. De- 
velopment work to extend the spraying 
methods to extremely high melting ma- 
terials is in progress. For the time 
being, coatings are applied as protection 
against attack by liquid metals (17), 
but applications as high temperature 
protection are under test. The coatings 
are usually porous which, however, has 
proved to be non-objectionable in many 
cases. Another offspring of this de- 


velopment is the production of large | 


pieces by spraying the cermet material 
against a thin sheet metal form (18). 
The material can be deposited to any 
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THESE ENGINES ARE 
NOW IN PRODUCTION 


AT CAE. 


... the 7T-37 Twin Jet Trainer 
with 
C.A.£. Turbine Power 


Air Corps flight training routine took 
a significant step forward recently, 
when the T-37 twin jet trainer entered 
its Phase VIII testing at Bainbridge 
Air Base, Georgia. Twenty hand- 
picked officers embarked on a course 
known as PROJECT PALM, with the 
two-way goal of training for them, 
and suitability testing for the plane. 
This new high-performance ship 
advances the jet phase of fliers’ 
training to an earlier stage in the 
training schedule, speeding the tran- 
sition from propeller-driven planes to 
jets, with gains in both safety and 
economy. Twin J69-T-9 turbines by 
C.A.E. provide the power. 
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GILMORE 


sy 
SIX COMPONENT 
THRUST 
MEASURING 
SYSTEM* 


Completes First Year 
of Missile Testing 


Installed at a West Coast missile testing 
facility, in August, 1956, this Gilmore weight 
and thrust measuring system indicates and 
records values of missile weight, propellant 
weight, thrust, gimballing motor forces, and 
thrust mis-clignment forces. Forces are trans- 
lated into. values represented by terms (a) 
weight and thrust, (b) pitch moment, (c) yaw 

t, (d) roll t, (e) XX axis side- 
load, (f) YY axis sideload. The system com- 
pensates for the effect of sideload with 
respect to wind and with respect to inter- 
action with moments. 


Computing equip t isolates all fixed com- 
ponents. A portable calibration unit, cali- 
brated by the National Bureau of Standards, 
is included in the installation. System accuracy 
is 0.1%. 

*Patents Pending 


For details on this and other Gilmore Systems 
Write Dept. J-1 


Instrumentation Systems 
For Industry and Science 


GILMORE INDUSTRIES ,ING. 


5713 Euclid Ave. @ Cleveland 3, Ohio 
West Coast Office: 5245 King St., Riverside, Calif. 


desired thickness, and a strong part is 
formed which has no dimensional limita- 
tions. The sheet metal form can later 
be removed by chemical means if de- 
sired, so that a pure cermet piece is 
obtained as an end-product. The sur- 
face of the pre-form is reproduced in 
every detail and a smooth surface can 
readily be obtained. Ribs or other 
stiffeners can be incorporated without 
difficulty. The method is still in its 
infancy, but it is expected to open new 
applications for cermet materials. Nose 
cones, leading edge parts and similar 
components which were formerly pro- 
duced from sheet metal seem to be 
naturals for this method. The strength 
of the material as produced by this 
“Spray-Forming” method should be 
satisfactory for most applications. If 
desired, however, an additional sintering 
treatment can be applied which im- 
proves the properties to some degree. 

Even though the developments of the 
last two years have not shown a con- 
siderable increase in the application of 
cermets, it seems to be clear that these 
materials are still of extreme interest for 
all high temperature applications, both 
in the aircraft and in the nuclear indus- 
tries. Designing for cermet materials is 
still very difficult, and a lot has to be 
learned to make good use of them. 
However, as the requirements for higher 
temperatures are increasing constantly, 
we will soon have reached a point where 
no other materials but cermets, hard 
metals and oxide ceramics will have a 
chance for proper performance. 
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